SEARCH · Engineering Papers

Results for “hyperbolic systems”

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.

396 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Barth, Timothy↗

NASA NTRS1975

Hyperbolic cosmic dust - Its origin and its astrophysical significance

The heliocentric radial distribution of the flux of hyperbolic cosmic dust particles, as measured by the Pioneer 8 and 9 spacecraft, is closely related to the radial variation of the spatial density of source or parent meteoroids. Within the limits of the experimental and theoretical uncertainties the spatial density of parent meteoroids, as deduced from the hyperbolic cosmic dust data, is found to be increasing with increasing heliocentric distance in the neighborhood of one AU. Other recent experimental evidence confirms this result. The new results also suggest that the ratio of the areal density of submicron-sized craters to the areal density of millimeter-sized craters will be less on the north-south faces of lunar rocks than on the east-west faces of the same rocks. The change in ratio is not as large as previously thought. It is noted that the solar system is not presently contributing significant amounts of dust to the interstellar medium though it may once have done so.

Zook, H. A.↗

NASA NTRS1981

Trends in unmanned planetary entry systems

Entry systems used in unmanned planetary exploration are discussed in terms of future missions and more complex demands made upon the design configurations of orbiters and atmospheric space probes. The presently on-going Galileo project involving a Jupiter orbiter and a probe of the Jovian atmosphere is mentioned, and the difficult entry requirements at a velocity of 48 km/sec and an angle of minus seven to minus 10 degrees into an atmosphere composed largely of hydrogen, are described. Re-entry parameters determined by the respective entry environment, type of entry approach (hyperbolic or elliptical orbit) entry velocity, and guidance technology are given attention. Aerobraking, involving relatively slow circularization of an initial elliptical orbit by repeated passes through the atmosphere, is investigated. Aerocapture, by which vehicles using controlled flight paths through the atmosphere, can go directly from a hyperbolic flyby trajectory to low circular orbit is described. Major attractions of aerobraking and aerocapture are analyzed in terms of modest technology needs and relatively small impact on spacecraft design.

French, J. R.↗

NASA NTRS1997

The effects of atmospheric [CO2] on carbon isotope fractionation and magnesium incorporation into biogenic marine calcite

The influences of atmospheric carbon dioxide on the fractionation of carbon isotopes and the magnesium incorporation into biogenic marine calcite were investigated using samples of the calcareous alga Amphiroa and benthic foraminifer Sorites grown in the Biosphere 2 Ocean system under variable atmospheric CO2 concentrations (approximately 500 to 1200 ppm). Carbon isotope fractionation was studied in both the organic matter and the skeletal carbonate. Magnesium analysis was to be performed on the carbonate removed during decalcification. These data have not been collected due to technical problems. Carbon isotope data from Amphiroa yields a linear relation between [CO2] and Delta(sup 13)C(sub Corg)values suggesting that the fractionation of carbon isotopes during photosynthesis is positively correlated with atmospheric [CO2]. [CO2] and Delta(sup 13)C(sub Corg) values for Sorites produce a relation that is best described by a hyperbolic function where Delta(sup 13)C(sub Corg) values increase between 300 and 700 ppm and decrease from 700 to 1200 ppm. Further investigation of this relation and Sorites physiology is needed.

Vieira, Veronica↗

NASA NTRS2008

Titan Flagship Mission 3-Degree-of-Freedom Simulation Analysis

A NASA flagship mission to Titan, the largest moon of Saturn and the only moon in the solar system with a significant atmosphere, has been designed that uses three separate spacecraft, each requiring significant interaction with the atmosphere. The first vehicle is a Titan lander for lower-atmosphere and surface science. The second is an aerial vehicle for aerial science at approximately 10 km altitude with an expected lifetime of one year. This spacecraft will use the natural winds of Titan to cover a large area over its lifetime. The third vehicle is a Titan orbiter that will interact with the atmosphere in two respects. The first atmospheric interaction is the orbital insertion maneuver that will be accomplished using aerocapture, during which time the hyperbolic approach of 6.5 km/s will be reduced to 1.6 km/s over 41 minutes with an exit periapsis altitude of 130 km. The second atmospheric interaction occurs after a propulsive maneuver has raised the periapsis after aerocapture to 1170 km, where the atmosphere will be sampled over several months. This is the first phase of aerosampling that covers southern latitudes. After a 3.3-year circular science phase at an altitude of 1700 km, a second phase of additional aerosampling is performed sampling northern latitudes. The atmospheric trajectory analysis for these three spacecraft will be discussed throughout this paper.

