TY - GEN
T1 - Parallel fully-implicit computation of magnetohydrodynamics acceleration experiments
AU - Wan, Tian
AU - Candler, Graham
PY - 2010
Y1 - 2010
N2 - A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.
AB - A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.
KW - Electron-temperature
KW - Joule heating
KW - fully-implicit method
UR - http://www.scopus.com/inward/record.url?scp=77955750872&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=77955750872&partnerID=8YFLogxK
U2 - 10.1063/1.3452140
DO - 10.1063/1.3452140
M3 - Conference contribution
AN - SCOPUS:77955750872
SN - 9780735407787
T3 - AIP Conference Proceedings
SP - 1553
EP - 1558
BT - ISCM II and EPMESC XII - Proc. of the 2nd Int. Symposium on Computational Mechanics and the 12th Int. Conf. on the Enhancement and Promotion of Computational Methods in Engineering and Science
T2 - 2nd International Symposium on Computational Mechanics, ISCM II, and the 12th International Conference on the Enhancement and Promotion of Computational Methods in Engineering and Science, EPMESC XII
Y2 - 30 November 2009 through 3 December 2009
ER -