R-T instability model of magnetic fluid and its numerical simulations

Abstract

The Rayleigh-Taylor(R-T) instability of ferrofluid has been the subject of recent research, because of its implications on the stability of stellar. By neglecting the viscosity and rotation of magnetic fluid, and assuming that the magnetic particles are irrotational and temperature insensitive, we obtain a simplified R-T instability model of magnetic fluid. For the interface tracing, we use five-order weighted essentially non-oscillatory(WENO) scheme to spatial direction and three-order TVD R-K method to time direction on the uniform mesh, respectively. If the direction of the external magnetic field is the same as that of gravity, the velocities of the interface will be increased. But if the direction of the external magnetic field is in opposition to the direction of gravity, the velocities of the interface will be decreased. When the direction of the external magnetic field is perpendicular to the direction of gravity, the symmetry of the interface will be destroyed. Because of the action which is produced by perpendicular external magnetic field, there are other bubbles at the boudaries which parallel the direction of gravity. When we increase the magnetic susceptibility of the magnetic fluids, the effects of external magnetic fields will be more distinct for the interface tracing.

Keywords

magnetic fluid / Rayleigh-Taylor instability / weighted essentially non-oscillatory scheme / TVD R-K method

References

[1] YamaguchiH.Engineering fluid mechanics [M], 2008, Netherlands, Springer

[2] RosensweigR. E.Ferrohydrodynamics [M], 1985, New York, Combridge University Press

[3] ShliomisM. I.. Efeective viscosity of magnetc fluid suspensions [J]. Soviet Phys JETP, 1972, 34(6): 1291-1294

[4] OdenbachS.Magnetoviscous effects in ferrofluids [M], 2002, Tokyo, Springer

[5] GotohK., YamadaM.. Thermal convective in a horizontal layer of magnetic fluids [J]. J Phys Soc, 1982, 51(9): 3042-3048

[6] BERKOVSKY B M, MEDVEDEV V F, KRAKOV M S. Magnetic fluid, engineering application [M]. Oxford Univ Press, 1993.

[7] BakerG. R., CaflischR., SiegelM.. Singularity formation during Rayleigh-Taylor instability [J]. J Fluid Mech, 1993, 252: 477-501

[8] MOORE D W. Avortex method applied to interfacial waves [M]. HORNUNG H G, MULLER E A, eds. Vortex Motion, Vieweg & Sons, 1982.

[9] XieX.-ming.. Existence and uniqueness of analytic solution for Rayleigh-Taylor problem [J]. J Differential Equations, 2007, 237: 116-132

[10] LiXu.Nonliear evolution of R-T instability and SPH simulation [D], 2002, Beijing, China Academy of Engineering Physics

[11] QIN Cheng-sen, WANG Pei. Comparison beween Rayleigh-Taylor instability of compressible and incompressible fluid [J]. CNIC-01769, IAPCM-0042. (in Chinese)

[12] AnuchinaN. N., VolkovV. I., GordeychukV. A., Es’kovN. S., IlyutinaO. S., KozyrevO. M.. Numerical simulations of Tayleigh-Taylor anf Richtmyer-Meshkov instability using MAH-3 code [J]. Journal of Computational and Applied Mathematics, 2004, 168: 11-20

[13] ShiJ., ZhangY.-t., ChuC.-wang.. Resolution of high order WENO schemes for complicated flow structures [J]. Journal of Computational Physics, 2003, 186: 690-696

[14] Anjiali DeviS. P., HemamaliniP. T.. Nonlinear Rayleigh-Taylor instability of two superposed magnetic fluids under parallel rotation and a normal magnetic field [J]. Journal of Magnetism and Magnetic Materials, 2007, 314: 135-139

[15] KorovinV. M., KubasovA. A.. Tangential magnetic field induced structure in a thin layer of viscous magnetic fluid when developing Rayleigh-Taylor instability [J]. Journal of Magnetism and Magnetic Materials, 1999, 202: 547-553

[16] MohamedA. A., el ShehaweyE. F., El-dibY. O.. Electroviscoelastic Rayleigh-Taylor instability of Maxwell fluids (I): Effect of a constant tangential electric field [J]. J Phys A, 1994, 27: 3937

[17] SHU Chi-wang. Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws [R]. ICASE Report No.97-65.