No Arabic abstract
The Tayler instability is a kink-type flow instability which occurs when the electrical current through a conducting fluid exceeds a certain critical value. Originally studied in the astrophysical context, the instability was recently shown to be also a limiting factor for the upward scalability of liquid metal batteries. In this paper, we continue our efforts to simulate this instability for liquid metals within the framework of an integro-differential equation approach. The original solver is enhanced by multi-domain support with Dirichlet-Neumann partitioning for the static boundaries. Particular focus is laid on the detailed influence of the axial electrical boundary conditions on the characteristic features of the Tayler instability, and, secondly, on the occurrence of electro-vortex flows and their relevance for liquid metal batteries.
Lithium metal cells are key towards achieving high specific energy and energy density for electrification of transportation and aviation. Anode-free cells are the limiting case of lithium metal cells involving no excess lithium and the highest possible specific energy. In addition, anode-free cells are easier, cheaper and safer as they avoid handling and manufacturing of lithium metal foils. Issues related to dendrite growth and poor cycling are magnified in anode-free cells due to lack of excess lithium. Electrolyte and current collector surface play a crucial role in affecting the cycling performance of anode-free cells. In this work, we have computationally screened for candidate current collectors that can nucleate lithium effectively and allow uniform growth. These are determined by the free energy of lithium adsorption and lithium surface diffusion barrier on candidate current collectors. Using density functional theory calculations, we show that Li-alloys possess ideal characteristics for Li nucleation and growth. These can lead to vastly improved specific energy compared to current transition metal current collectors.
We present results for the equilibrium statistics and dynamic evolution of moderately large ($n = {mathcal{O}}(10^2 - 10^3)$) numbers of interacting point vortices on the unit sphere under the constraint of zero mean angular momentum. We consider a binary gas consisting of equal numbers of vortices with positive and negative circulations. When the circulations are chosen to be proportional to $1/sqrt{n}$, the energy probability distribution function, $p(E)$, converges rapidly with $n$ to a function that has a single maximum, corresponding to a maximum in entropy. Ensemble-averaged wavenumber spectra of the nonsingular velocity field induced by the vortices exhibit the expected $k^{-1}$ behavior at small scales for all energies. The spectra at the largest scales vary continuously with the inverse temperature $beta$ of the system and show a transition from positively sloped to negatively sloped as $beta$ becomes negative. The dynamics are ergodic and, regardless of the initial configuration of the vortices, statistical measures simply relax towards microcanonical ensemble measures at all observed energies. As such, the direction of any cascade process measured by the evolution of the kinetic energy spectrum depends only on the relative differences between the initial spectrum and the ensemble mean spectrum at that energy; not on the energy, or temperature, of the system.
In this work, the electrohydrodynamic (EHD) instability induced by a unipolar charge injection is extended from a single-phase dielectric liquid to a two-phase system that consists of a liquid-air interface. A volume of fluid (VOF) model based two-phase solver was developed with simplified Maxwell equations implemented in the open-source platform OpenFOAMtextsuperscript. The numerically obtained critical value for the linear stability matches well with the theoretical values. To highlight the effect of the slip boundary at interface, the deformation of the interface is ignored. A bifurcation diagram with hysteresis loop linking the linear and finite amplitude criteria, which is Uf = 0.059, was obtained in this situation. It is concluded that the lack of viscous effect at interface leads to a significant increase in the flow intensity, which is the reason for the smaller instability threshold in two-phase system. The presence of interface also changes the flow structure and makes the flow vortices shift closer to the interface.
Velocity gradient is the basis of many vortex recognition methods, such as Q criterion, $Delta$ criterion, $lambda_{2}$ criterion, $lambda_{ci}$ criterion and $Omega$ criterion, etc.. Except the $lambda_{ci}$ criterion, all these criterions recognize vortices by designing various invariants, based on the Helmholtz decomposition that decomposes velocity gradient into strain rate and spin. In recent years, the intuition of no vortex in straight flows has promoted people to analyze the vortex state directly from the velocity gradient, in which vortex can be distinguished from the situation that the velocity gradient has couple complex eigenvalues. A specious viewpoint to adopt the simple shear as an independent flow mode was emphasized by many authors, among them, Kolar proposed the triple decomposition of motion by extracting a so-called effective pure shearing motion; Li et al. introduced the so-called quaternion decomposition of velocity gradient and proposed the concept of eigen rotation; Liu et al. further mined the characteristic information of velocity gradient and put forward an effective algorithm of Liutex, and then developed the vortex recognition method. However, there is another explanation for the increasingly clear representation of velocity gradient, that is the local streamline pattern based on critical-point theory. In this paper, the tensorial expressions of the right/left real Schur forms of velocity gradient are clarified from the characteristic problem of velocity gradient. The relations between the involved parameters are derived and numerically verified. Comparing with the geometrical features of local streamline pattern, we confirm that the parameters in the right eigen-representation based on the right real Schur form of velocity gradient have good meanings to reveal the local streamline pattern. Some illustrative examples from the DNS data are presented.
In a cylindrical container filled with an eutectic GaInSn alloy, an electro-vortex flow (EVF) is generated by the interaction of a non-uniform current with its own magnetic field. In this paper, we investigate the EVF phenomenon numerically and experimentally. Ultrasound Doppler Velocimetry (UDV) is applied to measure the velocity field in a cylindrical vessel. Second, we enhance an old numerical solver by taking into account the effect of Joule heating, and employ it for the numerical simulation of the EVF experiment. Special focus is laid on the role of the magnetic field, which is the combination of the current induced magnetic field and the external geomagnetic field. For getting a higher computational efficiency, the so-called parent-child mesh technique is applied in OpenFOAM when computing the electric potential, the current density and the temperature in the coupled solid-liquid conductor system. The results of the experiment are in good agreement with those of the simulation. This study may help to identify the factors that are essential for the EVF phenomenon, and for quantifying its role in liquid metal batteries.