A numerical study of a bubble pair rising side by side in external magnetic fields

2021 ◽  
Vol 926 ◽  
Author(s):  
Jie Zhang ◽  
Ming-Jiu Ni

The motion of a pair of bubbles rising side by side under the influence of external magnetic fields is numerically examined. Through solving the fully three-dimensional Navier–Stokes equations, the results reveal that the bubble interactions are rather sensitive to the field direction and strength. At first, we identify that, in a hydrodynamic flow, whether the two bubbles will bounce or coalesce depends on the developments of the counter-rotating streamwise vortices during the collision. In particular, for an originally bouncing bubble pair, a streamwise magnetic field tends to promote their coalescence by weakening the strengths of the standing streamwise vortices, and such a weakening effect is caused by the asymmetric distribution of the Lorentz force in the presence of another bubble such that a torque is induced to offset the original streamwise vortices. Under a horizontal magnetic field, on the other hand, the influences are highly dependent on the angle between the bubble centroid line and the field: a transverse field or a moderate spanwise field always leads the bubble pair to coalescence while a strong spanwise field has the opposite effect. This anisotropic effect comes from the Lorentz force induced flow diffusion along the magnetic field, which not only produces two pairs of streamwise vortices at the bubble rear, but also homogenizes the pressure along the magnetic lines. As the competition between the two mechanisms varies with the magnetic direction and strength, the interaction between the bubble pair also changes. We show that the external magnetic fields control the bubble interaction through reconstructing the vortex structures, and hence the core mechanisms are identified.

Processes ◽  
2021 ◽  
Vol 9 (5) ◽  
pp. 841
Author(s):  
Yuzhen Jin ◽  
Huang Zhou ◽  
Linhang Zhu ◽  
Zeqing Li

A three-dimensional numerical study of a single droplet splashing vertically on a liquid film is presented. The numerical method is based on the finite volume method (FVM) of Navier–Stokes equations coupled with the volume of fluid (VOF) method, and the adaptive local mesh refinement technology is adopted. It enables the liquid–gas interface to be tracked more accurately, and to be less computationally expensive. The relationship between the diameter of the free rim, the height of the crown with different numbers of collision Weber, and the thickness of the liquid film is explored. The results indicate that the crown height increases as the Weber number increases, and the diameter of the crown rim is inversely proportional to the collision Weber number. It can also be concluded that the dimensionless height of the crown decreases with the increase in the thickness of the dimensionless liquid film, which has little effect on the diameter of the crown rim during its growth.


Author(s):  
Alexander Vakhrushev ◽  
Abdellah Kharicha ◽  
Ebrahim Karimi-Sibaki ◽  
Menghuai Wu ◽  
Andreas Ludwig ◽  
...  

AbstractA numerical study is presented that deals with the flow in the mold of a continuous slab caster under the influence of a DC magnetic field (electromagnetic brakes (EMBrs)). The arrangement and geometry investigated here is based on a series of previous experimental studies carried out at the mini-LIMMCAST facility at the Helmholtz-Zentrum Dresden-Rossendorf (HZDR). The magnetic field models a ruler-type EMBr and is installed in the region of the ports of the submerged entry nozzle (SEN). The current article considers magnet field strengths up to 441 mT, corresponding to a Hartmann number of about 600, and takes the electrical conductivity of the solidified shell into account. The numerical model of the turbulent flow under the applied magnetic field is implemented using the open-source CFD package OpenFOAM®. Our numerical results reveal that a growing magnitude of the applied magnetic field may cause a reversal of the flow direction at the meniscus surface, which is related the formation of a “multiroll” flow pattern in the mold. This phenomenon can be explained as a classical magnetohydrodynamics (MHD) effect: (1) the closure of the induced electric current results not primarily in a braking Lorentz force inside the jet but in an acceleration in regions of previously weak velocities, which initiates the formation of an opposite vortex (OV) close to the mean jet; (2) this vortex develops in size at the expense of the main vortex until it reaches the meniscus surface, where it becomes clearly visible. We also show that an acceleration of the meniscus flow must be expected when the applied magnetic field is smaller than a critical value. This acceleration is due to the transfer of kinetic energy from smaller turbulent structures into the mean flow. A further increase in the EMBr intensity leads to the expected damping of the mean flow and, consequently, to a reduction in the size of the upper roll. These investigations show that the Lorentz force cannot be reduced to a simple damping effect; depending on the field strength, its action is found to be topologically complex.