Prince, Jill L.↗

NASA NTRS2019

Aerodynamic and Aerothermal Simulations of Mars Concept Vehicles Using Overset DPLR

NASA's Descent System Studies (DSS) Program is studying various concept vehicles to enable landing of heavy payloads on the surface of Mars. While it is desirable to run high-fidelity CFD simulations to accurately assess the aerodynamic and aerothermal effects of various design changes during EDL, it is usually difficult to quickly generate high-quality grids suitable for such analyses. One approach to address this bottleneck in mesh generation is through the use oversetting grids. Although the overset approach is efficient and powerful in solving partial differential equations on complex geometries, new users often find it challenging to apply overset concepts for their simulations. For example, generating hyperbolic grids with sufficient overlap; priority in hole-cutting on multiple overlapping grids; and fixes to assemble overlapping viscous grids at the body surface. The objective of this presentation is to introduce a simple process that combines the advantages of near-body, point-matched, structured grids with oversetting background grids suitable for grid alignment. This approach allows for grids that can be sequenced, reclustering of mesh spacing at the wall, and grid alignment with the bow shock. The current methodology is tested on a Mid-L/D configuration using the overset DPLR code.

Overset Grids↗

NASA NTRS1973

The starting transient of solid propellant rocket motors with high internal gas velocities

A comprehensive analytical model which considers time and space development of the flow field in solid propellant rocket motors with high volumetric loading density is described. The gas dynamics in the motor chamber is governed by a set of hyperbolic partial differential equations, that are coupled with the ignition and flame spreading events, and with the axial variation of mass addition. The flame spreading rate is calculated by successive heating-to-ignition along the propellant surface. Experimental diagnostic studies have been performed with a rectangular window motor (50 cm grain length, 5 cm burning perimeter and 1 cm hydraulic port diameter), using a controllable head-end gaseous igniter. Tests were conducted with AP composite propellant at port-to-throat area ratios of 2.0, 1.5, 1.2, and 1.06, and head-end pressures from 35 to 70 atm. Calculated pressure transients and flame spreading rates are in very good agreement with those measured in the experimental system.

Peretz, A.↗

NASA NTRS1995

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Yee, H. C.↗

NASA NTRS1983

Aerocapture and aeroassisted orbital transfer - A high performance leverage space technology

Single pass aerocapture and aeroassisted orbital transfer designs are discussed for modifying spacecraft orbits using planetary atmospheres. Aerodynamic drag can be used to shift a hyperbolic flyby trajectory to a planetary orbit by having the spacecraft altering its in-plane L/D ratio to maintain constant drag until the desired velocity is attained. The spacecraft would then move out of the atmosphere, sometimes with a small propulsive assist. Sample missions to Uranus, Jupiter, and Titan are mentioned, with mission scenarios commencing with a Shuttle launch. Biconic-shaped spacecraft are cited as having the best L/D for insertion accuracy, a low hypersonic ballistic coefficient, and volumeric efficiency. Aeroassisted orbital transfer involves atmospheric braking in order to shed velocity on a return trip from higher orbits. No decision has yet been made as to the spacecraft shape, but calculations indicate that aeroassist braking before return to LEO for rendezvous with the Shuttle has the potential for the lowest cost/mass to GEO of any delivery system currently under consideration.

Cruz, M. I.↗

NASA NTRS2003

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗

NASA NTRSDate not supplied

