Next Article in Journal / Special Issue
Isolation Properties of Low-Profile Magnetorheological Fluid Mounts
Previous Article in Journal
Euler’s Equation of Continuity: Additional Terms of High Order of Smallness—An Overview
Previous Article in Special Issue
Oscillating Magnetohydrodynamic Stokes Flow between Porous Plates with Spatiotemporally Periodic Reabsorption
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Frequency Power Spectra of Global Quantities in Unsteady Magnetoconvection

Department of Physics, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
*
Author to whom correspondence should be addressed.
Fluids 2021, 6(4), 163; https://doi.org/10.3390/fluids6040163
Submission received: 5 March 2021 / Revised: 5 April 2021 / Accepted: 8 April 2021 / Published: 19 April 2021
(This article belongs to the Special Issue Fluids in Magnetic/Electric Fields)

Abstract

:
We present the results of direct numerical simulations of power spectral densities for kinetic energy, convective entropy, and heat flux for unsteady Rayleigh–Bénard magnetoconvection in the frequency space. For larger values of frequency, the power spectral densities for all the global quantities vary with frequency (f) as f 2 . The scaling exponent is independent of Rayleigh number, Chandrasekhar’s number, and thermal Prandtl number.

1. Introduction

The temporal fluctuations [1] of spatially averaged (or, global) quantities are of interest in several fields of research, including turbulent flows [2,3,4,5,6,7], nanofluids [8], biological fluids [9,10], geophysics [11,12], and phase transitions [13,14]. The probability density function (PDF) of the temporal fluctuations of thermal flux in turbulent Rayleigh–Bénard convection (RBC) was found to have normal distribution with slight asymmetries near the tails. Auma i ^ tre and Fauve [5] also showed that the statistical properties of locally measured and globally averaged temperature field are quite different. They found that the power spectra of Nusselt number ( Nu ), which is a measure of heat flux across the fluid layer in RBC, varies with frequency (f) as f 2 at higher frequencies. The scaling exponent for the temperature field that is measured near the lower and upper boundaries is different from the one that is measured in the central part of the experimental cell. The direct numerical simulations (DNS) of the Nusselt number ( Nu ) also showed the similar behaviour in the presence of Lorentz force [15] in water-based nanofluids ( Pr = 4.0 ). The power spectral density (PSD) of the thermal flux in the frequency (f) space [2,5,15,16] was found to vary as f 2 . The frequency spectrum for the temperature field measured near the horizontal plates in Rayleigh–Bénard magnetoconvection [17] shows the exponent to vary between 6 and 4 . The frequency spectra of kinetic energy and convective entropy for the problem of magnetoconvection are rarely available, either experimentally or numerically. The frequency spectrum for the Nusselt number [15] was recently computed only for one value of the thermal Prandtl number.
In this work, we present the results that were obtained by DNS for the frequency spectra of three global quantities: spatially averaged kinetic energy per unit mass (E), convective entropy per unit mass ( E Θ ), and Nusselt number ( Nu ) in unsteady Rayleigh–Bénard magnetoconvection (RBM) [18,19,20] for several values of Rayleigh number ( Ra ), Prandtl number ( Pr ), and Chandrasekhar number ( Q ). The objective is to study the statistical properties of the fluctuating global quantities in Rayleigh-Bénard magnetoconvection. The kinetic energy per unit mass as well as the convective entropy per unit mass are found to vary with frequency as f 2 at relatively higher frequencies. This behaviour does not depend on the Rayleigh number ( Ra ), Prandtl number ( Pr ), and Chandrasekhar’s number ( Q ).

2. Governing Equations

