Magnetically Induced Rotating Rayleigh-Taylor Instability

요약

We present a protocol for preparing a two-layer density-stratified liquid that can be spun-up into solid body rotation and subsequently induced into Rayleigh-Taylor instability by applying a gradient magnetic field.

초록

Classical techniques for investigating the Rayleigh-Taylor instability include using compressed gasses¹, rocketry² or linear electric motors³ to reverse the effective direction of gravity, and accelerate the lighter fluid toward the denser fluid. Other authors^{e.g. 4,5,6} have separated a gravitationally unstable stratification with a barrier that is removed to initiate the flow. However, the parabolic initial interface in the case of a rotating stratification imposes significant technical difficulties experimentally. We wish to be able to spin-up the stratification into solid-body rotation and only then initiate the flow in order to investigate the effects of rotation upon the Rayleigh-Taylor instability. The approach we have adopted here is to use the magnetic field of a superconducting magnet to manipulate the effective weight of the two liquids to initiate the flow. We create a gravitationally stable two-layer stratification using standard flotation techniques. The upper layer is less dense than the lower layer and so the system is Rayleigh-Taylor stable. This stratification is then spun-up until both layers are in solid-body rotation and a parabolic interface is observed. These experiments use fluids with low magnetic susceptibility, |χ| ~ 10^-6 - 10^-5, compared to a ferrofluids. The dominant effect of the magnetic field applies a body-force to each layer changing the effective weight. The upper layer is weakly paramagnetic while the lower layer is weakly diamagnetic. When the magnetic field is applied, the lower layer is repelled from the magnet while the upper layer is attracted towards the magnet. A Rayleigh-Taylor instability is achieved with application of a high gradient magnetic field. We further observed that increasing the dynamic viscosity of the fluid in each layer, increases the length-scale of the instability.

서문

A density stratified fluid system consisting of two layers can be arranged in a gravitational field in either a stable or an unstable configuration. If the dense heavy layer underlies the less dense, light layer then the system is stable: perturbations to the interface are stable, restored by gravity, and waves may be supported on the interface. If the heavy layer overlays the light layer then the system is unstable and perturbations to the interface grow. This fundamental fluid instability is the Rayleigh-Taylor instability^7,8. Exactly the same instability may be observed in non-rotating systems that are accelerated towards the heavier layer. Due to the fundamental nature of the instability it is observed in very many flows that also vary greatly in scale: from small-scale thin film phenomena⁹ to astrophysical scale features observed in, for example, the crab nebula¹⁰, where finger-like structures are observed, created by pulsar winds being accelerated through denser supernova remnants. It is an open question as to how the Rayleigh-Taylor instability can be controlled or influenced once the initial unstable density difference has been established at an interface. One possibility is to consider bulk rotation of the system. The purpose of the experiments is to investigate the effect of rotation on the system, and whether this may be a route to stabilization.

We consider a fluid system that consists of a two-layer gravitationally unstable stratification that is subject to steady rotation about an axis parallel to the direction of gravity. A perturbation to an unstable two-layer density stratification leads to baroclinic generation of vorticity, i.e., overturning, at the interface, tending to break-up any vertical structures. However, a rotating fluid is known to organize itself into coherent vertical structures aligned with the axis of rotation, so-called 'Taylor columns'¹¹. Hence the system under investigation undergoes competition between the stabilizing effect of the rotation, that is organizing the flow into vertical structures and preventing the two layers overturning, and the destabilizing effect of the denser fluid overlying the lighter fluid that generates an overturning motion at the interface. With increased rotation rate the ability of the fluid layers to move radially, with opposite sense to each other, in order to rearrange themselves into a more stable configuration, is increasingly inhibited by the Taylor-Proudman theorem^12,13: the radial movement is reduced and the observed structures that materialize as the instability develops are smaller in scale. Fig. 1 shows qualitatively the effect of the rotation on the eddies that form as the instability develops. In the left hand image there is no rotation and the flow is an approximation to classical non-rotating Rayleigh-Taylor instability. In the right hand image all experimental parameters are identical to the left hand image except that the system is being rotated about a vertical axis aligned with the center of the tank. It can be seen that the effect of the rotation is to reduce the size of the eddies that are formed. This, in turn, results in an instability that develops more slowly than the non-rotating counterpart.