2016 ◽  
Vol 791 ◽  
pp. 568-588 ◽  
Author(s):  
Andrew D. Gilbert ◽  
Joanne Mason ◽  
Steven M. Tobias

In the process of flux expulsion, a magnetic field is expelled from a region of closed streamlines on a $TR_{m}^{1/3}$ time scale, for magnetic Reynolds number $R_{m}\gg 1$ ($T$ being the turnover time of the flow). This classic result applies in the kinematic regime where the flow field is specified independently of the magnetic field. A weak magnetic ‘core’ is left at the centre of a closed region of streamlines, and this decays exponentially on the $TR_{m}^{1/2}$ time scale. The present paper extends these results to the dynamical regime, where there is competition between the process of flux expulsion and the Lorentz force, which suppresses the differential rotation. This competition is studied using a quasi-linear model in which the flow is constrained to be axisymmetric. The magnetic Prandtl number $R_{m}/R_{e}$ is taken to be small, with $R_{m}$ large, and a range of initial field strengths $b_{0}$ is considered. Two scaling laws are proposed and confirmed numerically. For initial magnetic fields below the threshold $b_{core}=O(UR_{m}^{-1/3})$, flux expulsion operates despite the Lorentz force, cutting through field lines to result in the formation of a central core of magnetic field. Here $U$ is a velocity scale of the flow and magnetic fields are measured in Alfvén units. For larger initial fields the Lorentz force is dominant and the flow creates Alfvén waves that propagate away. The second threshold is $b_{dynam}=O(UR_{m}^{-3/4})$, below which the field follows the kinematic evolution and decays rapidly. Between these two thresholds the magnetic field is strong enough to suppress differential rotation, leaving a magnetically controlled core spinning in solid body motion, which then decays slowly on a time scale of order $TR_{m}$.


Author(s):  
Lijun Liu ◽  
Koichi Kakimoto

In order to control the impurity distribution and remove defects in a crystal grown in Czochralski growth for high quality crystals of silicon, it is necessary to study and control the melt-crystal interface shape, which plays an important role in control of the crystal quality. The melt-crystal interface interacts with and is determined by the convective thermal flow of the melt in the crucible. Application of magnetic field in the Czochralski system is an effective tool to control the convective thermal flow in the crucible. Therefore, the shape of the melt-crystal interface can be modified accordingly. Numerical study is performed in this paper to understand the effect of magnetic field on the interface deflection in Czochralski system. Comparisons have been carried out by computations for four arrangements of the magnetic field: without magnetic field, a vertical magnetic field and two types of cusp-shaped magnetic field. The velocity, pressure, thermal and electromagnetic fields are solved with adaptation of the mesh to the iteratively modified interface shape. The multi-block technique is applied to discretize the melt field in the crucible and the solid field of silicon crystal. The unknown shape of the melt-crystal interface is achieved by an iterative procedure. The computation results show that the magnetic fields have obvious effects on both the pattern and strength of the convective flow and the interface shape. Applying magnetic field in the Czochralski system, therefore, is an effective tool to control the quality of bulk crystal in Czochralski growth process.


2001 ◽  
Author(s):  
Haim H. Bau

Abstract In this paper, I review some of our work on the use of magneto hydrodynamics (MHD) for pumping, controlling, and stirring fluids in microdevices. In many applications, one operates with liquids that are at least slightly conductive such as biological fluids. By patterning electrodes inside flow conduits and subjecting these electrodes to potential differences, one can induce electric currents in the liquid. In the presence of a magnetic field, a Lorentz force is generated in a direction that is perpendicular to both the magnetic and electric fields. Since one has a great amount of freedom in patterning the electrodes, one can induce forces in various directions so as to generate complex flows including “guided” flows in virtual, wall-less channels. The magnetic flux generators can be either embedded in the device or be external. Despite their unfavorable scaling (the magnitude of the forces is proportional to the fluid volume), MHD offers many advantages such as the flexibility of applying forces in any desired direction and the ability to adjust the magnitude of the forces by adjusting either the electric and/or magnetic fields. We provide examples of (i) MHD pumps; (ii) controlled networks of conduits in which each conduit is equipped with a MHD actuator and by controlling the voltage applied to each actuator, one can direct the liquid to flow in any desired way without a need for valves; and (iii) MHD stirrers including stirrers that exhibit chaotic advection.


