arxiv: v1 [physics.plasm-ph] 2 Nov 2018

Transcription

1 Coherent nonlinear oscillations in magnetohydrodynamic plasma at Alfvén resonance Rupak Mukherjee, Rajaraman Ganesh, and Abhijit Sen Institute for Plasma Research, HBNI, Bhat, Gandhinagar , India arxiv: v1 [physics.plasm-ph] 2 Nov 218 Single fluid magnetohydrodynamic (MHD) equations have been studied through direct numerical simulation using pseudo-spectral method in two as well as three spatial dimensions at Alfvén resonance and exchange of energy between kinetic and magnetic variables have been observed. For a wide range of initial flow speeds or Alfven Mach number it is shown that the coherent nonlinear oscillation persist. As the Alfven Mach number is increased further, a tendency to mean field dynamo is seen. The dependency of the frequency of energy exchange decreases as the Alfvén Mach number increased leading to dynamo like phenomena and increases as the initial wave number excited in the fluid flow in increased. A finite mode representation of the single fluid MHD equations has been carried out in two spatial dimensions and similar behavior is found to hold. I. INTRODUCTION Understanding kinetic energy cascade for different initial spectrum has been a longstanding challenge in the field of fluid dynamics. Of late several direct numerical simulation (DNS) techniques have helped to analyse the kinetic energy spectra with the time evolution of any fluid governed by Navier-Stokes equation. The fluid equations represented as an infinite dimensional system is hard to solve analytically [1]. The complexity increases manifold when the fluid is charged and hence becomes susceptible to magnetic field. The self generated magnetic field dynamically affects the fluid flows through the Lorentz force term in the momentum equation. Thus the energy accumulated in any mode may be expected to find new route via nonlinear cascade to other modes through transformation of kinetic to magnetic energy and vice versa. In the presence of weak resistivity, the MHD model is known to predict irreversible conversion of magnetic energy into fluid kinetic energy (i.e. reconnection) as well as conversion of kinetic energy into mean large scale magnetic field (i.e. dynamo). Therefore it is interesting to ask oneself, for a given fluid type and magnetic field strength, are there fluid flow profiles which do neither - that is, neither does the flow generate mean magnetic field nor the magnetic field energy is converted to flow energy; instead there are nearly reversible coherent nonlinear oscillations. For a wide range of initial flow speeds or Alfven Mach number it is shown that the coherent nonlinear oscillation persist. As the Alfven Mach number is increased further, a tendency to mean field dynamo is seen. In this paper, we report numerical observation of interplay of kinetic and magnetic energy in two as well as three dimensional single fluid magnetohydrodynamic simulation for different set of initial flow profiles both chaotic and non-chaotic. A finite mode representation of rupak@ipr.res.in; rupakmukherjee1@gmail.com ganesh@ipr.res.in abhijit@ipr.res.in the field variables resulting in a set of nonlinear coupled ordinary differential equations has been presented in support of the pseudo-spectral DNS results. Our numerical as well as analytical observations indicate that nonlinear reciprocity of kinetic and magnetic energy can take place within the frame-work of single fluid magnetohydrodynamic representation of plasma and strongly depending on the M A. We start with a set of well known velocity flows as initial conditions and numerically time evolve the fluid as well as the magnetic variables self-consistently. The magnetic field we choose is uniform in space to start with. We look into the regime where the Alfvén speed is equal to the sound speed of the system. This choice facilitates most efficient coupling between magnetic and kinetic variables through the Lorentz force term. We find a continuous nonlinear exchange of energy between fluid and magnetic variables. The key results that we obtain from our DNS study backed with the results from an analytical galerkin representation are: For decaying two dimensional Orszag-Tang (OT) flow and Cat s Eye flow at Alfven resonance (M A = 1), an initial uniform magnetic field profile leads to nonlinear coherent oscillation between kinetic and magnetic modes. For decaying three dimensional Roberts flow, Cat s Eye flow and Arnold-Beltrami-Childress (ABC) flow at Alfven resonance (M A = 1), an initial uniform magnetic field profile leads to nonlinear coherent oscillation between kinetic and magnetic modes. The parameter dependency (Alfven Mach number) of this oscillation for three dimensional ABC flow is studied indicating boundary between two nonlinear events - coherent oscillation and dynamo effect. For externally forced OT and ABC flows, the nearly ideal magnetohydrodynamic plasma acts as a forced-relaxed system. The paper is arranged as follows. Section II summarises the equations that are time evolved in the code

