Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Slosh”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 523 records · Page 29

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

Numerical studies of surface tensions

Liquid-vapor (bubble) interface disturbances caused by various types of accelerations, including centrifugal, lateral and axial impulses, gravity gradient and g-jitter accelerations associated with spinning and slew motion in microgravity, are reviewed. Understanding of bubble deformations and fluctuations is important in the development of spacecraft orbital and attitude control techniques to secure its normal operation. This review discusses bubble deformations and oscillations driven by various forces in the microgravity environment. The corresponding bubble mass center fluctuations and slosh reaction forces and torques due to bubble deformations are also reviewed.

Hung, R. J.↗

Test-particle simulations in increasingly strong turbulence

Quasi-linear theory supposes that the energy in resonant fluctuations is small compared to that in the mean magnetic field. This is evident in the fact that the zeroth-order particle trajectories are helices about a mean field B(sub o) that is spatially uniform over many correlation lengths. However, in the solar wind it is often the case that the fluctuating part of the field is comparable in magnitude to the mean part. It is generally expected that quasi-linear theory remains viable for particles that are in resonance with a region of the fluctuation spectrum having only small energy density, but even so, care must be taken when comparing simulations to theoretical predictions. We have performed a series of test-particle simulations to explore the evolution of ion distributions in turbulent situations with varying levels of magnetic fluctuations. As delta-B/B(sub o) is increased the distinctions among absolute pitch angle (defined relative to B(sub o)), local pitch angle (defined relative to B(x)), and magnetic moment become important, some of them exhibiting periodic sloshing unrelated to the nonadiabatic processes of interest. Comparing and contrasting the various runs illustrates the phenomena that must be considered when the premise underlying quasi-linear theory are relaxed.

Pontius, D. H., Jr.↗

Expendable Launch Vehicle Studies

Analytical support studies of expendable launch vehicles concentrates on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane linearized dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed. The design of a robust LQR controller based on the reduced order system is accomplished using the parameter perturbation technique. The ELV modeling and analysis team has spent the past three years working on the theoretical development and application of sensitivity analysis to solve eigenvalue problems associated with structural dynamics. Specific application areas include stochastic vibrations, viscously damped vibrations and nonlinear dynamics. In stochastic linear vibrations, sensitivity analysis methods have been developed which determined the expectation and variance of the response for given probability density functions of various stochastic parameters. Examples worked include beam vibrations with up to six stochastic parameters; and the eigensensitivity results differ little from those obtained, with much effect, from Monte Carlo techniques. For viscously damped vibrations, eigensensitivity analysis has given excellent results for homogeneous beams, modeled by either quadratic or quartic eigenvalue equations. Specific applications include Kevin and Maxwell-type viscoelastic beams.

Bainum, P. M.↗

Dynamic modelling and simulation study for the Galileo spacecraft pulsed-mode spinup/400 N main engine burn/spindown maneuvers

Two Galileo dynamic models were developed to simulate the spinup, the main engine burn and the spindown maneuvers for the critical events of the Jupiter orbit insertion and perijove raise maneuver. The dynamic interaction between the spin thruster pulsing frequency, the science magnetometer boom flexible modes and the propellant slosh modes were studied. The prediction of safe orbit insertion was validated in flight. For the perijove raise maneuver, the simulation results indicated that propellant unporting would occur if the original spin thruster duty cycle was not updated. Based on further simulation results, a new duty cycle was selected to prevent propellant unporting. Subsequent Galileo flight data indicated that the perijove raise maneuver was successfully executed.

Ih, Che-Hang Charles↗

Analysis of Hydrodynamic (Landau) Instability in Liquid-Propellant Combustion at Normal and Reduced Gravity