The physical system consists of a thin layer of a Boussinesq fluid (e.g., liquid metals, melts of some alloys (i.e., N a N O 3 melt), nanofluids, etc.) of density ρ 0 and electrical conductivity σ confined between two horizontal plates, which are made of electrically non-conducting, but thermally conducting materials. The lower plate is heated uniformly and the upper plate is cooled uniformly, so that an adverse temperature gradient β is maintained across the fluid layer. A uniform magnetic field B 0 is applied in the vertical upward direction, which is also considered the positive direction of the z-axis. The x and y axes are in the horizontal plane with origin of the coordinate system taken on the lower plate. In the basic state, the fluid conducts heat without any motion. The stratification of the steady temperature field [ T s ( z ) ] and the fluid density [ ρ s ( z ) ], in the conduction state [18], are given as:
T s ( z ) = T b β z ,
ρ s ( z ) = ρ 0 1 + α T b T s ( z ) ,
where T b and ρ 0 are the temperature and density of the fluid at the lower plate, respectively. The fluid pressure [ P s ( z ) ], in conductive state, is:
P s ( z ) = C ρ 0 g z + 1 2 α β z 2 + B 0 2 2 μ 0 ,
where g is the acceleration due to gravity and μ 0 is the permeability of free space. The fluid pressure in the conductive state consists of hydrostatic, thermal, and magnetic pressures. The constant of integration (C) may be determined if the value of pressure at the upper plate [ P s ( z = d ) ] is known. If we take P s ( z = d ) = P 0 , where P 0 is a constant (e.g., air pressure at the upper plate), the constant C turns out to be
C = P 0 + B 0 2 2 μ 0 + ρ 0 g d + 1 2 α β d 2 .
The fluid pressure, in the basic conductive state, then reads as:
P s ( z ) = P 0 + ρ 0 g d z + 1 2 α β d 2 z 2 .
As soon as the temperature gradient across the fluid layer is raised above a critical value β c for fixed values of all fluid parameters (kinematic viscosity ν , thermal diffusivity κ , and thermal expansion coefficient α ) and the externally imposed magnetic field ( B 0 ), the convection sets in. All the fields are perturbed due to magnetoconvection and they may be expressed as:
ρ s ( z ) ρ ˜ ( x , y , z , t ) = ρ s ( z ) + δ ρ ( x , y , z , t ) ,
T s ( z ) T ( x , y , z , t ) = T s ( z ) + θ ( x , y , z , t ) ,
P s ( z ) P ( x , y , z , t ) = P s ( z ) + p ( x , y , z , t ) ,
B 0 B ( x , y , z , t ) = B 0 + b ( x , y , z , t ) ,
where v ( x , y , z , t ) , p ( x , y , x , t ) , θ ( x , y , z , t ) and b ( x , y , z , t ) are the fluid velocity, perturbation in the fluid pressure, the convective temperature and the induced magnetic field, respectively, due to magnetoconvection. The perturbative fields are made dimensionless by measuring all lengths in units of the clearance d between two horizontal plates, which is also the thickness of the fluid layer. Time is measured in units of the free fall time τ f = 1 / α β g . The convective temperature field θ and induced magnetic field b are made dimensionless by β d and B 0 Pm , respectively. The magnetoconvective dynamics is then described by the following dimensionless equations:
D t v = p + Pr Ra 2 v + Q Pr Ra z b + Pm ( b · ) b + θ e 3 ,
Pm D t b = Pm ( b · ) v + Pr Ra 2 b + z v ,
D t θ = 1 Ra Pr 2 θ + v 3 ,
· v = · b = 0 ,
where D t t + ( v · ) is the material derivative. The unit vector e 3 is directed vertically upward. Because the magnetic Prandtl number is very small ( Pm 10 5 ) for all terrestrial fluids, we set Pm equal to zero in the above set of hydrodynamic equations. The Navier–Stokes equation Equation (10) then takes the form, as given below.
D t v = p + Pr Ra 2 v + Q Pr Ra z b + θ e 3 .
The induced magnetic field ( b ) is slaved to the velocity field ( v ), as Equation (11) is simplified to
2 b = Ra Pr z v .
The fluid flow due to magnetoconvection, in the limit of Pm 0 , is described by the set of Equations (12)–(15). We consider the idealized boundary (stress-free) conditions for the velocity field on the horizontal boundaries. Relevant boundary conditions [18,21] at the horizontal plates, which are located at z = 0 and z = 1 , are:
v 1 z = v 2 z = v 3 = b 1 = b 2 = b 3 z = θ = 0 .
All of the fields are considered periodic in the horizontal plane. The fluid dynamics, as Pm 0 , is controlled by three dimensionless parameters: (1) Rayleigh number ( Ra = α β g d 4 ν κ ), (2) Prandtl number ( Pr = ν κ ), and (3) Chandrasekhar’s number ( Q = σ B 0 2 d 2 ρ 0 ν ). The critical values of Rayleigh number [ Ra c ( Q ) ] and critical wave number [ k c ( Q ) ] are [18]:
Ra c ( Q ) = π 2 + k c 2 k c 2 ( π 2 + k c 2 ) 2 + π 2 Q ,
k c ( Q ) = π a + + a 1 2 , w h e r e
a ± ( Q ) = 1 4 1 2 + Q π 2 ± 1 2 + Q π 2 2 1 4 1 2 1 3 .
The kinetic energy (E) and convective entropy ( E Θ ) per unit mass are defined as:
E ( t ) = 1 V 0 2 π k c d x 0 2 π k c d y 0 1 1 2 v 2 d z ,
E Θ ( t ) = 1 V 0 2 π k c d x 0 2 π k c d y 0 1 1 2 θ 2 d z .
The Nusselt number ( Nu ), which is the ratio of total heat flux and the conductive heat flux across the fluid layer, is given as:
Nu ( t ) = 1 + Ra Pr V 0 2 π k c d x 0 2 π k c d y 0 1 v 3 θ d V .
The hydromagnetic system of Equations (12)–(15) presented here may also be used to investigate magnetoconvection in nanofluids with a low concentration of non-magnetic metallic nanoparticles [15] in water. A homogeneous suspension of nanoparticles in a viscous fluid works as a nanofluid. Because the properties of a nanofluid depend on those of the base fluid and the nanoparticles, the effective values of dimensionless parameters would depend on both. All of the fluid parameters are may be replaced by their effective values in the presence of nanoparticles. If ϕ is the volume fraction of the spherically shaped nanoparticles, then the effective form of the density and electrical conductivity of a nanofluid may be expressed as:
ρ n f = ( 1 ϕ ) ρ f + ϕ ρ p ,
σ n f = ( 1 ϕ ) σ f + ϕ σ p ,
where ρ f and σ f are the density and electrical conductivity of the base fluid, respectively. Here, ρ p is the density and σ p is the electrical conductivity of nanoparticles. The effective thermal conductivity ( K n f ) of a nanofluid [22] for smaller values of ϕ is expressed as:
K n f = K f ( K p + 2 K f ) 2 ϕ ( K f K p ) ( K p + 2 K f ) + ϕ ( K f K p ) ,
where K f and K p are the thermal conductivities of the base fluid and that of sphere shaped nanoparticles, respectively. Similarly, the product of density and specific heat capacity [ ( ρ c V ) n f ] may be expressed through the following relation [23]:
( ρ c V ) n f = ( 1 ϕ ) ( ρ c V ) f + ϕ ( ρ c V ) p .
The effective dynamic viscosity ( μ n f ) of a nano-fluid [24] may also be expressed as:
μ n f = μ f ( 1 ϕ ) 2.5 .
The kinematic viscosity ( ν n f ) and thermal diffusivity ( κ n f ) for nanofluids may then be estimated while using the formulae:
ν n f = μ n f ρ n f = μ f ( 1 ϕ ) 2.5 ( 1 ϕ ) ρ f + ϕ ρ p ,
κ n f = K n f ( ρ c V ) n f = K f ( K p + 2 K f ) 2 ϕ ( K f K p ) ( 1 ϕ ) ( ρ c V ) f + ϕ ( ρ c V ) p ( K p + 2 K f ) + ϕ ( K f K p ) .
The effective thermal Prandtl number ( Pr n f ) for nanofluids is then defined as Pr n f = ν n f κ n f . The definitions of Rayleigh number ( Ra n f ) and Chandrasekhar’s number ( Q n f ) are also modified for nanofluids, and they read as:
Ra n f = α β g d 4 ν n f κ n f , Q n f = σ n f B 0 2 d 2 ( ρ 0 ) n f ν n f .
Substituting the values of dimensionless parameters Pr , Ra and Q by Pr n f , Ra n f and Q n f , respectively, the set of Equations (12)–(15) describes magnetoconvection in nanofluids with low concentration of nanoparticles in a base fluid. The thermal Prandtl number ( Pr n f ) of water-based nanofluids may be varied between 4.0 and 6.5 by varying the volume fraction ( ϕ ) of the spherical copper nanoparticles [15] between 8 % to 0.2 % .

