Particle Transport and Deposition in a Turbulent Square Duct Flow With an Imposed Magnetic Field

2014 ◽  
Vol 136 (12) ◽  
Author(s):  
Rui Liu ◽  
Surya P. Vanka ◽  
Brian G. Thomas

In this paper, we study particle transport and deposition in a turbulent square duct flow with an imposed magnetic field using direct numerical simulations (DNS) of the continuous flow and Lagrangian tracking of particles. The magnetic field and the velocity induce a current and the interaction of this current with the magnetic field generates a Lorentz force that brakes the flow and modifies the flow structure. A second-order accurate finite volume method is used to integrate the coupled Navier–Stokes and magnetohydrodynamic (MHD) equations and the solution procedure is implemented on a graphics processing unit (GPU). Magnetically nonconducting particles of different Stokes numbers are continuously injected at random locations in the inlet cross section of the duct and their rates of deposition on the duct walls are studied with and without a magnetic field. Because of the modified instantaneous turbulent flow structures as a result of the magnetic field, the deposition rates and patterns on the walls perpendicular to the magnetic field are lower than those on the walls parallel to the magnetic field.

Author(s):  
Rui Liu ◽  
Surya P. Vanka ◽  
Brian G. Thomas

In this paper we study the particle transport and deposition in a turbulent square duct flow with an imposed magnetic field using Direct Numerical Simulations (DNS) of the continuous flow. A magnetic field induces a current and the interaction of this current with the magnetic field generates a Lorentz force which brakes the flow and modifies the flow structure. A second-order accurate finite volume method in time and space is used and implemented on a GPU. Particles are injected at the entrance to the duct continuously and their rates of deposition on the duct walls are computed for different magnetic field strengths. Because of the changes to the flow due to the magnetic field, the deposition rates are different on the top and bottom walls compared to the side walls. This is different than in a non-MHD square duct flow, where quadrant (and octant) symmetry is obtained.


2015 ◽  
Vol 137 (7) ◽  
Author(s):  
K. Jin ◽  
S. P. Vanka ◽  
B. G. Thomas

In this paper, we study the three-dimensional (3D) flow of an electrically conducting fluid in a cubic cavity with the top wall moving and subjected to an external magnetic field. The governing flow and electromagnetic field equations are integrated by a second-order space and time accurate numerical scheme, implemented on a graphics processing unit (GPU) with high parallel efficiency. Solutions for several Reynolds and Stuart numbers have been obtained on sufficiently fine grids to achieve grid independent solutions. As expected, the magnetic field significantly influences the circulation in the cavity and modifies the shape and locations of the primary and secondary eddies. The observed flow patterns are illustrated graphically as well as through selected line plots and tabulated data. With increasing magnetic field strength, the center of the primary eddy is seen to shift to the top right corner. Further, situations where the flow is unsteady in the absence of the magnetic field have become steady after a certain value of the magnetic interaction parameter.


2005 ◽  
Vol 128 (4) ◽  
pp. 413-417 ◽  
Author(s):  
Mohammed J. Al-Khawaja ◽  
Mohammed Selmi

A liquid metal forced-convection fully developed laminar flow inside a square duct, whose surfaces are electrically insulated and subjected to a constant temperature in a transverse magnetic field, is solved numerically using the spectral method. The axial momentum, induction, and nonlinear energy equations are solved by expanding the axial velocity, magnetic field, and temperature in double Chebyshev series and are collocated at Gauss points. The resulting system of equations is solved numerically by Gauss elimination for the expansion coefficients. The velocity and the magnetic field coefficients are directly solved for, while the temperature coefficients are solved for iteratively. Results show that the velocity profile is flattened in the direction of the magnetic field, but it is more round in the direction normal to it, in a similar fashion to the case of circular tube studied previously. The powerful spectral method resolves the sharp velocity gradient near the duct walls very well leading to accurate calculation of friction factor and Nusselt number. These parameters increase with the strength of the magnetic field due to the increasing flatness of the velocity profile. Comparison with the results for the circular tube shows that the effect of magnetic field on square duct flow is slightly lower from that one for circular pipe flow.


1998 ◽  
Vol 21 (2) ◽  
pp. 299-305 ◽  
Author(s):  
Fengxin Chen ◽  
Ping Wang ◽  
Chaoshun Qu

In this paper we study the system governing flows in the magnetic field within the earth. The system is similar to the magnetohydrodynamic (MHD) equations. For initial data in spaceLp, we obtained the local in time existence and uniqueness ofweak solutions of the system subject to appropriate initial and boundary conditions.


Author(s):  
Aaron F. Shinn ◽  
S. P. Vanka

A semi-implicit pressure based multigrid algorithm for solving the incompressible Navier-Stokes equations was implemented on a Graphics Processing Unit (GPU) using CUDA (Compute Unified Device Architecture). The multigrid method employed was the Full Approximation Scheme (FAS), which is used for solving nonlinear equations. This algorithm is applied to the 2D driven cavity problem and compared to the CPU version of the code (written in Fortran) to assess computational speed-up.


2020 ◽  
Vol 90 (3) ◽  
pp. 482
Author(s):  
Н.М. Горшунов ◽  
Е.П. Потанин

