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 »
- Publication Date:
- Grant/Contract Number:
- AC02-09CH11466
- Type:
- Accepted Manuscript
- Journal Name:
- IEEE Transactions on Plasma Science
- Additional Journal Information:
- Journal Volume: 45; Journal Issue: 9; Journal ID: ISSN 0093-3813
- Publisher:
- IEEE
- Research Org:
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Sponsoring Org:
- USDOE
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; computational fluid dynamics (CFD); fusion; magnetohydrodynamics (MHD)
- OSTI Identifier:
- 1416157
Khodak, Andrei. Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications. United States: N. p.,
Web. doi:10.1109/TPS.2017.2734106.
Khodak, Andrei. Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications. United States. doi:10.1109/TPS.2017.2734106.
Khodak, Andrei. 2017.
"Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications". United States.
doi:10.1109/TPS.2017.2734106. https://www.osti.gov/servlets/purl/1416157.
@article{osti_1416157,
title = {Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications},
author = {Khodak, Andrei},
abstractNote = {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 a 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.},
doi = {10.1109/TPS.2017.2734106},
journal = {IEEE Transactions on Plasma Science},
number = 9,
volume = 45,
place = {United States},
year = {2017},
month = {8}
}