Rocket Lab Venus - Enabling Low-Cost Interplanetary Missions

Rocket Lab’s mission to Venus [1], launching in January 2025, aims to demonstrate that small launch vehicles, such as Electron, and high-energy small spacecraft, such as Photon, can enable a new paradigm of regular, low-cost interplanetary missions. The Electron launch vehicle will deliver the Photon spacecraft carrying the entry probe into low earth orbit, at which point Photon will take over and put both itself and the attached probe on an Earth-Venus transfer trajectory, with assistance from a Lunar flyby. The primary science objective of the mission is to sample the Venusian cloud layers between atitudes of 45 and 60 km above the surface with an autofluorescing nephelometer. The ~1kg scientific instrument will be carried by a small ~17 kg direct entry probe which will separate from the Photon spacecraft approximately 30 minutes prior to atmospheric entry. The probe outer mould line is a scaled down version of the Deep Space 2 probes [2]. NASA’s newly developed 3D Woven Carbon Phenolic, a derivative of the Heatshield for Extreme Entry Environment Technology (HEEET) material [3], is an enabling technology for this mission, and it will comprise the probe forebody heat shield.The probe will communicate directly back to the Earth during its descent via a small on-board radio frequency (RF) antenna. The hyperbolic Venus encounter will see the probe enter the Venusia atmosphere at approximately 11 km/s with an entry flight path angle (EFPA) of about 10°. This relatively low entry angle was selected to increase the duration of the descent phase, reduce the peak deceleration, and reduce the sensitivity of the nominal entry environments to navigational errors. For reference, the Pioneer Venus probes entered at EFPAs varying from about 25° to 70°, resulting in peak decelerations between around 220 and 450 g’s and peak stagnation point heat fluxes ranging from about 40 to 75 MW/m2 [4]. In comparison, the Rocket Lab Venus probe will experience a lower preak deceleration of around 60 g’s and a lower peak stagnation point heat flux of about 15 MW/m2, at the cost of a greater total heat load and correspondingly higher thermal protection system (TPS) mass fraction. This paper recounts the development history of the TPS for the Rocket Lab Venus probe, with a focus on the unique requirements of the mission. In particular, the effect on the TPS design of the mission’s tight budget, short development timeline, slim mass and RF link margins, and mild but long duration entry heating. Rocket Lab’s initial investigations into a ‘hybrid’ TPS [5] based on a scaled version of the Pioneer Venus full-density carbon phenolic forebody shield, and the abandonment of this approach in favour of NASA’s 3D Woven Carbon Phenolic is presented. Additionally, the leveraging of previous testing under the HEEET project at higher conditions and the associated TPS margining strategy that was adopted to avoid the need for extensive mission-specific arc jet testing is discussed. Finally, this paper concludes with a brief consideration of the feasibility of a future low-cost aerocapture technology demonstrator mission based on the 2025 Rocket Lab Venus mission and entry probe.

Lyle Campbell↗

NASA NTRS1993

User's manual for the HYPGEN hyperbolic grid generator and the HGUI graphical user interface

The HYPGEN program is used to generate a 3-D volume grid over a user-supplied single-block surface grid. This is accomplished by solving the 3-D hyperbolic grid generation equations consisting of two orthogonality relations and one cell volume constraint. In this user manual, the required input files and parameters and output files are described. Guidelines on how to select the input parameters are given. Illustrated examples are provided showing a variety of topologies and geometries that can be treated. HYPGEN can be used in stand-alone mode as a batch program or it can be called from within a graphical user interface HGUI that runs on Silicon Graphics workstations. This user manual provides a description of the menus, buttons, sliders, and typein fields in HGUI for users to enter the parameters needed to run HYPGEN. Instructions are given on how to configure the interface to allow HYPGEN to run either locally or on a faster remote machine through the use of shell scripts on UNIX operating systems. The volume grid generated is copied back to the local machine for visualization using a built-in hook to PLOT3D.

Chan, William M.↗

NASA NTRS1990

