3D MHD free surface fluid flow simulation based on magnetic-field induction equations

Abstract

The purpose of this paper is to present our recent efforts on 3D MHD model development and our results based on the technique derived from induced-magnetic-field equations. Two important features are utilized in our numerical method to obtain convergent solutions. First, a penalty factor is introduced in order to force the local divergence free condition of the magnetic fields. The second is that we extend the insulating wall thickness to ensure that the induced magnetic field at its boundaries is null. These simulation results for lithium film free surface flows under NSTX outboard mid-plane magnetic field configurations have shown that 3D MHD effects from a surface normal field gradient cause return currents to interact with surface normal fields and produce unfavorable MHD forces. This leads to a substantial change in flow pattern and a reduction in flow velocity, with most of the flow spilling over one side of the chute. These critical phenomena can not be revealed by 2D models. Additionally, a design which overcomes these undesired flow characteristics is obtained. #

Numerical simulation of dielectric bubbles coalescence under the effects of uniform magnetic field

“…Like other published studies in this field which have solved magnetohydrodynamic problems numerically [35,40], the CPU run-time by developed code in this research is also so long. The high run-time is caused by two major factors.…”

The Multicell Volume of Fluid (MC-VOF) Method for the Free Surface Simulation of MFD Flows. Part II: Experimental Verifications and Results

“…The induction B s is calculated by the r.m.s. value of the inductor line current density, J l , as follows: 5,6) .................. (3) in which, k B represents a reduction factor considering the reaction of the eddy currents induced in the Wood's metal and the nonmagnetic steel vessel, determined for the effective magnetic Reynolds number 6) ... (4) where, Cend denotes the mean value in axial direction of the end effect factor C, computed by the effective radius of the secondary, Reff, in a rotational two-pole magnetic field: 6,7) ........ (5) In Eq. (3), m = 3 is the number of phases, nw the number of series turns per phase of the winding, kw represents the winding factor, I the r.m.s.…”

The Multicell Volume of Fluid (MC-VOF) Method for the Free Surface Simulation of MFD Flows. Part I: Mathematical Model

“…The use for the S cells of a reduced conductivity, e.g. Fσ < σ, 12) yields smaller Lorentz forces driving unrealistic, too weak near-surface tangential flows, and as a result, a wrong free surface form can be calculated.…”

“…6(a), with a radius, re, and belonging to the vertical MC block of a MHS ( Fig. 4(a)), is given by .......................... (12) where the upwinded VOF function, Fe, is set to Fe = Fi,m, if the outgoing velocity, ue, is positive. If ue < 0, then Fe = Fi+1,m when the cell (i + 1, m) is an F, E or an S cell on an MHS ( Fig.…”