Equations are obtained that describe the characteristics of the azimuthal motion and radial expansion of a plasma jet under the action of a rotating transverse magnetic field of a dipole configuration in a longitudinal static magnetic field. The analysis was carried out both in the multicomponent approximation and on the basis of MHD equations taking into account the Hall effect. Based on the obtained dependences of the azimuthal and radial ion velocities on the magnetic field values, the separation characteristics of the direct-flow plasma centrifuge are estimated for the separation of a two-component binary mixture simulating spent nuclear fuel. It was shown that the concentration of the heavy uranium-plutonium component in the product flow can be increased from the initial 96 to 99.8% with a fuel component extraction of 0.87.


1989 ◽  
Vol 12 (1) ◽  
pp. 159-174
Author(s):  
B. D. Aggarwala ◽  
P. D. Ariel

In this paper, we investigate the flow of a viscous, Incompressible, electrically conducting fluid through a rectangular duct in the presence of a magnetic field, when one of the boundaries perpedicular to the magnetic field is partly conducting and partly Insulating, by a modified Boundary Integral Method.Three problems are considered (i) flow through an infinite channel, (ii) flow through a rectangular duct when the conducting part is symmetrically situated, and (iii) flow through a rectangular duct when the conducting part is arbltrarily positioned.Such problems have been studied before by asymptotic means for large values of M, the Hartmann number. Hoverer, the present modification of the Boundary Integral Method renders the problem computationally efficient and provides a reliable numerical solution for all values of M. For large M, our coputation time decreases significantly.


2008 ◽  
Vol 130 (9) ◽  
Author(s):  
K. M. Guleren ◽  
I. Afgan ◽  
A. Turan

The laminarization phenomenon for the flow under the combined effect of strong curvature and rotation is discussed based on numerical predictions of large-eddy simulation (LES). Initially, the laminarization process is presented for the fully developed flow inside a spanwise rotating straight square duct. LES predictions over a wide range of rotation numbers (Ro=0–5) show that the turbulent kinetic energy decreases monotonically apart from 0.2<Ro<0.5. Subsequently, a spanwise rotating U-duct flow is considered with Ro=±0.2. The interaction of curvature and Coriolis induced secondary flows enhances the turbulence for the negative rotating case, whereas this interaction ensues strong laminarization for the positive rotating case. Finally, the laminarization is presented in the impeller of a typical centrifugal compressor, rotating at a speed of Ω=1862rpm(Ro=0.6). The resulting LES predictions are observed to be better than those of Reynolds-averaged Navier-Stokes (RANS) in the regions where turbulence is significant. However, for the regions dominated by strong laminarization, RANS results are seen to approach those of LES and experiments.


2020 ◽  
Vol 633 ◽  
pp. A87 ◽  
Author(s):  
L. Griton ◽  
F. Pantellini

Context. As proven by measurements at Uranus and Neptune, the magnetic dipole axis and planetary spin axis can be off by a large angle exceeding 45°. The magnetosphere of such an (exo-)planet is highly variable over a one-day period and it does potentially exhibit a complex magnetic tail structure. The dynamics and shape of rotating magnetospheres do obviously depend on the planet’s characteristics but also, and very substantially, on the orientation of the planetary spin axis with respect to the impinging, generally highly supersonic, stellar wind. Aims. On its orbit around the Sun, the orientation of Uranus’ spin axis with respect to the solar wind changes from quasi-perpendicular (solstice) to quasi-parallel (equinox). In this paper, we simulate the magnetosphere of a fictitious Uranus-like planet plunged in a supersonic plasma (the stellar wind) at equinox. A simulation with zero wind velocity is also presented in order to help disentangle the effects of the rotation from the effects of the supersonic wind in the structuring of the planetary magnetic tail. Methods. The ideal magnetohydrodynamic (MHD) equations in conservative form are integrated on a structured spherical grid using the Message-Passing Interface-Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). In order to limit diffusivity at grid level, we used background and residual decomposition of the magnetic field. The magnetic field is thus made of the sum of a prescribed time-dependent background field B0(t) and a residual field B1(t) computed by the code. In our simulations, B0(t) is essentially made of a rigidly rotating potential dipole field. Results. The first simulation shows that, while plunged in a non-magnetised plasma, a magnetic dipole rotating about an axis oriented at 90° with respect to itself does naturally accelerate the plasma away from the dipole around the rotation axis. The acceleration occurs over a spatial scale of the order of the Alfvénic co-rotation scale r*. During the acceleration, the dipole lines become stretched and twisted. The observed asymptotic fluid velocities are of the order of the phase speed of the fast MHD mode. In two simulations where the surrounding non-magnetised plasma was chosen to move at supersonic speed perpendicularly to the rotation axis (a situation that is reminiscent of Uranus in the solar wind at equinox), the lines of each hemisphere are symmetrically twisted and stretched as before. However, they are also bent by the supersonic flow, thus forming a magnetic tail of interlaced field lines of opposite polarity. Similarly to the case with no wind, the interlaced field lines and the attached plasma are accelerated by the rotation and also by the transfer of kinetic energy flux from the surrounding supersonic flow. The tailwards fluid velocity increases asymptotically towards the externally imposed flow velocity, or wind. In one more simulation, a transverse magnetic field, to both the spin axis and flow direction, was added to the impinging flow so that magnetic reconnection could occur between the dipole anchored field lines and the impinging field lines. No major difference with respect to the no-magnetised flow case is observed, except that the tailwards acceleration occurs in two steps and is slightly more efficient. In order to emphasise the effect of rotation, we only address the case of a fast-rotating planet where the co-rotation scale r* is of the order of the planetary counter-flow magnetopause stand-off distance rm. For Uranus, r*≫ rm and the effects of rotation are only visible at large tailwards distances r ≫ rm.


Sign in / Sign up

Export Citation Format

Share Document