Abstract:Thermal convection in a horizontal layer with uniform internal heating and stress-free constant-temperature boundaries is analyzed numerically. The work is motivated by the questions arising in development of liquid metal batteries, in which convection is induced by the Joule heating of electrolyte. It is demonstrated that three-dimensional convection cells exist at subcritical Rayleigh numbers.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is whether there is subcritical convection in a horizontally - oriented fluid layer that is heated uniformly internally. Specifically, the focus of the research is on the convective behavior in the electrolyte layer caused by Joule heating in liquid metal batteries (LMBs).
### Problem Background
Liquid metal batteries are devices for short - term large - scale electrical energy storage. Their structure usually consists of three layers of liquid: a heavy metal layer (such as bismuth, antimony or lead - tin alloy) at the bottom, a light metal layer (such as sodium, lithium or magnesium) at the top, and a molten - salt electrolyte sandwiched in between. During the charging and discharging processes, a strong vertical current passes through these layers, resulting in Joule heating in the electrolyte layer and thus triggering thermal convection.
### Research Motivations
1. **Convection Caused by Joule Heating**: Recent numerical simulations have shown that Joule heating in the electrolyte layer can lead to thermal convection. This convection has been observed in laboratory prototypes and is expected to become more intense and turbulent in large - scale commercial batteries.
2. **Magnetohydrodynamic Effects**: Although the current in the battery generates a magnetic field, research has shown that, unless the battery size is very large, the influence of this magnetohydrodynamic effect on thermal convection can be negligible.
3. **Existence of Subcritical Convection**: For certain configurations, it is known that subcritical convection may exist in the pure conduction state within the range of Rayleigh number \( Ra_E < Ra < Ra_L \), but at present, there is a lack of specific description and confirmation of this convection.
### Research Objectives
The main objective of this paper is to explore, through numerical simulation, whether subcritical convection exists under given conditions (i.e., a single - layer electrolyte model with free - stress, constant - temperature boundaries) and to study the characteristics of these convections. Specifically:
- Verify the existence and stability of subcritical convection.
- Describe the flow - field structure and heat - transfer characteristics of subcritical convection.
- Explore the influence of different parameters (such as Rayleigh number, Prandtl number, etc.) on subcritical convection.
### Main Findings
- **Existence of Subcritical Convection**: The research has found that, in the case where the Rayleigh number is lower than the linear instability threshold \( Ra_L \), there do exist stable three - dimensional convection cells.
- **Convection Characteristics**: These convection cells are weak, resulting in only a slight increase in the Nusselt number, and the temperature distribution is close to the base state.
- **Spatial Structure**: The convection cells run through the entire height of the fluid layer, including the upper part with unstable stratification and the lower part with stable stratification. The downwelling regions are approximately circular in the horizontal plane, and the maximum downwelling velocity is significantly greater than the upwelling velocity.
### Conclusions
This research has successfully proven the existence of subcritical convection and provided a detailed description of its characteristics. This provides an important basis for understanding the thermal convection behavior in liquid metal batteries and also lays the foundation for further research on more complex convection phenomena.