The magnetic effects that modify the stress tensor in the fluid may be regarded as acting in the same way as a modified gravitational field. We are therefore able to create a gravitationally stable stratification and spin it up into solid body rotation. The magnetic body forces generated by imposing the gradient magnetic field then mimic the effect of modifying the gravitational field. This renders the interface unstable such that the fluid system behaves, to a good approximation, as a classical Rayleigh-Taylor instability under rotation. This approach has been previously attempted in two dimensions without rotation^14,15. For an applied gradient magnetic field with induced magnetic field B, the body force applied to a fluid of constant magnetic volume susceptibility χ is given by f = grad(χB²/µ₀), where B = |B| and µ₀ = 4π × 10^-7 N A^-2 is the magnetic permeability of free-space. We may therefore consider the magnet to manipulate the effective weight of each fluid layer, where the effective weight per unit volume of a fluid of density ρ in a gravitational field of strength g is given by ρg - χ (∂B²/∂z)/(2 µ₀).

프로토콜

NOTE: The experimental apparatus is shown schematically in Fig. 2. The main part of the apparatus consists of a rotating platform (300 mm × 300 mm) mounted on a copper cylinder (55 mm diameter) that descends under its own weight into the strong magnetic field of a superconducting magnet (1.8 T) with a room temperature vertical bore. The platform is made to rotate via an off-axis motor that turns a slip-bearing with a keyhole orifice. The copper cylinder is attached to a key-shaped drive shaft that simultaneously rotates, and descends once the holding-pin is removed.

1. Preparation of Non-standard Equipment

1. Flotation boat
   1. Make the size of the boat such that it fits comfortably within the experimental tank without touching the sides.
      NOTE: The flotation boat (see Fig. 3) consists of polystyrene walls and a sponge base.
   2. Protect the sponge with a layer of strong tissue paper.
      NOTE: The purpose of the tissue paper is to dissipate as much vertical momentum from the fluid poured into the boat as possible.

2. Preparation of Experiment

1. Preparation of liquid layers
   1. Allow distilled water to come up to laboratory temperature (22 ± 2 °C). Approximately 650 mL is required for each experimental realization.
      NOTE: Allowing the mixture to equilibrate prevents formation of bubbles in the experiment due to exsolving air.
   2. Separate the distilled water into equal volumes in two separate containers, A and B, which will be used to prepare liquid for the dense lower layer and light upper layer respectively.
   3. Ex-situ preparation of dense lower layer. To the contents of container A:
      1. Add NaCl to achieve a concentration of 0.43 mol NaCl per liter of water (approximately 25 g of NaCl per liter of water will be required);
      2. Add 0.33 g red and blue water-tracing dyes to the lower layer container (e.g., Cole-Parmer 00295-16 & -18);
      3. Add 0.1 g L^-1 fluorescein sodium.
         NOTE: The lower layer will be now be opaque in appearance and have a density of approximately 1012.9 ± 1.2 kg m^-3.
   4. Ex-situ preparation of light upper layer. To the contents of container B:
      1. Add MnCl₂ salt to achieve a concentration of 0.06 mol MnCl₂ per liter of water (approximately 12 g of MnCl₂ per liter of water).
         NOTE: The upper layer will be transparent in appearance and have a density of approximately 998.2 ± 0.5 kg m^-3.
   5. To vary the viscosity of the fluid layers, add glycerol C₃H₈O₃ in equal amounts to each layer until the desired viscosity is attained. Typical viscosities lie in the range 1.00 × 10^-3 - 21.00 × 10^-3 Pa s. The viscosity of each layer is the same.
      NOTE: The mixtures may be safely stored in their separate containers until required.
   6. Ex-situ preparation of density stratification.
      1. Add 300 mL of the contents of container A to the cylindrical inner tank (see Fig. 2).
      2. Immerse the flotation boat's sponge in fluid from container B.
         NOTE: After (2.1.6.2) the procedure is time sensitive, so do not carry out any further steps until all the magnet and the lighting, recording and mechanical mechanisms are ready.
      3. Lift the flotation boat out of the container B and, when it has stopped dripping, carefully place the flotation boat on top of the layer of dense fluid in the inner cylindrical tank.
      4. Begin to add light-layer fluid from container B to the flotation boat at a flow rate of 3 mL/min. Gradually increase this flow rate as the flotation boat lifts away from the interface between the two layers. Maintain a slow enough flow rate that the interface is not disturbed by the increased momentum of the fluid flow, but fast enough that this process takes no more than 20 min. Keep filling until the upper layer contains 320 mL of fluid.
         NOTE: The lower layer will be at a depth of approximately 33 mm, and the upper layer will be at a depth of approximately 39 mm.
      5. Carefully lower the lucite lid into the upper layer such that the layer depths of each layer are equal. Allow fluid and air to flow through the bleed holes, ensuring that no air is trapped beneath. Observe a layer (approx. 6 mm) of clear light layer liquid on top of the lucite lid.
         NOTE: If the process has been successful there will be two layers of liquid of equal depth with a sharp interface between them. The thickness of the diffusion layer at the interface will be less than 2 mm at this stage.
   7. Fill the outer tank with clear distilled water to a height 6 mm above the lucite lid of the inner tank. Upon observing square-on there will be no curvature-induced parallax resulting from the inner cylindrical tank.
      NOTE: Since the liquids in each layer are continuously diffusing across the interface at this point, proceed immediately to the following steps.
