1. Introduction
The instability of a uniform horizontal porous layer that is heated from below is the direct analogue of the Rayleigh-Bénard problem and is known as the Darcy-Bénard problem, the porous Bénard problem, or the Horton–Rogers–Lapwood problem. The first works on this linear instability problem were by Horton and Rogers [
1], and Lapwood [
2]. Since then researchers have extended these pioneering works, not only into the nonlinear regime, but also by adding further complications such as modified boundary conditions, internal heating, the presence of local thermal nonequilibrium, more than one diffusing species, non-Newtonian fluids, rotation, and the like. Given that the monograph by Nield and Bejan [
3] devotes more than 100 pages to this, we shall refer the reader to their review of the state of the art.
The topic of the present work restricts the porous layer (which is often assumed to be unbounded horizontally) to that of a finite cavity, further restricts it to admitting solely two-dimensional flow, but allows the cavity to be inclined. Once more, Nield and Bejan [
3] devote more than six pages to a review of the general topic of inclined layers, and once more, we refer the reader to their monograph. Rather, we shall review briefly those papers which are of very direct relevance to our work.
The first paper which considered the onset of two-dimensional convection in an inclined porous layer and which presented an exhaustive account of the linearized stability theory was by Rees and Bassom [
4]. Using a finite difference method to convert the linear stability equations into a matrix eigenvalue, problem they obtained the neutral curves for various inclinations of the layer. They found that there is a maximum inclination above which the basic state is linearly stable. Thus, we may expect convection to arise in an unbounded inclined layer when the inclination satisfies
, but the basic state is otherwise linearly stable.
A much more recent work by Wen and Chini [
5] considered nonlinear convection at moderate values of the Darcy–Rayleigh number and presented solutions when the inclination is as high as
, which is above the linear threshold given in [
4]. Further work by Wen and Chini [
6], which concentrated mainly on convection at very high values of the Darcy–Rayleigh number, used an energy method to show that nonlinear convection can occur at all inclinations as long as the Darcy–Rayleigh number exceeds
. Clearly, the source of these nonlinear solutions owes little to linearized stability theory, and this motivated Rees and Barletta [
7] to investigate how nonlinear solutions may be generated as the inclination increases past
. It was suspected that it could be explained only by the presence of subcritical instabilities arising somewhere within the range of
, and therefore, a weakly nonlinear theory was used to look for clues. It was found that the onset of convection is supercritical when
, but is otherwise subcritical. The analysis was confirmed by some strongly nonlinear finite difference computations for
, where for any chosen value of
within that range, strongly nonlinear flows were found at subcritical values of the Darcy–Rayleigh number. Further work on this aspect of the layer analogue of the present cavity problem is in progress.
The works of Wen and Chini [
5,
6] were concerned with porous layers. At roughly the same time, Guerrero-Martínez et al. [
8] used a finite difference method to determine the effect of inclination on the convective flows within a porous cavity as a function of the Darcy–Rayleigh number, the inclination, and the aspect ratio. They too discovered strongly nonlinear flows at
, and this fact alone forms the main motivation for the present study. So, one possible question to be answered is: is there a maximum inclination for which nonlinear convection exists, by which is meant something other than a buoyancy-induced single cell state? The paper by Guerrero-Martínez et al. [
8] presents results in the form of the variation of the Nusselt number with the inclination, but the solution curves often just terminate. This phenomenon tends to be indicative of the unseen presence of unstable steady states that cannot be computed using unsteady methods. It is also indicative of fold bifurcations. Therefore, we have selected to use a spectral method, which reduces the governing steady-state equations to a system of nonlinear algebraic equations which are then solved using the Newton-Raphson method. The presence of fold bifurcations also dictated the need to use a pseudo-arclength continuation scheme.
It is important to note that a similar type of analysis was undertaken by Riley and Winters [
9] using a continuation scheme, but where the basic solver involved a finite element discretization. They confined their analysis to a tilted unit cavity, and displayed their results in terms of increasing values of the Darcy–Rayleigh number.
In the following section, we present the governing equations and then we describe the numerical scheme in
Section 3.
Section 4 contains the result of our numerical validation exercise, together with a comparison with one case shown in [
8] in order to demonstrate the utility of the present numerical scheme. Detailed remarks are also given about modal exchange mechanisms in terms of the direction of circulation of the convection cells.
Section 5 forms our main results section where we show various examples of modal interaction for different cavity aspect ratios and values of Darcy–Rayleigh numbers. The general aim in that section is to acquire some intuitive understanding of what interactions are likely to happen in general. Finally, we present some conclusions in
Section 6.
2. Governing Equations
We will consider nonlinear two-dimensional convection in a porous cavity that has aspect ratio,
L, and is inclined at the angle
to the horizontal. It is assumed that the solid matrix and the fluid are in local thermal equilibrium, that Darcy’s law and the Boussinesq approximation apply, and that the fluid density is a linear function of the temperature. Subject to these assumptions, the dimensional governing equations are
and
The coordinates are
and
, which are, respectively, the coordinate up the layer and the coordinate that is perpendicular to the bounding surfaces, while the corresponding Darcy velocities are
and
. The temperature and the dynamic pressure are
T and
. Here,
K is the permeability,
the dynamic viscosity,
the coefficient of cubical expansion,
g the acceleration due to gravity,
the ratio of the heat capacity of the saturated porous medium to that of the saturating fluid,
the reference temperature, and
the fluid density at the reference temperature. Finally,
is the ratio of the average thermal conductivity of the saturated medium and the product of the density and the specific heat of the fluid.
The boundary conditions are as follows,
and
where
h is the thickness of the porous cavity,
L is its aspect ratio, and
. We shall restrict our attention to two-dimensional steady flows, and therefore, we shall set
and eliminate all
- and
-derivatives.
The full system of equations and boundary conditions may now be rendered dimensionless by using the following scalings,
The governing equations and boundary conditions become,
and
In the above, the Darcy–Rayleigh number,
, is given by,
Given that the flow is now taken to be two-dimensional, we may define the streamfunction
, using
which automatically satisfies the equation of continuity, Equation (
8), and hence, we obtain the system
Figure 1 shows the precise boundary conditions for
and
that need to be satisfied. In addition, the figure also depicts the isotherms for a weak free-convective circulation where the hot lower surface induces a buoyant flow up that surface, and the cold upper surface causes a corresponding flow down that surface.
4. Comparison with Guerrero-Martínez et al. [8] and Some Flow Patterns
It also seemed apposite to undertake a final check of the present encodings by comparing our results with those of Guerrero-Martínez et al. [
8] in order to show the utility of the pseudo-arclength continuation in the present context. In what follows, the value
n corresponds to the number of distinct cells within the cavity. Our mathematical definition of
n is one more than the number of sign changes in
at
within the range of
.
Figure 2 displays the variation of the Nusselt number with inclination for the case of
and
.
Figure 2a shows the finite difference computations of [
8] where their definition of the Nusselt number, which is
, varies with
. The separate (and differently coloured) curves correspond to different numbers of cells within the cavity (
,
and
), while the vertical dotted lines correspond to when their converged solution for one value of
has jumped suddenly to a different convection pattern as
has been increased or decreased incrementally.
Figure 2b shows the solution curves that were computed using the present spectral method with pseudo-arclength continuation. Those curves that are shown in
Figure 2a are reproduced in
Figure 2b with no discernible difference, and therefore, our computations agree with those in [
8]. However, it is evident that
Figure 2b also displays a large number of extra curves, all of which are fully converged steady-state solutions of the governing equations. These extra curves are either very difficult to be computed with an unsteady solver, or (which is much more likely to true) they represent solutions that are unstable.
The uppermost curve corresponds to a three-cell convection pattern which appears to approach a fold bifurcation as
increases because the slope of the Nusselt number curve becomes very large. In
Figure 2a, the solution then drops down to the red curve, which is a one-cell pattern. In
Figure 2b the corresponding curve turns back towards smaller inclinations. Finally, the curve undergoes a second fold bifurcation to the one-cell pattern that is characteristic of flows at high inclinations.
For this choice of parameters,
Figure 3 shows the streamlines and isotherms corresponding to the four representative examples that correspond to the black disks in
Figure 2b.
Figure 3a,b represent three-cell flows for
and
. In both cases, the outer cells are stronger than the middle one because their anticlockwise circulations are aided by the effect of buoyancy forces along the boundaries; we will call these natural circulations. By contrast, the inner cell is weaker because its circulation is inhibited by the action of buoyancy forces along the boundaries; this is an anti-natural circulation. This increasing inhibition of the middle cell as
increases not only leads to weaker circulations, but also to a decreased width.
As the curve is traversed, it passes a fold bifurcation when descends, transitions from representing a three-cell flow to representing a one-cell flow (see the change in colour), and encounters a second fold bifurcation at . Subsequently, the curve continues to represent a single cell until and beyond when the layer becomes vertical.
The process of transition from a three-cell state to a one-cell state may be seen by comparing
Figure 3b–d. The middle cell seen in
Figure 3b is then squeezed out of existence, leaving behind the last vestiges of the naturally circulating cells, as seen in
Figure 3c. In this case, the streamfunction is single-signed within the computational domain, and therefore, this is regarded as a single cell. The precise point along the
-curve where three cells become one is again marked by the change in the colour. Finally, at larger inclinations, all of the fluid circulates about the center of the cavity, as shown in
Figure 3d.
This type of transition between patterns with odd values of n is quite ubiquitous, and it is also possible to have a similar transition between two even values of n.
Figure 4 shows the two
cases, which are shown as yellow disks in
Figure 2. The inclination is much shallower than those used in
Figure 3, and therefore, the small difference between the streamlines may only be seen upon close inspection, but the underlying isotherm patterns are very different from one another. In
Figure 4a, the outer cells display natural circulations, while in
Figure 4b, they display anti-natural circulations. Therefore, it is to be expected that the former flow has a larger Nusselt number than the latter.
Figure 5 is an example of a transition between a pattern with an odd number of cells and a pattern with an even number. Again, the parameters are the same:
,
, and
, but they correspond to the green disks in
Figure 2. When one has four cells then one of the outer cells will be anti-natural; in
Figure 5a, which corresponds to the uppermost green disk in
Figure 2, it is the right-hand cell that plays this role. Indeed, it is quite clearly narrower and weaker than the remaining cells.
As the
curve is traversed towards where the middle green disk is, the right-hand cell continues to weaken and it is now barely present in
Figure 5b.
When one continues along the
curve, one eventually passes through
and towards negative values. This will bring one to the location of
, which is the mirror image of the lowest green disk in
Figure 2c. The streamline and isotherm pattern for that point will be precisely the reflection of
Figure 5c about
. However, the pattern we see represents
, and this point is very close to where the
solution bifurcates from the middle branch of the
(blue) curve. The bifurcation point itself marks when the former
solution has finally removed the right-hand anti-natural circulation. This evolution provides a mechanism whereby one cell is either lost or gained at a bifurcation point.
5. Results Using the Spectral Method
It is clearly impossible to present a comprehensive account of how the pattern of convection changes as a function of
,
L, and
, despite this being just a three-parameter set. Even with the results shown in
Figure 2, it is also clear that there exist many different solutions corresponding to different numbers of cells, even when two of the three parameters have been fixed. That said, it is possible to be comprehensive when
L is relatively small and
is not too large. Therefore, we shall present a selection of different cases with the overall aim of attempting to acquire a good intuition about how this cavity problem behaves in general. We shall split the presentation into subsections where the different aspect ratios are considered in each.
Before that, though, it is good to recall critical values of
as a function of
L for when the cavity is horizontal. It is well-known that the critical Rayleigh number is
for the horizontal porous layer [
10], where
k is the wave number. For an unbounded layer,
k may take any real value, but for a cavity of width
L, the wavenumber takes the discrete values
, where
n is the number of cells within the cavity. Therefore, Equation (
28) translates into
In all of the subsequent figures, the solution curves will be coloured to indicate the parity of the flow, i.e., the number of cells that exist within the cavity:
So, modes with an odd parity have been assigned rustic colours and modes that have an even parity when
are formed as dashed lines. We note that
Figure 2 in this paper has adopted the colouring scheme used in [
8].
5.1. Solution Curves for
This is the unit aspect ratio cavity. Given that the critical wave number for the unbounded horizontal layer is
, which corresponds to a length of 2 or, equivalently, a unit cell width, we would expect a single-cell flow to dominate the nonlinear dynamics for the unit cavity. If we define
to be the critical Darcy–Rayleigh number for an
n-cell flow when the cavity is horizontal, then we may use Equation (
29) to provide the following values:
where the numerical values are correct to three decimal places. These form the context for the following figures.
Figure 6 shows the evolution as
increases in the single-cell
curves through a whole rotation of the cavity through 360 degrees. When
, the heating is from above and there is no flow at all, and hence,
which corresponds to pure conduction. At other inclinations, a flow is generated and this provides an increased Nusselt number.
When
increases from zero, a relatively weak buoyancy drives a slow flow and the Nusselt number increases. Eventually, the Nusselt number reaches a maximum and decreases back to
when
. All of the curves shown in
Figure 6a correspond to natural circulations. However, once
exceeds
, there exists the possibility of having a thermoconvective instability and hence, a strong convection when
. We see this in
Figure 6b, where each of the solution curves exhibits a loop. It is perhaps easiest to focus on the
curve, where as
decreases from
, we eventually obtain a non-unit value of
when
, and this is a strongly nonlinear solution. Within this range of inclinations, the convection is a natural circulation, but once
becomes negative along the same curve, then the corresponding convection cell consists of an anti-natural circulation where buoyancy forces along the the boundaries again serve to inhibit the convective instability. Eventually, the curve passes through
where
and there is no flow. Everything that happens thereafter is the mirror image of what has already been described.
Figure 7 shows some solution curves where the values of
have been chosen to show the evolution of the system as
increases. When
, we see immediately that there is a short solution branch corresponding to two cells,
. This is consistent with the fact that
; see Equation (
31). This particular value of
is not strongly supercritical, and therefore, the
flow is weak and
is only slightly larger than 1. The
branch terminates at a bifurcation from the
curve, which is of the same nature as was shown in
Figure 5; in this case, the two-cell flow that exists at
gradually loses the cell with the anti-natural circulation, and hence, an
state is eventually obtained as
increases. We also note that the red
curve is symmetric about the vertical axis.
As
increases further, the
curve develops a cusp at
(see
) and then a loop (see
and beyond). In the meantime, a small section of curve close to
and
is being identified as being an
solution and is coloured in green. At
, a cusp may then be seen when
, and this marks the onset of the convection of a three-cell mode, since
is now almost exactly the value of
; see Equation (
31). When
, the
cusp has now developed into a loop, reflecting the fact that we now have a three-cell solution of moderate amplitude.
At , the curve has become more complicated, while the Nusselt number for has increased to such an extent that it is only just smaller than that for the solution.
In all of the above, we have merely computed the solution curves, and the present code/method, despite its complexity, is incapable of providing information on the stability of the computed solutions. It is also worth mentioning that the solution curves become increasingly intricate and difficult to obtain as increases further. Thus far, only the and 3 solutions have appeared, but we would need to consider larger values of in order to see some solutions, since when .
5.2. Solution Curves for
For
, we may use Equation (
29) again to provide the following critical values for
for different values of
n:
Naturally, for a cavity of length, , we would expect the primary mode of instability to be formed from two cells. Clearly, an pattern is expected to appear at a slightly higher Darcy–Rayleigh number, and then both the and patterns arise subsequently at the same value of .
Figure 8 shows some subcritical cases for
, 20, 30, and 40. Qualitatively, these single-cell solutions have exactly the same variation with
as was seen in
Figure 6a for
.
Figure 9 shows the solution curves for six chosen values of
. When
, which is just above the critical value for the
pattern, we see an
curve within a small range of inclinations centered on
. This weak nonlinear flow bifurcates from the unicellular state when
, and the
curve then forms the unique solution for larger inclinations.
At the slightly higher , which is above the critical value of for and , we see that the solution curves have already evolved quite substantially. The black curve that is present for large inclinations eventually transforms into the green labeling as decreases, and a loop is formed near to . So, when the layer is horizontal, there are only two nonlinear solutions, one with two cells and one with three. We also note that the red curve now bifurcates from a green curve.
As
increases still further, the evolution of the solution curves becomes increasingly complicated. Thus, as
increases from 60 through
, the black
curve centered on
develops a cusp that then transforms into a loop with a nonlinear solution at
. At the same values,
, a blue
solution first appears, and a weak nonlinear solution may be seen when
. These curves continue to evolve and develop further loops as
continues to increase, and when
, we also see the appearance of the orange
solution, given that
is now above the value of
given in Equation (
32).
5.3. Solution Curves for
Figure 10 is the
counterpart to
Figure 6a and
Figure 8, and it has been plotted using the same vertical scale for the sake of comparison. The shapes of these curves offer no surprises, except that the
curve is shaded in green when the layer is close to being horizontal. For such a cavity, the
mode has
as its critical Darcy–Rayleigh number when
, and therefore, the value of
is slightly supercritical. Clearly, the
solution evolves into a three-cell solution as
tends towards zero because of the presence of thermoconvective instability.
For
we may again use Equation (
29) to determine the following critical values of
when the layer is horizontal:
Figure 11 shows the detailed solution curves for
, 50, 70, and 100. The subfigure corresponding to
is a close-up of the uppermost curve in
Figure 10, but there is a very narrow green loop centered on
that shows that
is supercritical and that a nonlinear
solution exists, albeit a fairly weak one.
A small increase in
to 50 brings with it the potential for three patterns to exist:
,
, and
; see Equation (
33). However, first, we note the
curve, which is the sole solution for large inclinations, transforms into an
curve as
decreases, and as
decreases still further towards zero, a final transformation to a five-cell state takes place. The orange colour is a herald of the fact that the critical value,
, is only very slightly higher than 50.
We do indeed see the presence of the two-cell and four-cell solutions, as predicted, where the latter is stronger when
, and this is possibly due to the fact that it has the lower critical Darcy–Rayleigh number. Both the
and
solutions bifurcate from the central loop, which transitions between being
and an
curve. A more complicated scenario may be seen for
. The
and
solutions are now much stronger, and the
solution, hinted at when
, is also well into the nonlinear regime. It is of interest to note that the straightforward bifurcation of the
(blue) solution from the
(green) solution that we saw for
now assumes a more complicated form with the presence of a very tight loop. Indeed the presence of tight loops turns out to be quite ubiquitous when either
or
L increases. We also see the relatively weak presence of an
(purple) pattern which bifurcates from the
branch at the bottom of the subfigure; its presence somewhere was expected, given the value of
presented in Equation (
33).
When
, we obtain the case that we considered in
Figure 2. Given Equation (
33), we could have predicted the presence of modes
through to
, and the newly added modes
(brown) and
(grey) are quite weak. Mode
bifurcates from an
branch in a manner that is now quite familiar. On the other hand, the way in which the
mode appears is quite unusual. That said, one of main characteristics of this subfigure is the number of loops that occur. All three solutions with even parity exhibit a loop centered about
, while another arises close to
and
. There are also two very sharp changes in direction that are difficult to follow in the present code, despite it being designed to do so. These are at
and
.
5.4. The Effect of Varying L When
Thus far, we have chosen to consider three different aspect ratios of the cavity, and to concentrate on how the solution curves vary as
increases. All of these aspect ratios correspond to integer values of
L, and therefore, the onset of convection arises when
when the cavity is horizontal, and
forms the most unstable mode of convection. It is well-known, and easily derivable from Equation (
29), that the most unstable mode has
n cells when
L satisfies
Therefore, this subsection is devoted to illustrating how inclination and nonlinearity affect this conclusion when
L takes noninteger values. We shall allow
L to increase from
to
in steps of
. This range spans the full transition from an odd parity pattern through to an even parity, and on to the next odd parity pattern. We have chosen to consider just the one case,
, which represents a moderately supercritical Darcy–Rayleigh number. Nine different cases may be seen in
Figure 12, and all of the subfigures use the same vertical scale.
When
, we are well within the range of values of
L for which we expect one cell when
, but an
solution bifurcates from the
branch and is quite weak. The large-
solution curves are exactly like those seen in
Figure 6b, for example. When
, we are just below
, which is when the
and
modes have the same critical Darcy–Rayleigh number. Nevertheless, the two-cell solution is slightly stronger, at least in terms of the Nusselt number, over a narrow range of inclinations. A green colouring at the bottom of the subfigure corresponds once more to the imminent appearance of an
solution for which the critical Darcy–Rayleigh number is
.
The green pattern may be seen in the , , and subfigures, and these flows increase in strength as L increases within that range. However, the pattern now dominates all other modes over an increasing range of inclinations centered on . When this range is roughly . In all of the subfigures, up to and including , the solution itself bifurcates from the solution.
We note that the onset of an pattern when is when , which is less than 65, and indeed, a close examination of the lowest part of the curve (i.e., under quite strong magnification) shows the initial appearance of the blue solution. We have examined this in detail, and it displays the same qualitative transition via the red state to a bifurcation from the black as the subfigure shows the red curve transitions via the green state to a bifurcation from the (continuous) curve.
As
L increases further via
, we see a strengthening of the
flow relative to the
flow. The
pattern also continues to grow, which is opposite to that of the
flow. These blue curves are incomplete and we have been unable to compute them. However,
Figure 11 (
, a close case) shows what we believe to be the shape which will be adopted here.
Our conclusion based on this sample of computations is (i) that when L is close to being an odd integer, m, then the pattern dominates for quite a wide range of inclinations; and (ii) that when L is close to being an even integer, m, then the pattern dominates, but over a relatively narrow range of inclinations; that when L approaches an odd integer value, m, then the strength of the pattern weakens while that of the pattern increases (e.g., comparing the behaviours of the red and blue curves between and ). We presume that a similar observation may be made concerning when L approaches an even integer value.
5.5. Some Solutions for Yet Larger Aspect Ratios
It is perhaps no surprise that the computational difficulties increase when either or L are increased. Larger values of the Darcy–Rayleigh number allows for further modal patterns to fit within the cavity, and for an accompanying increase in the intricacy for how those different patterns interact. While the same observations are also true when L increases, not only does the number of computations increase, but the length of each computation itself increases as the third power of , and also the number of spectral modes used in the x-direction.
In this subsection, we shall confine ourselves to
and to the following cavity aspect ratios:
, 5, 6, and 7. The aim will be to determine a qualitative understanding of the modal exchanges that take place in a shallow cavity as the aspect ratio increases.
Figure 13 depicts the corresponding solution curves, although none of the individual subfigures is complete. Our initial numerical experimentation showed that an astonishing degree of complexity is obtained in the region fairly closely centered on
, and for smaller values of
. This has already been seen in Figure for
. Therefore,
Figure 13 has retained only those solution curves that exist for the larger inclinations (usually those with odd parity) and a selection of even parity solutions where
n and
L are close to one another.
If attention is confined to the even-parity modes shown in
Figure 13, then we see that there is a natural tendency for curves corresponding to a fixed value of
n to rise as
L increases towards the value of
n, and to descend as
L increases further. The presence of even parity cells is generally confined to within approximately
of the horizontal, although there is a tendency for that limit to increase slightly as
L increases. It is also interesting to note that the
mode transfers more heat than any other mode when
and
is close to zero, and the same may also be said for the
mode when
.
If we now confine attention to the larger inclinations where only those solutions with an odd parity exist, then we see the presence of many fold bifurcations. Generally, the right-most structure takes an S-shape and it connects a pair of fold bifurcations. Other structures at smaller inclinations take a distinctive downward-pointing tongue-like shape, which is illustrated well for the case.
Interactions and couplings between the solution curves as
L increases are many and of different kinds, and it is difficult to generalize or even to describe them easily. For example, the green
loop for
has become a cusp when
, and although we do not display it here, it becomes a smooth fold when
. Further, the
curve, which forms a unique nonlinear solution at high inclinations, eventually transforms into the green
pattern when
, but it transforms into the orange
pattern when
. In this regard, the transitional case,
, is not clear in
Figure 13, although the numerical data show that the transition from
is to
, as it is for
. We illustrate how this happens in
Figure 14 where we have used the aspect ratios,
and
. For both the cases presented, we note that each consists of just one continuous solution curve. The computation itself is displayed as a thick line, while its mirror image is displayed using a thin line in order to clarify what happens to the solution trajectory when
L increases from
to
. There is clearly more than one qualitatve difference between the two trajectories.
5.6. A Maximum Inclination for Convective Instability?
Finally, and given that one of the original motivations for this study was the observation by Guerrero-Martínez [
8] that buoyant instability is manifested in a porous cavity at an inclination that is well above that given by the linearized stability theory for the unconfined layer, it is worth presenting some initial data on this aspect.
The preceding figures demonstrate that when
or larger, then the right-most solution curve displaying the
to
modal transition is the one that provides the maximum inclination for the existence of buoyancy-induced instability. This inclination may be obtained by the simple expedient of fitting a suitable parabola to the ‘nose’ of such a curve and finding its extremum.
Table 2 shows the result of this process, and also of the second ‘nose’, the one with the
to
transition; these data are also displayed in
Figure 15.
A brief glance at
Figure 15a shows that there is an increasing freedom for strongly nonlinear convection to exist as the aspect ratio increases, but it also suggests that there could be a limiting value of the maximum inclination when
L is asymptotically large. When plotted against
in
Figure 15b, there is a suspicion that
could be a linear function of
, and that the limiting inclination when
could therefore be very close to
for the green disks. Intriguingly, the same limiting value appears to apply for the orange disks. At present, we do not know whether this is a significant observation or not, but we shall deem it to be worthy of further investigation later.
An alternative viewpoint is provided by
Figure 16, which shows solution curves for
and for values of
in steps of 10 between
and
. Here, we show only a section of the right-most solution curve in each case, i.e., that part which is close to the maximum inclination for cellular convection. The locus of the respective maximum inclinations are given by the dashed red line, and it appears that the variation of
is not too far from being linear. Clearly, it is impossible to predict whether the maximum inclination will continue to increase in this way as
increases further, or whether there will be an maximum inclination that is independent of
for a fixed value of
L, or whether there is a global maximum inclination that is independent of both
and
L. Again, we intend that this aspect will be the subject of subsequent research, but we can certainly conclude that nonlinear convection is possible when
is as large as
from the horizontal.
6. Conclusions
This has been a three-parameter system, something for which a fairly comprehensive set of results may usually be obtained, and from which a definitive understanding of the problem is often acquired. It is clear that these desirable outcomes are not possible for the inclined Darcy-Bénard cavity. Part of the reason for this is the possibility of having more than one potential convection pattern. Another reason is the potential for different solution trajectories to interact and to change their morphologies as a governing parameter changes. We have uncovered certain principles for these changes and some features that are ubiquitous:
Even parity solutions have a relationship that is even about .
Odd parity solutions consist of two different sets where one set consists of flows with one more cell having a natural circulation than having an anti-natural circulation. The other set has one more cell with anti-natural circulations. The former persists to larger inclinations and transfers more heat.
Even parity solutions can become of odd parity by either losing or gaining a cell at a bifurcation.
Odd parity solutions can transform into other patterns with odd parity by losing or gaining an even number of cells. The same mechanism applies when even parity solutions transform gradually into other even parity solutions.
At sufficiently high inclinations, the flow is unique and consists of just one circulating cell that is driven by the natural buoyancy forces on the hot and cold boundaries.
When the aspect ratio is sufficiently large, then S-shaped solution trajectories are obtained, as are tongue-like shapes. These have been found to interact as the aspect ratio changes.
For a fixed value of , the evidence so far suggests that there is a maximum inclination angle above which no nonlinear convection arises, apart from the single-cell circulation.
For a fixed aspect ratio, the evidence so far shows that it is possible to have strong convection when is larger than , and this suggests strongly that larger inclinations may be achieved.
Whilst we have presented a carefully curated and very large set of results involving many hundreds of hours of computation, there remains much more that could be done. One possible desire would be the creation of an animation which marches slowly through the equivalents of
Figure 12 and
Figure 13 as
L increases in small increments. This would give a sense of understanding that a few printed figures do not. The same would be true for a fixed value of
L, but where
increases slowly, and we would then be able to witness how the flows might evolve as the physically variable parameter,
, varies. We would anticipate that hysteretical effects will arise. Other possible questions include those of distinguishing between solutions that are stable (i.e., computable with a time-dependent solver) and those that are not. Another one would query how precisely does the maximum inclination for strongly nonlinear convections vary with both
and
L, and can the large-
L limit be subject to some type of combined analytical/numerical treatment? Additionally, whilst the present numerical scheme has allowed us to follow the twists and turns of the solution trajectories, it has not been sufficiently able to follow every curve easily, particularly when the change in direction is very rapid, or when there is a neighbouring but distinct trajectory that will sometimes cause the Newton–Raphson scheme to slow down considerably or not to converge at all.