3. Direct Numerical Simulations

Direct numerical simulations of the magnetoconvective flows for different values of the dimensionless parameters are carried out using pseudo-spectral method. All of the perturbative fields are expanded such that they are consistent with the boundary conditions. Perturbations are expanded as:
Ψ ( x , y , z , t ) = l , m , n Ψ l m n ( t ) e i k c ( l x + m y ) cos ( n π z ) ,
Φ ( x , y , z , t ) = l , m , n Φ l m n ( t ) e i k c ( l x + m y ) sin ( n π z ) ,
where Ψ ( x , y , z , t ) = [ v 1 , v 2 , p ] and Φ ( x , y , z , t ) = [ v 3 , θ ] . The time dependent Fourier amplitudes of these fields are denoted by Ψ l m n ( t ) = [ U l m n , V l m n , P l m n ] and Φ l m n ( t ) = [ W l m n , Θ l m n ] , where l, m, and n are integers. The horizontal wave vector of all the perturbations is k H = ( l e 1 + m e 2 ) k c ( Q ) , where e 1 and e 2 are the unit vectors along positive directions of the x- and y-axes, respectively. All of the numerical simulations are carried out in a three dimensional simulation box of size L × L × 1 , where L = 2 π / k c ( Q ) . The value of k c is computed using the expression for the critical wave number k c ( Q ) [see Equations (18) and (19)]. The continuity equations decide the possible values of the integers l , m , n . They can take values that satisfy the following equation.
i l k c ( Q ) U l m n + i m k c ( Q ) V l m n + n π W l m n = 0 .
The computation of non-linear terms ( v · ) v and v · θ are computed while using Fast Fourier Transformation (FFT). It is done using the following steps:
(i)
Real space variables v ( x , y , z ) = ( v 1 , v 2 , v 3 ) and θ ( x , y , z ) are computed at a given time t using inverse FFT of v ( k ) and θ ( k ) , where k = l k c e 1 + m k c e 2 + n π e 3 .
(ii)
The multiplication of field variables v i ( x , y , z ) v j ( x , y , z ) and v i ( x , y , z ) θ ( x , y , z ) for ( i , j = 1 , 2 , 3 ) are done at each grid point of the simulation box.
(iii)
FFT[ v i ( x , y , z ) v j ( x , y , z ) ] and FFT[ v i ( x , y , z ) θ ( x , y , z ) ] are computed using the package FFTW.
(iv)
Subsequently, the terms i k j × F F T [ v i ( x , y , z ) v j ( x , y , z ) ] and i k j × F F T [ v i ( x , y , z ) θ ( x , y , z ) ] with j = 1 , 2 as well as i k j × F F T [ v i ( x , y , z ) v j ( x , y , z ) ] and i k j × F F T [ v i ( x , y , z ) θ ( x , y , z ) ] with j = 3 are computed.
The aliasing error is removed using 2 / 3 -rule [25]. The integration in time is performed using a standard fourth order Runge–Kutta (RK4) scheme. The time step of RK4 integration scheme was monitored, so that the Courant–Friedrichs–Lewy (CFL) condition was satisfied for all times. The time step of 10 3 (in dimensionless units) is used for integration. The data points for the temporal signal of all global quantities are recorded at all time steps. The grid size was chosen, such that the smallest dissipative (Kolmogorov) scale is resolved. The thermal dissipative scale was the smallest for our simulations. A resolutions of 128 × 128 × 128 was good enough for several of the simulations that are presented here. Each of the simulations was usually carried out for more than 600 dimensionless time units. For P r > 1 , we carried the simulations with higher resolutions of 256 × 256 × 256 . Das and Kumar [15] test the code for a magnetoconvection problem, which give more details about the simulations. Once a simulation for fixed values of dimensionless parameters Ra , Pr and Q is completed, the same procedure is repeated with a new set of Ra , Pr , and Q .