Author(s):  
Jean Franc¸ois Sigrist ◽  
Christian Laine ◽  
Dominique Lemoine ◽  
Bernard Peseux

This paper is related to the study of a nuclear propulsion reactor prototype for the French Navy. This prototype is built on ground and is to be dimensioned toward seismic loading. The dynamic analysis takes the coupled fluid structure analysis into account. The basic fluid models used by design engineers are inviscid incompressible or compressible. The fluid can be described in a bidimensional by slice or a three-dimensional approach. A numerical study is carried out on a generic problem for the linear FSI dynamic problem. The results of this study are presented and discussed. As a conclusion, the three-dimensional inviscid incompressible fluid appears to be the best compromise between the description of physical phenomena and the cost of modeling. The geometry of the reactor is such that large displacements of the structure in the fluid can occur. Therefore, the linearity hypothesis might not be longer valid. The case of large amplitude imposed oscillating motion of a cylinder in a confined fluid is numerically studied. A CFD code is used to investigate the fluid behavior solving the NAVIER-STOKES equations. The forces induced on the cylinder by the fluid are computed and compared to the linear solution. The limit of the linear model can then be exhibited.


2018 ◽  
Author(s):  
Jiajun Chen ◽  
Yue Sun ◽  
Hang Zhang ◽  
Dakui Feng ◽  
Zhiguo Zhang

Mixing in pipe junctions can play an important role in exciting force and distribution of flow in pipe network. This paper investigated the cross pipe junction and proposed an improved plan, Y-shaped pipe junction. The numerical study of a three-dimensional pipe junction was performed for calculation and improved understanding of flow feature in pipe. The filtered Navier–Stokes equations were used to perform the large-eddy simulation of the unsteady incompressible flow in pipe. From the analysis of these results, it clearly appears that the vortex strength and velocity non-uniformity of centerline, can be reduced by Y-shaped junction. The Y-shaped junction not only has better flow characteristic, but also reduces head loss and exciting force. The results of the three-dimensional improvement analysis of junction can be used in the design of pipe network for industry.


2021 ◽  
Author(s):  
Leily Abidi

A three dimensional numerical simulation of the effect of an axial magnetic field on the fluid flow, heat and mass transfer within the solvent of GE0.98Si0.02 grown by the travelling solvent method is presented. The full steady state Navier-Stokes equations, as well as the energy, continuity and the mass transport equations, were solved numerically using the finite element technique. It is found that a strong convective flow exists in the solvent, which is known to be undesirable to achieve a uniform crystal. An external axial magnetic field is applied to suppress this convection. By increasing the magnetic induction, it is observed that the intensity of the flow at the centre of the crucible reduces at a faster rate than near the wall. This phenomenon creates a stable and flat growth interface and the silicon distribution in the horizontal plane becomes relatively homocentric. The maximum velocity is found to obey a power law with respect to the Hartmann number Umax Ha⁻⁷/⁴


Author(s):  
N. Lymberopoulos ◽  
K. Giannakoglou ◽  
I. Nikolaou ◽  
K. D. Papailiou ◽  
A. Tourlidakis ◽  
...  

Mechanical constraints dictate the existence of tip clearances in rotating cascades, resulting to a flow leakage through this clearance which considerably influences the efficiency and range of operation of the machine. Three-dimensional Navier-Stokes solvers are often used for the numerical study of compressor and turbine stages with tip-clearance. The quality of numerical predictions depends strongly on how accurately the blade tip region is modelled; in this respect the accurate modelling of tip region was one of the main goals of this work. In the present paper, a 3-D Navier-Stokes solver is suitably adapted so that the flat tip surface of a blade and its sharp edges could be accurately modelled, in order to improve the precision of the calculation in the tip region. The adapted code solves the fully elliptic, steady, Navier-Stokes equations through a space-marching algorithm and a pressure correction technique; the H-type topology is retained, even in cases with thick leading edges where a special treatment is introduced herein. The analysis is applied to two different cases, a linear cascade and a compressor rotor, and comparisons with experimental data are provided.


Sign in / Sign up

Export Citation Format

Share Document