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
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 (
), which is a measure of heat flux across the fluid layer in RBC, varies with frequency (
f) as
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 (
) also showed the similar behaviour in the presence of Lorentz force [
15] in water-based nanofluids (
). The power spectral density (PSD) of the thermal flux in the frequency (
f) space [
2,
5,
15,
16] was found to vary as
. The frequency spectrum for the temperature field measured near the horizontal plates in Rayleigh–Bénard magnetoconvection [
17] shows the exponent to vary between
and
. 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 (
), and Nusselt number (
) in unsteady Rayleigh–Bénard magnetoconvection (RBM) [
18,
19,
20] for several values of Rayleigh number (
), Prandtl number (
), and Chandrasekhar number (
). 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
at relatively higher frequencies. This behaviour does not depend on the Rayleigh number (
), Prandtl number (
), and Chandrasekhar’s number (
).
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.,
melt), nanofluids, etc.) of density
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
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 [
] and the fluid density [
], in the conduction state [
18], are given as:
where
and
are the temperature and density of the fluid at the lower plate, respectively. The fluid pressure [
], in conductive state, is:
where
g is the acceleration due to gravity and
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 [
] is known. If we take
, where
is a constant (e.g., air pressure at the upper plate), the constant
C turns out to be
The fluid pressure, in the basic conductive state, then reads as:
As soon as the temperature gradient across the fluid layer is raised above a critical value
for fixed values of all fluid parameters (kinematic viscosity
, thermal diffusivity
, and thermal expansion coefficient
) and the externally imposed magnetic field (
), the convection sets in. All the fields are perturbed due to magnetoconvection and they may be expressed as:
where
,
,
and
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
. The convective temperature field
and induced magnetic field
are made dimensionless by
and
, respectively. The magnetoconvective dynamics is then described by the following dimensionless equations:
where
is the material derivative. The unit vector
is directed vertically upward. Because the magnetic Prandtl number is very small (
) for all terrestrial fluids, we set
equal to zero in the above set of hydrodynamic equations. The Navier–Stokes equation Equation (
10) then takes the form, as given below.
The induced magnetic field (
) is slaved to the velocity field (
), as Equation (
11) is simplified to
The fluid flow due to magnetoconvection, in the limit of
, 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
and
, are:
All of the fields are considered periodic in the horizontal plane. The fluid dynamics, as
, is controlled by three dimensionless parameters: (1) Rayleigh number (
), (2) Prandtl number (
), and (3) Chandrasekhar’s number (
). The critical values of Rayleigh number [
] and critical wave number [
] are [
18]:
The kinetic energy (
E) and convective entropy (
) per unit mass are defined as:
The Nusselt number (
), which is the ratio of total heat flux and the conductive heat flux across the fluid layer, is given as:
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:
where
and
are the density and electrical conductivity of the base fluid, respectively. Here,
is the density and
is the electrical conductivity of nanoparticles. The effective thermal conductivity (
) of a nanofluid [
22] for smaller values of
is expressed as:
where
and
are the thermal conductivities of the base fluid and that of sphere shaped nanoparticles, respectively. Similarly, the product of density and specific heat capacity [
] may be expressed through the following relation [
23]:
The effective dynamic viscosity (
) of a nano-fluid [
24] may also be expressed as:
The kinematic viscosity (
) and thermal diffusivity (
) for nanofluids may then be estimated while using the formulae:
The effective thermal Prandtl number (
) for nanofluids is then defined as
. The definitions of Rayleigh number (
) and Chandrasekhar’s number (
) are also modified for nanofluids, and they read as:
Substituting the values of dimensionless parameters
,
and
by
,
and
, respectively, the set of Equations (
12)–(
15) describes magnetoconvection in nanofluids with low concentration of nanoparticles in a base fluid. The thermal Prandtl number (
) of water-based nanofluids may be varied between
and
by varying the volume fraction (
) of the spherical copper nanoparticles [
15] between
to
.
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:
where
and
. The time dependent Fourier amplitudes of these fields are denoted by
and
, where
l,
m, and
n are integers. The horizontal wave vector of all the perturbations is
, where
and
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
, where
. The value of
is computed using the expression for the critical wave number
[see Equations (
18) and (
19)]. The continuity equations decide the possible values of the integers
. They can take values that satisfy the following equation.
The computation of non-linear terms and are computed while using Fast Fourier Transformation (FFT). It is done using the following steps:
- (i)
Real space variables and are computed at a given time t using inverse FFT of and , where
- (ii)
The multiplication of field variables and for () are done at each grid point of the simulation box.
- (iii)
FFT[] and FFT[] are computed using the package FFTW.
- (iv)
Subsequently, the terms and with as well as and with are computed.
The aliasing error is removed using
-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
(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
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
, we carried the simulations with higher resolutions of
. 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
,
and
is completed, the same procedure is repeated with a new set of
,
, and
.
4. Results and Discussion
The simulations are done for several values of thermal Prandtl number (
). These values of
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
, while the Chandrasekhar’s number is varied in the range
. In numerical simulations, the value of
is always mapped to unity. The total temperature field [
], in the dimensionless units, then reads as:
.
Figure 1 shows the combined plot of the total temperature field [
] and the velocity field [
] simultaneously at two different instants for a water based nanofluids (
) for
and
. The dimensions of the simulation box are
, where
for
. 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
and
and for two different values of Chandrasekhar number (
).
Figure 2a shows the total temperature field [
] and the velocity field [
] in the simulation at a randomly chosen time for
. The flow structure that is shown in
Figure 2b is for
.
Figure 3 shows the temporal variations of three global quantities for
,
and for two different values of
: (1) the kinetic energy per unit mass (
E), (2) the convective entropy per unit mass (
), and (3) the Nusselt number (
). 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
. The red curve is for
and the blue curve is for
. The energy signal varies irregularly with time. They show appreciable fluctuations. The mean of kinetic energy decreases with increase in
. The fluctuations of kinetic energy also decreases with an increase in
. Curves in the third and fourth rows (from the top) show the temporal variation of
, while curves in the fifth and sixth rows display the temporal signal of
. 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
. Their fluctuations also decrease, as
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
,
and
. The PSDs of the fluid speed [
] 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,
), the slope of curves
on the log-log scale varies between
to
. The energy spectra [
] 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 (
and
) 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
is found to have insignificant noise for
. The frequency spectrum of energy per unit mass [
] 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 [
] of the kinetic energy scales with frequency (
f) almost as
with
. The scaling behaviour is valid for more than two decades. The PSD shows a clear scaling behaviour for
. The scaling exponent is independent of
,
and
in this frequency window.
Table 1 gives the numerically computed values of the exponent (
) for different values of
,
, and
. A similar scaling law [
] was also observed in rotating Rayleigh–Bénard convection [
16].
Figure 4b shows the PSDs of the convective entropy [
] of the fluid in the frequency space for different values of
,
and
. Its power spectra is also noisy in the dimensionless frequency range
. The slope on the log-log scale varies between
and
. However, for
,
scales with frequency as
with
.
Table 1 lists the numerically computed values of the exponent
. Interestingly, the power spectra of the temperature fluctuations are also found to vary as
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
for
. Even for very small values of
(
), 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,
for several values of
,
, and
. 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
, if
. The scaling exponent varies between
to
in this frequency window. However, for higher values of dimensionless frequency (
), the spectra for thermal flux [
] also shows very clear scaling:
, with
.
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
, which is different than the scaling behaviour in the hydrodynamic systems discussed here.
The scaling law showing the variation of the power spectra as
starts at a critical frequency (
) for different values of
.
Figure 5 shows the variation of
for
,
, and
with
for two different sets of
and
. The critical frequency
decreases, as
is increased (see
Figure 5a).
Figure 5b,c show the variations of
and
, respectively, with
. The values of critical frequencies are different for
,
and
. However they all decrease with an increase in
for given values of
and
. 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
, as
is raised.