2 2 MHD3D described in section III. In section IV we provide all the parameters in which the code is run and section V provides all the direct numerical simulation (DNS) result we have obtained in two and three dimensions. Section VI describes an analytical two dimensional finite mode representation of the DNS results and section VII summarises all the results and analysis done in this paper. II. GOVERNING EQUATIONS The basic equations governing the dynamics of the magnetohydrodynamic fluid are as follows: ρ t + (ρ u) = (1) [ ) ] (ρ u) + t ρ u u + (P + B2 I B 2 B = µ 2 u + ρf (2) B t + ( u B B ) u = η 2 B (3) P = Cs 2 ρ (4) In this system of equations, ρ, u, P and B are the density, velocity, kinetic pressure and magnetic field of a fluid element respectively. µ and η denote the coefficient of kinematic viscosity and magnetic resistivity. We assume µ and η are constants over space and time. The symbol represents the dyadic between two vector quantities. The kinetic Reynolds number (Re) and magnetic Reynolds number (Rm) are defined as Re = UL µ and Rm = UL η where U is the maximum velocity of the fluid system to start with and L is the system length. We also define the sound speed of the fluid as C s = U M s, where, M s is the sonic Mach number of the fluid. The Alfven speed is determined from V A = U M A where M A is the Alfven Mach number of the plasma. The initial magnetic field present in the plasma is determined from the relation B = V A ρ, where, ρ is the initial density profile of the fluid. III. DETAILS OF MHD3D CODE In order to simulate the plasma dynamics governed through the above-mentioned MHD equations, a DNS code, MHD3D, has been developed in-house at Institute for Plasma Research. MHD3D is an OpenMP parallel three dimensional weakly compressible, viscous, resistive magnetohydrodynamic code using pseudo-spectral technique to simulate a general scenerio of three dimensional magnetohydrodynamic turbulence problem. The pesudo-spectral technique, one of the most accurate CFD techniques available till today, uses the FFTW libraries [2] which is one of the fastest Fourier Transform libraries developed recently. This technique is applied to calculate the spatial derivatives and to evaluate the non-linear terms involved in the basic underlying equations with a standard de-aliasing by 2/3 truncation rule. The time derivative is solved using multiple time solvers viz. Adams-Bashforth, Runge-Kutta 4 and Predictor-Correcter algorithms and all the solvers have been checked to provide identical results. For two-dimensional hydro-dynamic case, we reproduced the results obtained earlier by Drazin [3], Ray [4] and Keppens et al [5] for a Kelvin-Helmholtz like scenario in both incompressible and compressible limits. The compressible magnetodydrodynamic solver in two dimensions is benchmarked with Keppens et al [5] and Orszag et al [6]. The three dimensional weakly compressible neutral fluid solver is benchmarked with the result of Samtaney et al[7] for the root mean square of velocity divergence and the skewness of velocity field for a decaying turbulence case. The three dimensional kinetic dynamo effect is matched with Galloway et al [8] in the case of ABC flow for a grid resolution of More benchmarking details and specifications about the code are given in previous literature [9, 1]. IV. PARAMETER DETAILS We choose L = 2π periodic box as our simulation domain and fix ρ = 1 and keep these parameters identical throughout our simulation. The initial magnitude of density (ρ ) is known to affect the dynamics and growth rate of an instability in a compressible neutral fluid [11, 12]. However, in this report we keep the initial density fixed (ρ = 1) for all the runs unless stated otherwise. For two dimensional runs we choose a grid resolution of N x = N y = 128 and a time stepping interval δt = 1 5. The kinetic and magnetic Reynolds numbers are chosen as Re = Rm = 1 4. We run our code for the incompressible limit (M s =.1). For three dimensional runs we choose a grid resolution of N x = N y = Nz = 64 and a time stepping interval δt = 1 4. The kinetic and magnetic Reynolds numbers are chosen as Re = Rm = 45. For all cases we run our code for the incompressible limit (M s =.1). For both two and three dimensional runs we check our code with higher grid resolution (N) and smaller time stepping (δt) and found no significant variation in the results of our test runs. The velocity profiles we choose are equilibrium stable solutions of Navier-Stokes equation and hence if there were no Lorentz force term present in the momentum equation, the flows would not have been disturbed even after a long time evolution. But back coupling of the magnetic variables through Lorentz force term affects the flow profile thereby giving rise to completely new dynamics in the system.

3 FIG. 1. evolution of shifted kinetic [u 2 (t) u 2 ()], magnetic [B 2 (t) B 2 ()] and total [u 2 (t) u 2 () + B 2 (t) B 2 ()] energy per grid for two dimensional Orszag-Tang flow with N = 128 2, U = 1, k = 1, M s =.1, M A = 1 and Re = Rm = 1 4 in the absence of external forcing. The energy keeps on switching between kinetic and magnetic modes keeping the total energy conserved. It is found that the time period of oscillation is T = We choose M A = 1 i.e. the mean velocity (U ) is equal to the Alfven speed of the system (V A ) (i.e. Alfven resonance) for both the two as well as the three dimensional runs. For two dimensional flows the initial magnetic field profile is chosen as B x = B y = B and that for three dimensional flows as B x = B y = B z = B. V. SIMULATION RESULTS A. Results for Decaying 2D Orszag-Tang Flow For two dimensional Orszag-Tang flow the velocity profile is chosen as u x = A sin(k y) u y = +A sin(k x) (5) with A = 1 and k = 1. The interchange of energy between kinetic and magnetic variables is plotted in Fig. 1. The x-axis represents the time evolution and the y-axis represents the shifted kinetic and magnetic energy where the shift indicates that the initial values of the corresponding variables are subtracted from the time evolution data. A long time evolution of the system with several cycles of oscillations with N = 256 2, δt = 1 6, U = 1, k = 1, M s =.1, M A = 1 and Re = Rm = 1 5 found to show numerical convergence. From figure in it evident that the time period of oscillation is T = FIG. 2. evolution of shifted kinetic, magnetic and total energy per grid for two dimensional Cat s Eye flow with N = 128 2, U = 1, A =.5, k = 1, M s =.1, M A = 1 and Re = Rm = 1 4 in the absence of external forcing. The energy keeps on switching between kinetic and magnetic modes keeping the total energy conserved. The decay of total energy is solely due to viscous and resistive effects. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. The time period of oscillation is T = B. Results for Decaying 2D Cat s Eye Flow For two dimensional Cat s Eye flow the velocity profile is chosen as u x = + sin(k x) cos(k y) A cos(k x) sin(k y) u y = cos(k x) sin(k y) + A sin(k x) cos(k y) (6) with A =.5, k = 1. The interchange of energy between kinetic and magnetic variables is plotted in Fig. 2. It is found that the time period of oscillation is T = C. Results for Forced 2D Orszag-Tang Flow [ ] A sin(kf y) f = α +A sin(k f x) (7) We choose, α =.1, A = 1 and k f = 1. The kinetic and magnetic energy continue to oscillate but the magnitudes of the peaks varies with time. The total (kinetic + magnetic) energy for the decaying cases were constant and the tiny decay was solely due to viscous and resistive effects. But for the forced case though the forcing acts for all the time over the whole system, the total energy oscillates with half the frequency of oscillation of the kinetic and magnetic energy [Fig. 3].

4 FIG. 3. evolution of shifted kinetic, magnetic and total energy per grid for two dimensional Orszag-Tang flow with N = 128 2, U = 1, k = 1, M s =.1, M A = 1 and Re = Rm = 1 4 in the presence of external forcing. The energy keeps on switching between kinetic and magnetic modes. The total energy also oscillates around the same frequency and the whole system acts as a forced-relaxed system. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables FIG. 4. evolution of shifted kinetic, magnetic and total energy per grid for three dimensional Roberts flow with N = 64 3, U =.1, A = B = C = 1, k = 1, M s =.1, M A = 1 and Re = Rm = 45 in the absence of external forcing. The energy keeps on switching between kinetic and magnetic modes keeping the total energy conserved. The decay of total energy is solely due to viscous and resistive effects. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. The time period of oscillation is T = 3.68 D. Results for Decaying 3D Roberts Flow For three dimensional Roberts flow the velocity profile is chosen as u x = U [A sin(k z)] u y = U [B sin(k x)] u z = U [C sin(k y)] (8) with A = B = C = 1, k = 1. The interchange of energy between kinetic and magnetic variables is plotted in Fig. 4. It is found that the time period of oscillation is T = F. Results for Decaying 3D ABC Flow For three dimensional Arnold-Beltrami-Childress flow we choose the velocity profile as u x = U [A sin(k z) + C cos(k y)] u y = U [B sin(k x) + A cos(k z)] u z = U [C sin(k y) + B cos(k x)] (1) with A = B = C = 1. The interchange of energy between kinetic and magnetic variables is plotted in Fig. 6. It is found that the time period of oscillation is T = E. Results for Decaying 3D Cat s Eye Flow For three dimensional 3D Cat s Eye flow the velocity profile is chosen as u x = U [B sin(k y)] u y = U [A sin(k x)] u z = U [A cos(k x) B cos(k y)] 3 5 (9) with A =, B = 2A. The interchange of energy between kinetic and magnetic variables is plotted in Fig. 5. It is found that the time period of oscillation is T = G. Analysis of Decaying 3D ABC Flow 1. Effect of Initial Wave-number We choose k = 1, 2, 4, 8, 16 keeping M A = 1. The frequency of oscillation increases with the increase of k value. However the viscous (µ 2 u) as well as the resistive terms (η 2 B) in Equation (12 and 3) affect the dynamics heavily when a higher wavenumber is excited. It is because the second derivative in the viscous and the resistive terms becomes more prominant when the velocity and magnetic field variables contain a higher magnitude of k [Table I, Fig. 9].

5 FIG. 5. evolution of shifted kinetic, magnetic and total energy per grid for three dimensional Cat s Eye flow with N = 64 3, U =.1, A =, B = 2A, k = 1, Ms =.1, MA = and Re = Rm = 45 in the absence of external forcing. The energy keeps on switching between kinetic and magnetic modes keeping the total energy conserved. The decay of total energy is solely due to viscous and resistive effects. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. The time period of oscillation is T = 3.51 k Period of Oscillation (T) TABLE I. Frequency of oscillation (ω) of kinetic and magnetic energy with mode number of excitation (k ) 2. Effect of Alfven Speed We study the effect of Alfven speed (V A ) on the energy exchange process by changing the Alfven Mach number M A for k f = 1. Fig.(1) represents the oscillation of kinetic energy for M A =.1,.2,.3,.4,.5 and Fig.(11) represents the oscillation of kinetic energy for M A =.5, 1., 1.5. From both the figures it is evident that the frequency of oscillation decreases linearly with the increase of M A. However if we further increase M A to 1, 1 and 1 the oscillation vanishes and saturation of magnetic energy is observed after a linear growth of magnetic energy at the cost of kinetic energy indicating a dynamo effect taking place [Fig.13]. We primarily believe these oscillations are due to excitation of Alfven waves. We know the time period oscillation of Alfven wave (T Alfven ) varies linearly with B which is directly proportional to V A = U M A. Thus if, M A is increased, T Alfven decreases linearly Fig.(1, 11). We analyse the frequency of oscillation for a series of FIG. 6. evolution of shifted kinetic, magnetic and total energy per grid for three dimensional Arnold-Beltrami- Childress flow with N = 64 3, U =.1, A = B = C = 1, k = 1, M s =.1, M A = 1 and Re = Rm = 45 in the absence of external forcing. The energy keeps on switching between kinetic and magnetic modes keeping the total energy conserved. The decay of total energy is solely due to viscous and resistive effects. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. The time period of oscillation is T = M A values and find good agreement with this linearity prediction [Table II, Fig. 12]. M A Period of Oscillation (T) TABLE II. period of oscillation (T ) of kinetic and magnetic energy with M A. H. Results for Forced ABC Flow Now we turn on the forcing on the velocity field. The forcing term we apply for all the time in the system is as follows: f = α A sin(k f z) + C cos(k f y) B sin(k f x) + A cos(k f z) (11) C sin(k f y) + B cos(k f x) We choose, α =.1, A = B = C = 1 and k f = 1. Similar to the externally forced two dimensional system [Fig.(3)], the total energy oscillates with time Fig.(14). Initially [Fig.(15)] the frequency of oscillation

6 k = 1; Kinetic k = 1; Magnetic k = 1; Total k = 2; Kinetic k = 2; Magnetic k = 2; Total k = 4; Kinetic k = 4; Magnetic k = 4; Total k = 8; Kinetic k = 8; Magnetic k = 8; Total k = 16; Kinetic k = 16; Magnetic k = 16; Total FIG. 7. evolution of shifted kinetic (solid), magnetic (dashed) and total (dotted) energy per grid for three dimensional Arnold-Beltrami-Childress flow with N = 64 3, U =.1, A = B = C = 1, k = 1 (red), 2 (blue), 4 (magenta), M s =.1, M A = 1 and Re = Rm = 45 in the absence of external forcing. Drag forces are found to affect the kinetic modes more than the magnetic modes as expected and explained in the text. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. of kinetic and magnetic energy were same though every second peak of magnetic energy had much lesser amplitude. The frequency of total energy is half of the frequency of kinetic energy. In intermediate time of evolution [Fig.(16)], the every second weak peak of magnetic energy smoothes out and in the late time evolution [Fig.(17)] we observe that the frequency of oscillation of magnetic energy becomes half of that of kinetic energy thereby becoming similar to that of total energy. We also plot the Lissajous curve of Kinetic and magnetic energy for the forced system in Fig. (18) FIG. 8. evolution of shifted kinetic (solid), magnetic (dashed) and total (dotted) energy per grid for three dimensional Arnold-Beltrami-Childress flow with N = 64 3, U =.1, A = B = C = 1, k = 8 (red), 16 (blue), M s =.1, M A = 1 and Re = Rm = 45 in the absence of external forcing. Drag forces are found to affect the kinetic modes more than the magnetic modes as expected and explained in the text. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. Period of Oscillation(T) k FIG. 9. Variation of time period of oscillation (T ) of kinetic and magnetic energy with k. [see Table I] VI. GALERKIN REPRESENTATION IN TWO DIMENSIONS In order to gain further understanding of the process, we resort to the finite mode representation for the field variables. The analytical theory is constructed for two dimensional decaying cases and is believed to hold for three dimensional decaying phenomena also. In two dimensions, we construct stream function (ψ) defined as u = ẑ ψ and magnetic vector potential ( A), defined as, B = A. In two dimensions the normalised incom- pressible MHD equations turns out as t 2 ψ + (ẑ ψ ) 2 ψ = (( B) B) (12) A t + (ẑ ψ )A = (13) We represent the finite mode expansion as, ψ(x, y) = ψ sin kx + e iky (ψ 1 + ψ 3 cos kx) + e iky (ψ 1 + ψ 3 cos kx) (14) A(x, y) = A sin kx + e iky (A 1 + A 3 cos kx) + e iky (A 1 + A 3 cos kx) (15)

7 7 Shifted M A =.1 M A =.2 M A =.3 M A =.4 M A =.5 Period of Oscillation (T) FIG. 1. evolution of shifted kinetic energy per grid for three dimensional Arnold-Beltrami-Childress flow with N = 64 3, U =.1, A = B = C = 1, k = 1, M s =.1, M A =.1 (red),.2 (greed),.3 (blue),.4 (magenta),.5 (black) and Re = Rm = 45 in the absence of external forcing. The frequency of oscillation decreases with the increase of M A. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. Shifted M A =.5 M A = 1. M A = M A FIG. 12. Variation of time period of oscillation (T ) of kinetic and magnetic energy with M A. [see Table II] Energy e-5 MA = 1; Kinetic MA = 1; Magnetic 1e-6 MA = 1; Kinetic MA = 1; Magnetic 1e-7 MA = 1; Kinetic MA = 1; Magnetic 1e FIG. 13. Saturation of magnetic energy with higher values of M A indicating dynamo effect FIG. 11. evolution of shifted kinetic energy per grid for three dimensional Arnold-Beltrami-Childress flow with N = 64 3, U =.1, A = B = C = 1, k = 1, M s =.1, M A =.5 (red), 1 (blue), 1.5 (magenta) and Re = Rm = 45 in the absence of external forcing. The frequency of oscillation decreases with the increase of M A. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables. where k is the wave vector corresponding to the largest scale length of the system. Substituting Equation 14 and 15 into Equation 12 and 13 we obtain the time evolution of the galerkin truncated modes. The evolution of various modes of stream function and magnetic vector potential can thus be written as, dψ = ik 2 (ψ dt 3ψ 1 ψ 3 ψ1) (16) dψ 1 = i dt 2 k2 (ψ ψ 3 + A A 3 ) (17) dψ 3 = dt (18) da = ik 2 (A 3 ψ1 A dt 3ψ 1 + A 1 ψ3 A 1ψ 3 ) (19) da 1 = i dt 2 k2 (A ψ 3 A 3 ψ ) (2) da 3 = ik 2 (A ψ 1 A 1 ψ ) dt (21) We start with a case similar to that of our initial condition used for the DNS code (ψ = 1, ψ 1 = = ψ 3, A = A 1 = A 3 = 1 and k = 1) and time evolve the above set of ordinary differential equations using RK4 method. In Fig. 19 we show the time evolution of real part of ψ 1

8 FIG. 16. Intermediate time evolution obtained from Fig FIG. 14. evolution of shifted kinetic, magnetic and total energy per grid for three dimensional Arnold-BeltramiChildress flow with N = 643, U =.1, A = B = C = 1, k = 1, Ms =.1, MA = 1 and Re = Rm = 45 in the presence of external forcing. The energy keeps on switching between kinetic and magnetic modes. The total energy also oscillates around the same frequency and the whole system acts as a forced-relaxed system. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables FIG. 17. Late time evolution obtained from Fig FIG. 15. Initial time evolution obtained from Fig and A1. The variables (ψ1 and A1 ) are found to keep exchanging their values as the time evolves. We also find the increment of frequency with the increase of k in two dimensional system similar to that in three dimensional system [Fig 2,21. We further explore the parameter space by choosing lower values of A indicating the higher MA cases where the dynamo effect happened in the 3D DNS runs. No saturation of magnetic energy is observed in the low degree of freedom model, indicating that the dynamo effect requires more number of modes to occur which is beyond the scope of the present analytical model FIG. 18. evolution of kinetic and magnetic energy per grid for three dimensional Arnold-Beltrami-Childress flow with N = 643, U =.1, A = B = C = 1, k = 1, Ms =.1, MA =.1 and Re = Rm = 45 in the presence of external forcing. The shift denotes the initial values of the kinetic and magnetic energies (at time t = ) are subtracted from the time evolution data of the corresponding variables.

9 9 1 ψ 1 (t) - ψ 1 () A 1 (t) - A 1 () 1 ψ 1 (t) - ψ 1 () A 1 (t) - A 1 () [ψ 1 (t) - ψ 1 ()]; [A 1 (t) - A 1 ()] [ψ 1 (t) - ψ 1 ()]; [A 1 (t) - A 1 ()] (t) (t) FIG. 19. Galerkin truncated two dimensional MHD equations representing similar nonlinear coherent oscillations observed in DNS results with k = 1. FIG. 21. Galerkin truncated two dimensional MHD equations representing similar nonlinear coherent oscillations observed in DNS results with k = 4. VII. SUMMARY AND FUTURE WORK [ψ 1 (t) - ψ 1 ()]; [A 1 (t) - A 1 ()] ψ 1 (t) - ψ 1 () A 1 (t) - A 1 () In this paper we show the single fluid magnetohydrodynamic equations supports energy exchange between kinetic and magnetic variables over a long enough period in both two and three dimensions. Thus it is shown that the energy exchange phenomena takes place for both non-chaotic as well as chaotic fields. The time period of oscillation of the energy exchange is found to vary linearly with Alfven Mach number (M A ). Also frequency is seen to increase with the initial wave-number (k) excited within the flow. In presence of an external forcing identical to the initial velocity field, the system is found to act as a forced-relaxed system showing the periodic variation of the total energy of the system as well as exchange of energy between kinetic and magnetic modes. The three dimensional calculations are performed over 64 3 grid resolution. A better spatial resolution is expected to be capable of simulating higher fluid and magnetic Reynold s number runs. Thus the periodicity of energy exchange in three dimensional flows could be observed for larger number of cycles for high resolution runs (t) VIII. ACKNOWLEDGEMENT FIG. 2. Galerkin truncated two dimensional MHD equations representing similar nonlinear coherent oscillations observed in DNS results with k = 2. The development as well as benchmarking of MHD2D and MHD3D have been done at Udbhav and Uday clusters at IPR. [1] P. Morrison, Encyclopedia of mathematical physics pp (26). [2] M. Frigo and S. G. Johnson, Proceedings of the IEEE 93, 216 (25), special issue on Program Generation, Optimization, and Platform Adaptation. [3] P. Drazin, Journal of Fluid Mechanics 1, 571 (1961). [4] T. P. Ray, Monthly Notices of the Royal Astronomical Society 198, 617 (1982).

10 1 [5] R. Keppens, G. Tóth, R. Westermann, and J. Goedbloed, Journal of Plasma Physics 61, 1 (1999). [6] S. A. Orszag and C.-M. Tang, Journal of Fluid Mechanics 9, 129 (1979). [7] R. Samtaney, D. I. Pullin, and B. Kosović, Physics of Fluids 13, 1415 (21). [8] D. Galloway and U. Frisch, Geophysical & Astrophysical Fluid Dynamics 36, 53 (1986). [9] R. Mukherjee, A. Gupta, and R. Ganesh, arxiv preprint arxiv: (218). [1] R. Mukherjee, R. Ganesh, V. Saini, U. Maurya, N. Vydyanathan, and B. Sharma, arxiv preprint arxiv: (218). [11] B. Bayly, C. Levermore, and T. Passot, Physics of Fluids A: Fluid Dynamics 4, 945 (1992). [12] D. Terakado and Y. Hattori, Physics of Fluids 26, 8515 (214).