The burning of liquid propellants is a fundamental combustion problem that is applicable to various types of propulsion and energetic systems. The deflagration process is often rather complex, with vaporization and pyrolysis occurring at the liquid/gas interface and distributed combustion occurring either in the gas phase or in a spray. Nonetheless, there are realistic limiting cases in which combustion may be approximated by an overall reaction at the liquid/gas interface. In one such limit, the gas flame occurs under near-breakaway conditions, exerting little thermal or hydrodynamic influence on the burning propellant. In another such limit, distributed combustion occurs in an intrusive regime, the reaction zone lying closer to the liquid/gas interface than the length scale of any disturbance of interest. Finally, the liquid propellant may simply undergo exothermic decomposition at the surface without any significant distributed combustion, such as appears to occur in some types of HydroxylAmmonium Nitrate (HAN)-based liquid propellants at low pressures. Such limiting models have recently been formulated,thereby significantly generalizing earlier classical models that were originally introduced to study the hydrodynamic stability of a reactive liquid/gas interface. In all of these investigations, gravity appears explicitly and plays a significant role, along with surface tension, viscosity, and, in the more recent models, certain reaction-rate parameters associated with the pressure and temperature sensitivities of the reaction itself. In particular, these parameters determine the stability of the deflagration with respect to not only classical hydrodynamic disturbances, but also with respect to reactive/diffusive influences as well. Indeed, the inverse Froude number, representing the ratio of buoyant to inertial forces, appears explicitly in all of these models, and consequently, in the dispersion relation that determines the neutral stability boundaries beyond which steady, planar burning is unstable to nonsteady, and/or nonplanar (cellular) modes of burning. These instabilities thus lead to a number of interesting phenomena, such as the sloshing type of waves that have been observed in mixtures of HAN and TriEthanolAmmonium Nitrate (TEAN) with water. Although the Froude number was treated as an O(1) quantity in these studies, the limit of small inverse Froude number corresponding to the microgravity regime is increasingly of interest and can be treated explicitly, leading to various limiting forms of the models, the neutral stability boundaries, and, ultimately, the evolution equations that govern the nonlinear dynamics of the propagating reaction front. In the present work, we formally exploit this limiting parameter regime to compare some of the features of hydrodynamic instability of liquid-propellant combustion at reduced gravity with the same phenomenon at normal gravity.

Margolis, Stephen B.↗

Dynamics of Superfluid Helium in Low-Gravity

This report summarizes the work performed under a contract entitled 'Dynamics of Superfluid Helium in Low Gravity'. This project performed verification tests, over a wide range of accelerations of two Computational Fluid Dynamics (CFD) codes of which one incorporates the two-fluid model of superfluid helium (SFHe). Helium was first liquefied in 1908 and not until the 1930s were the properties of helium below 2.2 K observed sufficiently to realize that it did not obey the ordinary physical laws of physics as applied to ordinary liquids. The term superfluidity became associated with these unique observations. The low temperature of SFHe and it's temperature unifonrmity have made it a significant cryogenic coolant for use in space applications in astronomical observations with infrared sensors and in low temperature physics. Superfluid helium has been used in instruments such as the Shuttle Infrared Astronomy Telescope (IRT), the Infrared Astronomy Satellite (IRAS), the Cosmic Background Observatory (COBE), and the Infrared Satellite Observatory (ISO). It is also used in the Space Infrared Telescope (SIRTF), Relativity Mission Satellite formally called Gravity Probe-B (GP-B), and the Test of the Equivalence Principle (STEP) presently under development. For GP-B and STEP, the use of SFHE is used to cool Superconducting Quantum Interference Detectors (SQUIDS) among other parts of the instruments. The Superfluid Helium On-Orbit Transfer (SHOOT) experiment flown in the Shuttle studied the behavior of SFHE. This experiment attempted to get low-gravity slosh data, however, the main emphasis was to study the low-gravity transfer of SFHE from tank to tank. These instruments carried tanks of SFHE of a few hundred liters to 2500 liters. The capability of modeling the behavior of SFHE is important to spacecraft control engineers who must design systems that can overcome disturbances created by the movement of the fluid. In addition instruments such as GP-B and STEP are very sensitive to quasi-steady changes in the mass distribution of the liquid. The CFD codes were used to model the fluid's dynamic motion. Tests in one-g were performed with the main emphasis on being able to compute the actual damping of the fluid. A series of flights on the NASA Lewis reduced gravity DC-9 aircraft were performed with the Jet Propulsion Laboratory (JPL) Low Temperature Flight Facility and a superfluid Test Cell. The data at approximately 0.04g, lg and 2g were used to determine if correct fundamental frequencies can be predicted based on the acceleration field. Tests in zero gravity were performed to evaluate zero gravity motion.

Frank, David J.↗

Propellant Management in Booster and Upper Stage Propulsion Systems

A summary review of some of the technical issues which surround the design of the propulsion systems for Booster and Upper Stage systems are presented. The work focuses on Propellant Geyser, Slosh, and Orientation. A brief description of the concern is given with graphics which help the reader to understand the physics of the situation. The most common solutions to these problems are given with there respective advantages and disadvantages.

Fisher, Mark F.↗

MAP stability, design, and analysis

The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The design and analysis of the MAP attitude control system (ACS) have been refined since work previously reported. The full spacecraft and instrument flexible model was developed in NASTRAN, and the resulting flexible modes were plotted and reduced with the Modal Significance Analysis Package (MSAP). The reduced-order model was used to perform the linear stability analysis for each control mode, the results of which are presented in this paper. Although MAP is going to a relatively disturbance-free Lissajous orbit around the Earth-Sun L(2) Lagrange point, a detailed disturbance-torque analysis is required because there are only a small number of opportunities for momentum unloading each year. Environmental torques, including solar pressure at L(2), aerodynamic and gravity gradient during phasing-loop orbits, were calculated and simulated. Thruster plume impingement torques that could affect the performance of the thruster modes were estimated and simulated, and a simple model of fuel slosh was derived to model its effect on the motion of the spacecraft. In addition, a thruster mode linear impulse controller was developed to meet the accuracy requirements of the phasing loop burns. A dynamic attitude error limiter was added to improve the performance of the ACS during large attitude slews. The result of this analysis is a stable ACS subsystem that meets all of the mission's requirements.

Ericsson-Jackson, A. J.↗

MAP Stability, Design and Analysis

The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The design and analysis of the MAP attitude control system (ACS) have been refined since work previously reported. The full spacecraft and instrument flexible model was developed in NASTRAN, and the resulting flexible modes were plotted and reduced with the Modal Significance Analysis Package (MSAP). The reduced-order model was used to perform the linear stability analysis for each control mode, the results of which are presented in this paper. Although MAP is going to a relatively disturbance-free Lissajous orbit around the Earth-Sun L2 Lagrange point, a detailed disturbance-torque analysis is required because there are only a small number of opportunities for momentum unloading each year. Environmental torques, including solar pressure at L2, and aerodynamic and gravity gradient during phasing-loop orbits, were calculated and simulated. A simple model of fuel slosh was derived to model its effect on the motion of the spacecraft. In addition, a thruster mode linear impulse controller was developed to meet the accuracy requirements of the phasing loop burns. A dynamic attitude error limiter was added to improve the performance of the ACS during large attitude slews. The result of this analysis is a stable ACS subsystem that meets all of the mission's requirements.

Ericsson -Jackson, A.J.↗

X-33 Attitude Control System Design for Ascent, Transition, and Entry Flight Regimes

The Vehicle Control Systems Team at Marshall Space Flight Center, Systems Dynamics Laboratory, Guidance and Control Systems Division is designing under a cooperative agreement with Lockheed Martin Skunkworks, the Ascent, Transition, and Entry flight attitude control system for the X-33 experimental vehicle. Ascent flight control begins at liftoff and ends at linear aerospike main engine cutoff (NECO) while Transition and Entry flight control begins at MECO and concludes at the terminal area energy management (TAEM) interface. TAEM occurs at approximately Mach 3.0. This task includes not only the design of the vehicle attitude control systems but also the development of requirements for attitude control system components and subsystems. The X-33 attitude control system design is challenged by a short design cycle, the design environment (Mach 0 to about Mach 15), and the X-33 incremental test philosophy. The X-33 design-to-launch cycle of less than 3 years requires a concurrent design approach while the test philosophy requires design adaptation to vehicle variations that are a function of Mach number and mission profile. The flight attitude control system must deal with the mixing of aerosurfaces, reaction control thrusters, and linear aerospike engine control effectors and handle parasitic effects such as vehicle flexibility and propellant sloshing from the uniquely shaped propellant tanks. The attitude control system design is, as usual, closely linked to many other subsystems and must deal with constraints and requirements from these subsystems.

Hall, Charles E.↗

Closed Loop Software Control of the MIDEX Power System

