skip to main content


This content will become publicly available on August 21, 2018

Title: Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications

Here, the analysis of many fusion applications such as liquid-metal blankets requires application of computational fluid dynamics (CFD) methods for electrically conductive liquids in geometrically complex regions and in the presence of a strong magnetic field. A current state of the art general purpose CFD code allows modeling of the flow in complex geometric regions, with simultaneous conjugated heat transfer analysis in liquid and surrounding solid parts. Together with a magnetohydrodynamics (MHD) capability, the general purpose CFD code will be a valuable tool for the design and optimization of fusion devices. This paper describes an introduction of MHD capability into the general purpose CFD code CFX, part of the ANSYS Workbench. The code was adapted for MHD problems using a magnetic induction approach. CFX allows introduction of user-defined variables using transport or Poisson equations. For MHD adaptation of the code three additional transport equations were introduced for the components of the magnetic field, in addition to the Poisson equation for electric potential. The Lorentz force is included in the momentum transport equation as a source term. Fusion applications usually involve very strong magnetic fields, with values of the Hartmann number of up to tens of thousands. In this situation amore » system of MHD equations become very rigid with very large source terms and very strong variable gradients. To increase system robustness, special measures were introduced during the iterative convergence process, such as linearization using source coefficient for momentum equations. The MHD implementation in general purpose CFD code was tested against benchmarks, specifically selected for liquid-metal blanket applications. Results of numerical simulations using present implementation closely match analytical solutions for a Hartmann number of up to 1500 for a 2-D laminar flow in the duct of square cross section, with conducting and nonconducting walls. Results for a 3-D test case are also included.« less
ORCiD logo [1]
  1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Publication Date:
Grant/Contract Number:
Accepted Manuscript
Journal Name:
IEEE Transactions on Plasma Science
Additional Journal Information:
Journal Volume: 45; Journal Issue: 9; Journal ID: ISSN 0093-3813
Research Org:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Sponsoring Org:
Country of Publication:
United States
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; computational fluid dynamics (CFD); fusion; magnetohydrodynamics (MHD)
OSTI Identifier: