Numerical solution of the initial boundary value problem with vacuum boundary conditions for the magnetic field induction equation in a ball
The aim of this work is to develop and test an algorithm for numerical solution of the initialboundary value problem with vacuum boundary conditions for the equation of magnetic field induction in a ball based on a modification of finite-difference methods of computational fluid dynamics and electrodynamics for use in orthogonal curvilinear coordinates. The problem of the numerical solution of the magnetic field induction equation in the ball G is considered in this paper using the magnetohydrodynamics (MHD) model. It is assumed that outside the ball the magnetic field induction B = grad ψ, and the potential ψ is a harmonic function regular at infinity. At the boundary ∂G , the vacuum boundary conditions are set. They include the requirement that the magnetic field is continuous and the normal component of the current density is equal to zero. The discretization of the induction equation is carried out using the control-volume method and the FDTD (Finite Difference Time Domain) method modified for orthogonal curvilinear coordinates. When integrating over time, a completely implicit scheme is used. The total (convective and diffusion) flows are approximated according to the power-law scheme. The proposed algorithm consists in sequentially executing the following steps at each time layer: Step 1: The approximation of the radial component of the magnetic field at ∂G is determined from the discretized continuity equation for the boundary control volumes. Step 2: The external Neumann problem for the Laplace equation is solved using the Kelvin transformation. The approximation for the tangential components of the magnetic field at the ∂G is calculated using the found potential ψ. Step 3: The found components of B are used as boundary conditions for finding a solution to the induction equation. The steps 1 to 3 are repeated until convergence is achieved at a given moment in time. The algorithm was tested on the problem of the magnetic field diffusion in a conducting ball, which has an analytical solution. In this paper, the discretization of the magnetic field induction equation the solution of which satisfies the solenoidality condition is obtained in arbitrary orthogonal curvilinear coordinates. A new algorithm for numerical solution of the initial-boundary value problem with vacuum boundary conditions for the equation of magnetic field induction in a ball is developed. The comparison between the numerical solution and the analytical solution demonstrates the convergence of the constructed finite-difference scheme.
Keywords
уравнение индукции,
вакуумные граничные условия,
численное решение,
induction equation,
vacuum boundary conditions,
numerical solutionAuthors
| Bychin Igor V. | Surgut State University; Scientific Research Institute for System Analysis of the Russian Academy of Sciences | igor-bychin@yandex.ru |
| Gorelikov Andrey V. | Surgut State University; Scientific Research Institute for System Analysis of the Russian Academy of Sciences | gorelikov_a@list.ru |
| Ryakhovskii Alexei V. | Surgut State University | echo47@rambler.ru |
Всего: 3
References
Olson P.L., Glatzmaier G.A., Coe R.S. Complex polarity reversals in a geodynamo model // Earth Planet. Sci. Lett. 2011. V. 304. P. 168-179. DOI: 10.1016/j.epsl.2011.01.031.
Olson P., Dr.scoll P., Am.t H. Dipole collapse and reversal precursors in a numerical dynamo // Phys. Earth. Planet. Inter. 2009. V. 173. P. 121-140. DOI: 10.1016/j.pepi.2008.11.010.
Jouve L., Brun A.S., Arlt R., Brandenburg A., Dikpati M., Bonanno A., Kapyla P.J., Moss D., Rempel M., G.lman P., Korp. M.J., Kosov.chev A.G. A solar mean field dynamo benchmark // Astron. Astrophys. 2008. V. 483. P. 949-960. DOI: 10.1051/0004-6361:20078351.
Glatzmaier G.A. Computer simulations of Jupiter's deep internal dynamics help interpret what Juno sees // Proc. Nat. Acad. Sci. 2018. V. 115 Iss. 27. P. 6896-6904. DOI: 10.1073/pnas. 1709125115.
Jones C. A dynamo model of Jupiter's magnetic field // Icarus. 2014. V. 241. P. 148-159. DOI: 10.1016/j.icarus.2014.06.020.
Jones C.A. Convection-driven geodynamo models // Phil. Trans. R. Soc. Lond. A. 2000. V. 358. P. 873-897. DOI: 10.1098/rsta.2000.0565.
Jackson A.A. et al. A spherical shell numerical dynamo benchmark with pseudo-vacuum magnetic boundary conditions // Geophysical Journal International. 2014. V. 196 Iss. 2. P. 712-723. DOI: 10.1093/gji/ggt425.
Bandaru V., Boeck T., Krasnov D., Schumacher J. A hybrid finite difference/boundary element procedure for the simulation of turbulent MHD duct flow at finite magnetic Reynolds number // Journal of Computational Physics. 2016. V. 304. P. 320-339. DOI: 10.1016/j.jcp. 2015.10.007.
Glatzmaier G.A. Numerical simulations of stellar convective dynamos I. The model and method // Journal of Computational Physics. 1984. V. 55. P. 461-484. DOI: 10.1016/0021-9991(84)90033-0.
Reshetnyak M.Yu. Simulation in Geodynamo. Saarbrucken: Lambert, 2013. 180 p.
Kageyama A., Miyagoshi T., Sato T. Formation of current coils in geodynamo simulations //Nature. 2008. V. 454. P. 1106-1109. DOI: 10.1038/nature07227.
Галанин М.П., Попов Ю.П. Квазистационарные электромагнитные поля в неоднородных средах: Математическое моделирование. М.: Наука; Физматлит, 1995. 320 с.
Hejda P., Reshetnyak M. Control volume method for the dynamo problem in the sphere with the free rotating inner core // Studia geoph. et geod. 2003. V. 47. P. 147-159. DOI: 10.1023/A:1022207823737.
Куликовский А.Г., Погорелов Н.В, Семенов А.Ю. Математические вопросы численного решения гиперболических систем уравнений. М: Физматлит, 2001. 608 с.
Бетелин В.Б., Галкин В.А., Гореликов А.В. Алгоритм типа предиктор - корректор для численного решения уравнения индукции в задачах магнитной гидродинамики вязкой несжимаемой жидкости // Доклады Академии наук. 2015. T. 464. № 5. C. 525-528.
Колмычков В.В., Мажорова О.С., Федосеев Е.Э. Численный метод решения уравнений МГД // Препринты ИПМ им. М.В.Келдыша. 2009. № 30. 28 с.
Toth G. The ∇ ⋅ B =0 constraint in shock-capturing magnetohydrodynamics codes // Journal of Computational Physics. 2000. V. 161. P. 605-652. DOI: 10.1006/jcph.2000.6519.
Попов М.В., Елизарова Т.Г. Моделирование трехмерных МГД-течений в рамках магнитной квазигазодинамики // Препринты ИПМ им. М.В.Келдыша. 2013. № 23. 32 с.
Powell K.G., Roe P.L., Linde T.J., Gombosi T.I., De Zeeuw D.L. A Solution-Adaptive Upwind Scheme for Ideal Magnetohydrodynamics // Journal of Computational Physics. 1999. V. 154. p. 284-309. DOI: 10.1006/jcph.1999.6299.
Teyssier R., Commercon B. Numerical Methods for Simulating Star Formation // Frontiers in Astronomy and Space Sciences. 2019. V. 6. Article 51. DOI: 10.3389/fspas.2019.00051.
Ландау Л.Д., Лифшиц Е.М. Теоретическая физика: учеб. пособ. для вузов: в 10 т. Т. VIII. Электродинамика сплошных сред. М.: Физматлит, 2005. 656 с.
Куликовский А.Г., Любимов Г.А. Магнитная гидродинамика: учеб. пособие. 3-е изд. М.: Логос, 2011. 324 с.
Cowling T.G. Magnetohydrodynamics. Bristol: Adam Hilger, 1976. 136 p.
Patankar S. Numerical Heat Transfer and Fluid Flow. Washington DC: Hemisphere Publishing, 1980. 197 p.
Versteeg H.K., Malalasekera W. An Introduction to Computational Fluid Dynamics. Harlow: Pearson Education Limited, 2007. 503 p.
Ferziger J.H., Peric M. Computational Methods for Fluid Dynamics. Berlin: Springer, 2002. 423 p.
Taflove A., Hagness S. Computational Electrodynamics: The Finite-Difference Time-Domain Method, 2nd edition. Boston: Artech House, 2000. 866 p.
Григорьев А.Д. Методы вычислительной электродинамики М.: Физматлит, 2012. 432 c.
Ates A., Altun O., Kilicman A. On A Comparison of Numerical Solution Methods for General Transport Equation on Cylindrical Coordinates // Applied Mathematics & Information Sciences. 2017. V. 11. No. 2. P. 433-439. DOI: 10.18576/amis/110211.
Dewangan S.K., Sinha S.L. Comparison of various numerical differencing schemes in predicting non-newtonian transition flow through an eccentric annulus with inner cylinder in rotation // International Journal of Mechanical and Industrial Engineering. 2013. V. 3. Iss. 1. P. 119-129.
Arshad K., Rafique M., Majid A., Jabeen S. Analyses of MHD Pressure Drop in a Curved Bend for Different Liquid Metals // Journal of Applied Sciences. 2007. V. 7. Iss. 1. P. 72-78. DOI: 10.3923/jas.2007.72.78.
Тихонов А.Н., Самарский А.А. Уравнения математической физики: учебник для ун-тов. 4-е изд., испр. М.: Наука, 1972. 746 с.
Гореликов А.В. Аксиально-симметричное решение задачи о размагничивании проводящего шара // Вестник кибернетики. 2018. № 4(32). C. 9-15.