4. Results and Discussion

The simulations are done for several values of thermal Prandtl number ( 0.1 Pr 6.4 ). These values of Pr are relevant for Earth’s liquid outer core [26]. They are also relevant for the problem of crystal growth [27] and water-based nanofluids [8]. The Rayleigh number is varied in a range 7.0 × 10 4 Ra 3.04 × 10 6 , while the Chandrasekhar’s number is varied in the range 50 Q 10 3 . In numerical simulations, the value of Δ T = T b T u is always mapped to unity. The total temperature field [ T ˜ ( x , y , z , t ) ], in the dimensionless units, then reads as: T ˜ ( x , y , z , t ) = 1 z + θ ( x , y , z ) . Figure 1 shows the combined plot of the total temperature field [ T ˜ ( x , y , z , t ) ] and the velocity field [ v ( x , y , z , t ) ] simultaneously at two different instants for a water based nanofluids ( Pr n f = 6.4 ) for Ra n f = 5 × 10 5 and Q n f = 100 . The dimensions of the simulation box are L × L × 1 , where L = 2 π k c ( Q n f ) = 1.697 for Q n f = 100 . The colour bar that is given in each of the viewgraphs shows the temperature field at different points in the convecting fluid. Red, blue, and other colours stand for the hottest, coolest, and intermediate temperatures, respectively. The arrows show the flow directions at the outer surfaces of the simulation box. Thin thermal boundary layers near the lower and upper surfaces of the simulations box are clearly visible. There are generations of thermal plumes (red coloured patches) at the lower boundaries if the temperature gradient is large enough. They are of different sizes and they appear quite irregularly at different locations on the horizontal plate. This also contributes to the fluctuations of fields. There is no viscous boundary layer because of the use of stress-free boundary conditions on the velocity field. Two viewgraphs in Figure 1 confirm the unsteady magnetoconvection. The flow structure is changing with time drastically. Figure 2 shows the variation of flow structures for Pr n f = 6.4 and Ra n f = 5 × 10 5 and for two different values of Chandrasekhar number ( Q n f ). Figure 2a shows the total temperature field [ T ˜ ( x , y , z ) ] and the velocity field [ v ( x , y , z ) ] in the simulation at a randomly chosen time for Q n f = 100 . The flow structure that is shown in Figure 2b is for Q n f = 250 .
Figure 3 shows the temporal variations of three global quantities for Ra = 5.0 × 10 5 , Pr = 4.0 and for two different values of Q : (1) the kinetic energy per unit mass (E), (2) the convective entropy per unit mass ( E Θ ), and (3) the Nusselt number ( Nu ). All fields are recorded on all grid points and then the global quantities are computed by averaging over the three-dimensional simulation box. The first two sets of curves (from the top) show the variation of E with dimensionless time for two different values of Q . The red curve is for Q = 100 and the blue curve is for Q = 400 . The energy signal varies irregularly with time. They show appreciable fluctuations. The mean of kinetic energy decreases with increase in Q . The fluctuations of kinetic energy also decreases with an increase in Q . Curves in the third and fourth rows (from the top) show the temporal variation of E Θ , while curves in the fifth and sixth rows display the temporal signal of Nu . They also vary irregularly with time. The mean values of the convective entropy per unit mass as well as the Nusselt number decrease with an increase in Q . Their fluctuations also decrease, as Q is raised.
Figure 4 displays the power spectrum densities (PSD) for the global (spatially averaged) quantities in the frequency (f) space for several values of Ra , Pr and Q . The PSDs of the fluid speed [ E ( f ) = | v ( f ) | 2 ] are shown in Figure 4a. As the thermal energy is injected slowly in the fluid, the spectra show a small slope at low frequencies. In a small frequency window (approximately, 0.04 f 1.0 ), the slope of curves E ( f ) f on the log-log scale varies between 3.2 to 5.1 . The energy spectra [ E ( f ) ] have more noise in this frequency range. It may be due to irregularity in the size and frequency of thermal plumes from a particular position on the lower boundary. The spectra for nanofluids ( Pr n f = 4.0 and 6.4 ) show more noise. Smaller values of the thermal diffusivity of the fluid lead to an enhancement of the thermal noise, which makes the spectra nosier at lower frequencies. However, the E ( f ) is found to have insignificant noise for 1 < f < 200 . The frequency spectrum of energy per unit mass [ E ( f ) ] decays slowly in this regime. This may be due to the box averaging of the perturbative fields, and it is different than the behaviour near horizontal boundaries. The PSD [ E ( f ) ] of the kinetic energy scales with frequency (f) almost as f α with α 2 . The scaling behaviour is valid for more than two decades. The PSD shows a clear scaling behaviour for f > 1 . The scaling exponent is independent of Pr , Ra and Q in this frequency window. Table 1 gives the numerically computed values of the exponent ( α ) for different values of Ra , Pr , and Q . A similar scaling law [ E ( f ) f 2 ] was also observed in rotating Rayleigh–Bénard convection [16].
Figure 4b shows the PSDs of the convective entropy [ E Θ ( f ) = | θ ( f ) | 2 ] of the fluid in the frequency space for different values of Ra , Pr and Q . Its power spectra is also noisy in the dimensionless frequency range 0.04 < f < 1.0 . The slope on the log-log scale varies between 5.9 and 6.4 . However, for f > 1.0 , E Θ scales with frequency as f β with β 2 . Table 1 lists the numerically computed values of the exponent β . Interestingly, the power spectra of the temperature fluctuations are also found to vary as f 2 in the rotating RBC experiments [2,16]. The temperature of convecting fluid near the horizontal plates in RBM was measured in experiments on magnetoconvection [17]. The frequency spectra of the temperature field near the horizontal plates vary clearly as f 2 for Q = 0 . Even for very small values of Q ( Q 10 3 10 4 ), this scaling regime shrinks drastically. The data points for the temperature field were measured locally in this experiment. The frequency spectra of local quantities near the boundary are known to decay faster [5]. The convective entropy is a global quantity and it is computed by taking an average over the simulation box. That is why the power law for the temperature field recorded locally in the experiment [17] may reflect the behaviour in a thin boundary layer near the horizontal plates. In addition, we have used stress-free conditions on the velocity field on the horizontal plates. We only have thermal boundary layers. These may be reasons for disagreement in the power behaviour.
Figure 4c shows the power spectral densities of the thermal flux [Nusselt number, Nu ( f ) ] for several values of Ra , Pr , and Q . PSDs of the Nusselt number also show the scaling behaviour. The PSDs are noisy, as in the case of energy and entropy signals, for dimensionless frequencies 0.04 < f < 1.0 , if Pr > 1 . The scaling exponent varies between 4.5 to 6.4 in this frequency window. However, for higher values of dimensionless frequency ( 1 < f < 200 ), the spectra for thermal flux [ Nu ( f ) ] also shows very clear scaling: Nu ( f ) f γ , with γ 2 . Table 1 shows the values of the exponent γ that is computed in DNS. The measurements of the spectra of thermal flux in RBC also shows the similar scaling law [5]. The scaling behaviour for all global quantities for Rayleigh-Bénard magnetoconvection shows a similar power law. This is also observed in experiments in different systems [2,5] as well as in numerical simulations [15,16]. The power spectra of systems showing self-organised criticality [28] scale with frequency as f 1 , which is different than the scaling behaviour in the hydrodynamic systems discussed here.
The scaling law showing the variation of the power spectra as f 2 starts at a critical frequency ( f c ) for different values of Q . Figure 5 shows the variation of f c for E ( f ) , E Θ ( f ) , and Nu ( f ) with Q for two different sets of Ra and Pr . The critical frequency f c ( E ) decreases, as Q is increased (see Figure 5a). Figure 5b,c show the variations of f c ( E Θ ) and f c ( Nu ) , respectively, with Q . The values of critical frequencies are different for E ( f ) , E Θ ( f ) and Nu ( f ) . However they all decrease with an increase in Q for given values of Ra and Pr . A vertical magnetic field is known to delay the onset of fluid flow in RBC [18]. The external magnetic field suppresses the fluid motion, which leads to a decrease in f c ( Q ) , as Q is raised.

5. Conclusions

The results of DNS on Rayleigh–Bénard magnetoconvection show that frequency spectra of global quantities (the kinetic energy E ( f ) , convective entropy E Θ ( f ) , and the Nusselt number Nu ( f ) ) vary differently in different frequency windows. At low frequencies, the thermal energy enters slowly into the moving fluid. The frequency spectra, therefore have a small slope at low frequencies. In a narrow intermediate frequency range, the spectral densities of global quantities decrease sharply with frequency (f). They show power law behaviour. However, the scaling exponents are dissimilar for different global quantities. The energy spectrum [ E ( f ) ] varies with the frequency (f) as E ( f ) f ξ , with 3.2 ξ 5.1 . Scaling behaviour for the convective entropy and Nusselt number are E Θ ( f ) f η , with 5.9 η 6.4 and Nu ( f ) f ζ with 4.5 ζ 6.4 , respectively. For frequencies above a critical frequency f c , all global quantities show universal scaling behaviour. The critical values f c ( E ) , f c ( E Θ ) , and f c ( Nu ) are different for kinetic energy, convective entropy, and the Nusselt number, respectively. The spectra of all three global quantities scale with frequency as f 2 . The scaling exponent of 2 at higher frequencies are independent of Ra , Pr , and Q . The scaling behaviour of the global quantities at higher frequencies appears to be a universal feature in turbulent flows. The results that are presented here may also be relevant for some problems of geophysics, water-based nanofluids, crystal growth, and magnetohydrodynamics.