The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The design and analysis of the MAP attitude control system (ACS) have been refined since work previously reported. The full spacecraft and instrument flexible model was developed in NASTRAN, and the resulting flexible modes were plotted and reduced with the Modal Significance Analysis Package (MSAP). The reduced-order model was used to perform the linear stability analysis for each control mode, the results of which are presented in this paper. Although MAP is going to a relatively disturbance-free Lissajous orbit around the Earth-Sun L2 Lagrange point, a detailed disturbance-torque analysis is required because there are only a small number of opportunities for momentum unloading each year. Environmental torques, including solar pressure at L2, and aerodynamic and gravity gradient during phasing-loop orbits, were calculated and simulated. A simple model of fuel slosh was derived to model its effect on the motion of the spacecraft. In addition, a thruster mode linear impulse controller was developed to meet the accuracy requirements of the phasing loop burns. A dynamic attitude error limiter was added to improve the performance of the ACS during large attitude slews. The result of this analysis is a stable ACS subsystem that meets all of the mission's requirements.

Castell, Karen↗

Performance Evaluation Gravity Probe B Design

This final report documents the work done to develop a 6 degree-of-freedom simulation of the Lockheed Martin Gravity Probe B (GPB) Spacecraft. This simulation includes the effects of vehicle flexibility and propellant slosh. The simulation was used to investigate the control performance of the spacecraft when subjected to realistic on orbit disturbances.

Francis, Ronnie↗

Cryogenic Fuel Tank Draining Analysis Model

One of the technological challenges in designing advanced hypersonic aircraft and the next generation of spacecraft is developing reusable flight-weight cryogenic fuel tanks. As an aid in the design and analysis of these cryogenic tanks, a computational fluid dynamics (CFD) model has been developed specifically for the analysis of flow in a cryogenic fuel tank. This model employs the full set of Navier-Stokes equations, except that viscous dissipation is neglected in the energy equation. An explicit finite difference technique in two-dimensional generalized coordinates, approximated to second-order accuracy in both space and time is used. The stiffness resulting from the low Mach number is resolved by using artificial compressibility. The model simulates the transient, two-dimensional draining of a fuel tank cross section. To calculate the slosh wave dynamics the interface between the ullage gas and liquid fuel is modeled as a free surface. Then, experimental data for free convection inside a horizontal cylinder are compared with model results. Finally, cryogenic tank draining calculations are performed with three different wall heat fluxes to demonstrate the effect of wall heat flux on the internal tank flow field.

Greer, Donald↗

Multiple Core Galaxies

Nuclei of galaxies often show complicated density structures and perplexing kinematic signatures. In the past we have reported numerical experiments indicating a natural tendency for galaxies to show nuclei offset with respect to nearby isophotes and for the nucleus to have a radial velocity different from the galaxy's systemic velocity. Other experiments show normal mode oscillations in galaxies with large amplitudes. These oscillations do not damp appreciably over a Hubble time. The common thread running through all these is that galaxies often show evidence of ringing, bouncing, or sloshing around in unexpected ways, even though they have not been disturbed by any external event. Recent observational evidence shows yet another phenomenon indicating the dynamical complexity of central regions of galaxies: multiple cores (M31, Markarian 315 and 463 for example). These systems can hardly be static. We noted long-lived multiple core systems in galaxies in numerical experiments some years ago, and we have more recently followed up with a series of experiments on multiple core galaxies, starting with two cores. The relevant parameters are the energy in the orbiting clumps, their relative.masses, the (local) strength of the potential well representing the parent galaxy, and the number of cores. We have studied the dependence of the merger rates and the nature of the final merger product on these parameters. Individual cores survive much longer in stronger background potentials. Cores can survive for a substantial fraction of a Hubble time if they travel on reasonable orbits.

Miller, R.H.↗

Modeling the Restraint of Liquid Jets by Surface Tension in Microgravity

An axisymmetric phase field model is developed and used to model surface tension forces on liquid jets in microgravity. The previous work in this area is reviewed and a baseline drop tower experiment selected 'for model comparison. A mathematical model is developed which includes a free surface. a symmetric centerline and wall boundaries with given contact angles. The model is solved numerically with a compact fourth order stencil on a equally spaced axisymmetric grid. After grid convergence studies, a grid is selected and all drop tower tests modeled. Agreement was assessed by comparing predicted and measured free surface rise. Trend wise agreement is good but agreement in magnitude is only fair. Suspected sources of disagreement are suspected to be lack of a turbulence model and the existence of slosh baffles in the experiment which were not included in the model.