An investigation of design optimization using a 2-D viscous flow code with multigrid

Computational fluid dynamics (CFD) codes have advanced to the point where they are effective analytical tools for solving flow fields around complex geometries. There is also a need for their use as a design tool to find optimum aerodynamic shapes. In the area of design, however, a difficulty arises due to the large amount of computer resources required by these codes. It is desired to streamline the design process so that a large number of design options and constraints can be investigated without overloading the system. There are several techniques which have been proposed to help streamline the design process. The feasibility of one of these techniques is investigated. The technique under consideration is the interaction of the geometry change with the flow calculation. The problem of finding the value of camber which maximizes the ratio of lift over drag for a particular airfoil is considered. In order to test out this technique, a particular optimization problem was tried. A NACA 0012 airfoil was considered at free stream Mach number of 0.5 with a zero angle of attack. Camber was added to the mean line of the airfoil. The goal was to find the value of camber for which the ratio of lift over drag is a maximum. The flow code used was FLOMGE which is a two dimensional viscous flow solver which uses multigrid to speed up convergence. A hyperbolic grid generation program was used to construct the grid for each value of camber.

Doria, Michael L.↗

NASA NTRS1994

System for locating the sources of wideband dE/dt from lightning

A system has been developed to measure wideband electic field derivatives (dE/dt) at five ground stations in a 15 km x 15 km network at Kennedy Space Center. Individual station responses are normalized using digital filters. Pulse-timing resolution is improved to much less than 50-ns sample interval by interpolation using packing in the frequency domain. A time tag for each pulse is defined as the mean of the times of the rising-edge half peak, peak, and falling-edge half peak. The standard deviation in these times defines the timing error and is shown to be a function of noise and bandwidth rather than digitization rate. Each of the four unknowns for a pulse source location (x,y,z) and time of occurrence (t) is found from the five time-tag measurements using different weightings for all five combinations of the four-station hyperbolic equations. Weighting factors and errors in x,y,z and t are estimated using error propagation techniques.

Thomson, E. M.↗

NASA NTRS2016

Optimizing Mars Sphere of Influence Maneuvers for NASA's Evolvable Mars Campaign

NASA's Human Spaceflight Architecture Team is refining human exploration architectures that will extend human presence to the Martian surface. For both Mars orbital and surface missions, NASA's Evolvable Mars Campaign assumes that cargo and crew can be delivered repeatedly to the same destination. Up to this point, interplanetary trajectories have been optimized to minimize the total propulsive requirements of the in-space transportation systems, while the pre-deployed assets and surface systems are optimized to minimize their respective propulsive requirements separate from the in-space transportation system. There is a need to investigate the coupled problem of optimizing the interplanetary trajectory and optimizing the maneuvers within Mars's sphere of influence. This paper provides a description of the ongoing method development, analysis and initial results of the effort to resolve the discontinuity between the interplanetary trajectory and the Mars sphere of influence trajectories. Assessment of Phobos and Deimos orbital missions shows the in-space transportation and crew taxi allocations are adequate for missions in the 2030s. Because the surface site has yet to be selected, the transportation elements must be sized to provide enough capability to provide surface access to all landing sites under consideration. Analysis shows access to sites from elliptical parking orbits with a lander that is designed for sub-periapsis landing location is either infeasible or requires expensive orbital maneuvers for many latitude ranges. In this case the locus of potential arrival perigee vectors identifies the potential maximum north or south latitudes accessible. Higher arrival velocities can decrease reorientation costs and increase landing site availability. Utilizing hyperbolic arrival and departure vectors in the optimization scheme will increase transportation site accessibility and provide more optimal solutions.

Merrill, Raymond G.↗

NASA NTRS1995

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - 2

The global asymptotic nonlinear behavior of 1 1 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODES) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDES.

Yee, H. C.↗

NASA NTRS1995

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - Part 2

The global asymptotic nonlinear behavior of 11 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDEs.

Yee, H. C.↗

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↗