Author Contributions

Conceptualization: K.K.; Methodology: S.D. and K.K.; Simulations: S.D.; Analysis: K.K. and S.D.; Writing: K.K. and S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received partial funding from IIT Kharagpur, Kharagpur West Bengal, India.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available on request from SD.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Botet, R.; Ploszajczak, M. Universal Fluctuations (The Phenomenology of Hadronic Matter): World Scientific Lecture Notes in Physics; World Scientific: Singapore, 2002. [Google Scholar]
  2. Boubnov, B.M.; Golitsyn, G.S. Temperature and velocity field regimes of convective motions in a rotating plane fluid layer. J. Fluid Mech. 1990, 219, 215–239. [Google Scholar] [CrossRef]
  3. Fauve, S.; Laroche, C.; Castaing, B. Pressure fluctuations in swirling turbulent flows. J. Phys. II 1993, 3, 271–278. [Google Scholar] [CrossRef]
  4. Niemela, J.J.; Skrbek, L.; Sreenivasan, K.R.; Donnelly, R.J. Turbulent convection at very high Rayleigh numbers. Nature 2000, 404, 837–840. [Google Scholar] [CrossRef] [PubMed]
  5. Aumaître, S.; Fauve, S. Statistical properties of the fluctuations of the heat transfer in turbulent convection. Europhys. Lett. 2003, 62, 822–828. [Google Scholar] [CrossRef]
  6. Tanriverdi, V.; Tilgner, A. Global fluctuations in magnetohydrodynamic dynamos. New J. Phys. 2011, 13, 033019. [Google Scholar] [CrossRef]
  7. Gerolymos, G.; Vallet, I. Pressure, density, temperature and entropy fluctuations in compressible turbulent plane channel flow. J. Fluid Mech. 2014, 757, 701–746. [Google Scholar] [CrossRef] [Green Version]
  8. Kakac, S.; Pramuanjaroenkij, A. Review of convective heat transfer enhancement with nanofluids. Int. J. Heat Mass Transf. 2009, 52, 3187–3196. [Google Scholar] [CrossRef]
  9. Tsimring, L.S. Noise in biology. Rep. Prog. Phys. 2014, 77, 026601. [Google Scholar] [CrossRef] [Green Version]
  10. Deco, G.; Kringelbach, M.L.; Jirsa, V.K.; Ritter, P. The dynamics of resting fluctuations in the brain: Metastability and its dynamical cortical core. Sci. Rep. 2017, 7, 3095. [Google Scholar] [CrossRef] [Green Version]
  11. Damon, P.E.; Lerman, J.C.; Long, A. Temporal fluctuations of atmospheric 14C: Causal factors and implications. Ann. Rev. Earth Planet Sci. 1978, 6, 457–494. [Google Scholar] [CrossRef]
  12. Tyler, R.H.; Boyer, T.P.; Minami, T.; Zweng, M.M.; Reagan, J.R. Electrical conductivity of the global ocean. Earth Planet Space 2017, 69, 156. [Google Scholar] [CrossRef] [Green Version]
  13. Mobilia, M.; Georgiev, I.T.; Täuber, U.C. Phase Transitions and Spatio-Temporal Fluctuations in Stochastic Lattice Lotka–Volterra Models. J. Stat. Phys. 2007, 128, 447–483. [Google Scholar] [CrossRef] [Green Version]
  14. Brennecke, F.; Mottl, R.; Baumann, K.; Landig, R.; Donner, T.; Esslinger, T. Real-time observation of fluctuations at the driven-dissipative Dicke phase transition. Proc. Natl. Acad. Sci. USA 2013, 110, 11763–11767. [Google Scholar] [CrossRef] [Green Version]
  15. Das, S.; Kumar, K. Thermal flux in unsteady Rayleigh-Bénard magnetoconvection. Int. J. Heat Mass Transf. 2019, 142, 118413. [Google Scholar] [CrossRef] [Green Version]
  16. Pharasi, H.K.; Kumar, K.; Bhattacharjee, J.K. Frequency spectra of turbulent thermal convection with uniform rotation. Phys. Rev. E 2014, 90, 041004. [Google Scholar] [CrossRef] [Green Version]
  17. Burr, U.; Müller, U. Rayleigh–Bénard convection in liquid metal layers under the influence of a vertical magnetic fieldn. Phys. Fluids 2001, 13, 3247–3257. [Google Scholar] [CrossRef]
  18. Chandrasekhar, S. Hydrodynamic and Hydromagnetic Stability; Oxford University Press: London, UK, 1961. [Google Scholar]
  19. Fauve, S.; Laroche, C.I.; Libchaber, A.; Perrin, B. Chaotic Phases and Magnetic Order in a Convective Fluid. Phys. Rev. Lett. 1984, 52, 1774–1777. [Google Scholar] [CrossRef]
  20. Weiss, N.O.; Proctor, M.R.E. Magnetoconvection; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  21. Basak, A.; Raveendran, R.; Kumar, K. Rayleigh-Bénard convection with uniform vertical magnetic field. Phys. Rev. E 2014, 90, 033002. [Google Scholar] [CrossRef] [Green Version]
  22. Maxwell, J. A Treatise on Electricity and Magnetism; Oxford University Press: London, UK, 1873. [Google Scholar]
  23. Selimefendigil, F.; Öztop, H.F. Numerical study of MHD mixed convection in a nanofluid filled lid driven square enclosure with a rotating cylinder. Int. J. Heat Mass Transf. 2014, 78, 741–754. [Google Scholar] [CrossRef]
  24. Brinkman, H. The viscosity of concentrated suspensions and solutions. J. Chem. Phys. 1952, 20, 571–581. [Google Scholar] [CrossRef]
  25. Canuto, C. Spectral Methods in Fluid Dynamics: Springer Series in Computational Physics; Springer: Berlin/Heidelberg, Germany, 1988. [Google Scholar]
  26. Olson, P.; Glatzmaier, G.A. Magnetoconvection and thermal coupling of the Earth’s core and mantle. Phil. Trans. R. Soc. Lond. A 1996, 354, 1413–1424. [Google Scholar]
  27. Lan, C.W.; Kou, S. Heat transfer, fluid flow and interface shapes in floating-zone crystal growth. J. Crystal Growth 1991, 108, 351–366. [Google Scholar] [CrossRef]
  28. Bak, P.; Tang, C.; Wiesenfeld, K. Self-organized criticality. Phys. Rev. A 1988, 38, 364–374. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Plot of the total temperature field [ T ( x , y , z , t ) ] and velocity field [ v ( x , y , z , t ) ] at dimensionless time (a) t = t 0 and (b) t = t 0 + 5 in the simulation box for a water-based nanofluid [ Ra n f = 5.0 × 10 5 , Pr n f = 6.4 and Q n f = 100 ]. Colour bars describe the temperature distributions in the convecting fluid. Arrows show the directions of the fluid flow.