Chato, David J.↗

Development of a New Membrane Casting Apparatus for Studying Macrovoid Defects in Low-G

A new membrane-casting apparatus is developed for studying macrovoid defects in polymeric membranes made by the wet- and dry-casting process in low-gravity. Macrovoids are large (10-50 micron), open cavities interspersed among the smaller pores in the substructure under the gelled skin surface layer of the cast membrane. Although their occurrence is considered endemic to the wet- and dry-casting process since they can lead to compaction or skin rupture in the membrane process, recent studies suggest several useful applications such as transdermal and osmotic drug delivery systems, miniature bioreactors, etc. However, lack of knowledge about the macrovoid formation mechanism is an obstacle to further development of applications using them. An on-going debate is the role of the surface-tension-driven solutocapillary convection during macrovoid formation. The rapid growth of macrovoids within 1-5 seconds and the high polymer concentration in and near macrovoids make it difficult to explain the mechanism of macrovoid growth by diffusion alone, which is the widely accepted hypothesis proposed by Reuvers et al. The hypothesis advanced by our research group can explain this rapid growth via a mechanism that involves diffusion from the casting solution in the meta-stable region to the macrovoid enhanced by solutocapillary convection induced by the steep nonsolvent concentration gradient in the vicinity of the macrovoid. Since macrovoid growth is hypothesized to be the interplay of a solutocapillary-induced driving force counteracted by viscous drag and buoyancy, eliminate the latter provides a means for testing this hypothesis. Moreover, free convection mass transfer in the nonsolvent immersion bath used to cause phase-separation in membrane casting complicates developing a model for both the wet-casting process and macrovoid growth. The low-g environment minimizes gravitationally induced free convection thereby permitting a tractable solution to the ternary diffusion equations that characterize membrane formation. NASA's Parabolic Flight Research Aircraft provides a small window of low-g (approximately 25 s) that can be used to study macrovoid development in both wet- and dry-cast membranes if an appropriate casting apparatus is used. This casting apparatus should be able to cast the membrane in both low- and high-g in a manner so that essential one-dimensional mass transfer conditions are achieved to insure lateral uniformity in the membrane. The apparatus used in previous research on membrane casting in low-gravity was operated with the plunger driven mechanism. The spring-loaded plunger pushes the bottom block containing the polymer casting solution well directly under the absorbent chamber located in the upper stationary block. However, membranes made via this casting apparatus often displayed lateral nonuniformities that precluded obtaining quantitative information on the macrovoid growth process. Thus, it was necessary to determine the reason for these structural irregularities observed in the low-g casting apparatus. Both experimental as well as computer simulation studies of the low-g casting apparatus established that the impulsive action of the plunger caused the undesired structural nonuniformities. The simulation results showed that the width-to-depth aspect ratio of the shallow well that contains the casting solution in this apparatus was not an important factor in minimizing this problem. Even for a 40:1 (width : depth) aspect ratio, any convection induced by the horizontal motion of the interface of the casting solution will be damped out within 6.25x10(exp 4) seconds. However, the experimental studies revealed that the impulsive motion of the plunger caused a 'sloshing' of the casting solution that had to be eliminated. Therefore, the plungerdriven mechanism was changed to a cam-driven mechanism that did not cause any impulsive motion of the casting solution. Other refinements to this new membrane-casting apparatus include provision for removing the membranes from the casting wells in a less destructive manner. This was accomplished by using a slit geometry for the casting well that permitted disassembly for removal of the cast membrane. The materials used in the construction of this casting apparatus were chosen to insure wetting at the side walls and to maintain precise control of the thickness of the polymer solution in the casting well. An additional provision in this new casting apparatus is the ability to carry out both wet- as well as dry-casting. As such, this apparatus permitted the first studies of the wet-casting of polymeric membranes in low-g. Both wet- and dry-casting experiments on NASA's KC-135 research aircraft employing this new membrane-casting apparatus are scheduled in July 2002. The morphology of the resulting membranes will be characterized using an environmental scanning electron microscope (ESEM). The results of these low-g studies will be reported later.

Lee, Hanyong↗

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↗