2. Spin-up of the stratification
   1. Place the experimental tank on the platform.
   2. Position the arrangement with the copper cylinder in the bore of the magnet, the drive shaft through the keyhole orifice in the track and the holding pin in position. Ensure that the tank is far away (60 cm) from the magnet such that the magnetic forces on the liquids are negligible at this position.
      NOTE: Carrying the experimental tank containing the stratification presents few difficulties; long, low amplitude, sloshing waves set up by walking with the tank will decay away, having negligible effect on the quality of the interface achieved when floating the upper layer on.
   3. Turn on the motor, increasing the rate of rotation at 0.002 rad s^-2, spinning-up the fluid to the desired rotation rate. For the rotation rates in ¹⁶ the spin-up time was of the order 20 min - 60 min.
      NOTE: The fastest rotation rate used was 13.2 rad s^-1.

3. Execution of Experiment

1. Ensure that the magnet is indicating a field strength of 1.2 T, and that at the height at which the instability is initiated the field gradient is (grad B²)/2 = -14.3 T² m^-1, where B is the magnetic induction.
2. Ensure that the video camera is arranged such that when the drive shaft is in its lowest position either the side view of the experiment is in focus, or a plan view is in focus through a mirror placed above the experiment.
3. Ensure the ambient lighting is at the correct levels, such that none of the image captured by the camera is saturated, but that the full response is used (grayscale intensities in the range 0-255).
4. Begin video recording (240 fps). Use a remote control to prevent moving the camera while operating the record function.
5. Remove the holding pin, allowing the tank to descend, while rotating, into the magnetic field.

4. Reset Experiment

1. Reset experimental rig
   1. Use the remote control to stop the video recording.
   2. Save the movie file to disk.
   3. By hand, lower the voltage to the motor so that it slows to a standstill. Perform this gradually so as to prevent spillages.
   4. Remove experimental arrangement from magnet.
   5. Dispose of the mixed liquid layers appropriately (see Manganese Chloride Tetrahydrate MSDS).
   6. Rinse the tank with water (it does not need to be distilled), until all traces of salts have been washed away. Avoid direct skin contact with liquids.
   7. Dry the tank carefully with tissue paper to ensure that no residue is left that may contaminate subsequent experiments.

5. Image Processing

1. Extract the individual images from each movie frame and save in lossless .png format. Mask out any unwanted areas of each frame, for example the platform or copper cylinder.
2. Calculate the two-dimensional auto-correlation function¹⁶ of each image frame for 2 s after initiation of the instability using a discrete Fast Fourier Transform. Record the minimum, mean, and maximum value of the observed wavelength for the rotation rate of the experiment and the viscosity of the fluid layers.

결과

Fig. 4 shows the development of the Rayleigh-Taylor instability at the interface between the two fluids, for four different rotation rates: Ω = 1.89 rad s^-1 (top row), Ω = 3.32 rad s^-1, Ω = 4.68 rad s^-1, and Ω = 8.74 rad s^-1 (bottom row). The interface is shown evolving in time from t = 0 s (left hand column) with increments of 0.5 s to t = 3.0 s (right hand column). The right hand column therefore r...

토론

There are two critical steps within the protocol. The first is 2.1.6.4. If the light layer is floated on the dense layer too rapidly then irreversible mixing of the two miscible fluid layers takes place. It is essential that this is avoided and that a sharp (<2 mm) interface between the two layers is achieved. The second critical step is 3.1.5. If the experiment is released toward the magnet without being fully spun-up into solid body rotation or without the visualization and image capture apparatus in position and o...

자료

Name	Company	Catalog Number	Comments
Blue water tracing dye	Cole-Parmer	00295-18
Red water tracing dye	Cole-Parmer	00295-16
Sodium Chloride			>99% purity
Manganese Chloride Tetrahydrate			See MSDS
Fluorescein sodium salt
Magnet	Cryogenic Ltd. London