Figure 1. Plot of the total temperature field [ T ( x , y , z , t ) ] and velocity field [ v ( x , y , z , t ) ] at dimensionless time (a) t = t 0 and (b) t = t 0 + 5 in the simulation box for a water-based nanofluid [ Ra n f = 5.0 × 10 5 , Pr n f = 6.4 and Q n f = 100 ]. Colour bars describe the temperature distributions in the convecting fluid. Arrows show the directions of the fluid flow.
Fluids 06 00163 g001
Figure 2. Plot of the total temperature field [ T ( x , y , z , t ) ] and velocity field [ v ( x , y , z , t ) ] for different values of Chandrasekhar’s number: (a) Q n f = 100 and (b) Q n f = 250 in the simulation box for a water-based nanofluid [ Ra n f = 5.0 × 10 5 and Pr n f = 6.4 . Colour bars and arrows describe the distribution of total temperature field ( T ˜ ) and fluid velocity ( v ), respectively, at a given time t.
Figure 2. Plot of the total temperature field [ T ( x , y , z , t ) ] and velocity field [ v ( x , y , z , t ) ] for different values of Chandrasekhar’s number: (a) Q n f = 100 and (b) Q n f = 250 in the simulation box for a water-based nanofluid [ Ra n f = 5.0 × 10 5 and Pr n f = 6.4 . Colour bars and arrows describe the distribution of total temperature field ( T ˜ ) and fluid velocity ( v ), respectively, at a given time t.
Fluids 06 00163 g002
Figure 3. Temporal variations of the kinetic energy (E), entropy ( E Θ ) and Nusselt number ( Nu ) in a water-based nanofluid of Prandtl number Pr n f = 4.0 for the Rayleigh number Ra n f = 5.0 × 10 5 . The red curves are for the Chandrasekhar number Q n f = 100 , while the blue curves are for Q n f = 400 .
Figure 3. Temporal variations of the kinetic energy (E), entropy ( E Θ ) and Nusselt number ( Nu ) in a water-based nanofluid of Prandtl number Pr n f = 4.0 for the Rayleigh number Ra n f = 5.0 × 10 5 . The red curves are for the Chandrasekhar number Q n f = 100 , while the blue curves are for Q n f = 400 .
Fluids 06 00163 g003
Figure 4. Frequency power spectral densities (PSD) of (a) the energy per unit mass [ E ( f ) = | v ( f ) | 2 ], (b) the convective entropy per unit mass [ E Θ ( f ) = | θ ( f ) | 2 ], and (c) thermal flux [ Nu ( f ) ] in the frequency space for Earth’s liquid outer core ( Pr 0.1 , 0.2 ) and for water-based nanofluids with less than 8 % of spherical copper nanoparticles ( Pr n f 4.0 , 6.4 ) of for different values of Ra , Q , and Pr .
Figure 4. Frequency power spectral densities (PSD) of (a) the energy per unit mass [ E ( f ) = | v ( f ) | 2 ], (b) the convective entropy per unit mass [ E Θ ( f ) = | θ ( f ) | 2 ], and (c) thermal flux [ Nu ( f ) ] in the frequency space for Earth’s liquid outer core ( Pr 0.1 , 0.2 ) and for water-based nanofluids with less than 8 % of spherical copper nanoparticles ( Pr n f 4.0 , 6.4 ) of for different values of Ra , Q , and Pr .
Fluids 06 00163 g004
Figure 5. Variation of critical values of the dimensionless frequencies: (a) f c ( E ) , (b) f c ( E Θ ) , and (c) f c ( Nu ) for the frequency spectra of kinetic energy [ E ( f ) ], entropy spectra [ E Θ ( f ) ], and thermal flux [ Nu ( f ) ], respectively, with the Chandrasekhar number [ Q ] for Prandtl number ( Pr ) = 0.1 [red triangles] and 1.0 [blue circles].
Figure 5. Variation of critical values of the dimensionless frequencies: (a) f c ( E ) , (b) f c ( E Θ ) , and (c) f c ( Nu ) for the frequency spectra of kinetic energy [ E ( f ) ], entropy spectra [ E Θ ( f ) ], and thermal flux [ Nu ( f ) ], respectively, with the Chandrasekhar number [ Q ] for Prandtl number ( Pr ) = 0.1 [red triangles] and 1.0 [blue circles].
Fluids 06 00163 g005
Table 1. List of Prandtl number Pr , Chandrasekhar number Q , Rayleigh number Ra , exponents of Kinetic energy ( α ) , exponents of Entropy ( β ) , and exponents of Nusselt number ( γ ) .
Table 1. List of Prandtl number Pr , Chandrasekhar number Q , Rayleigh number Ra , exponents of Kinetic energy ( α ) , exponents of Entropy ( β ) , and exponents of Nusselt number ( γ ) .
Pr Ra Q Exponent α Exponent β Exponent γ
0.1 7.0 × 10 4 100 1.97 1.97 1.96
300 1.97 1.97 1.97
500 1.96 1.96 1.97
700 1.96 1.97 1.97
0.2 7.0 × 10 4 100 1.96 1.97 1.96
300 1.97 1.97 1.96
500 1.96 1.96 1.96
1.0 3.04 × 10 6 300 1.96 1.97 1.96
500 1.96 1.97 1.96
700 1.97 1.96 1.96
1000 1.96 1.97 1.97
2.0 3.04 × 10 6 500 1.96 1.96 1.96
700 1.96 1.96 1.96
1000 1.96 1.97 1.96
4.0 5.0 × 10 5 100 1.96 1.97 1.96
200 1.96 1.97 1.97
400 1.96 1.97 1.97
6.4 5.0 × 10 5 50 1.96 1.97 1.96
100 1.97 1.97 1.96
250 1.97 1.96 1.97
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Das, S.; Kumar, K. Frequency Power Spectra of Global Quantities in Unsteady Magnetoconvection. Fluids 2021, 6, 163. https://doi.org/10.3390/fluids6040163

AMA Style

Das S, Kumar K. Frequency Power Spectra of Global Quantities in Unsteady Magnetoconvection. Fluids. 2021; 6(4):163. https://doi.org/10.3390/fluids6040163

Chicago/Turabian Style

Das, Sandip, and Krishna Kumar. 2021. "Frequency Power Spectra of Global Quantities in Unsteady Magnetoconvection" Fluids 6, no. 4: 163. https://doi.org/10.3390/fluids6040163

APA Style

Das, S., & Kumar, K. (2021). Frequency Power Spectra of Global Quantities in Unsteady Magnetoconvection. Fluids, 6(4), 163. https://doi.org/10.3390/fluids6040163

Article Metrics

Back to TopTop