NUMERICAL MODELING OF THE TWO-TEMPERATURE PLASMA EQUATION IN THE PYTHON PROGRAMMING ENVIRONMENT
- 31 March 2026
- journal article
- Published by K. Zhubanov Aktobe Regional University in Научный журнал "Вестник Актюбинского регионального университета имени К. Жубанова"
- Vol. 83 (1) , 69-77
Abstract
This study presents numerical modeling and simulation of a two-temperature plasma system using the Python programming environment. This paper presents numerical modeling and simulation of a two-temperature plasma system using the Python programming environment. The model is based on a connected system of nonlinear equations of thermal conductivity describing the processes of energy exchange between electronic and ionic components. To obtain numerical solutions, an implicit finite-difference Crank–Nicholson scheme was implemented, providing second-order accuracy in time and space. The computational process was carried out using the NumPy, SciPy and SymPy libraries, and visualization of the results and user interaction were implemented using the Mathplotlib and Tkinter tools. The developed program allows you to interactively enter parameters, calculate temperature evolution and visualize the results in real time. Numerical experiments were carried out for the values of parameters corresponding to the natural and climatic conditions of the Republic of Karakalpakstan, where the density of the solar energy flux varies in the range of 1000-2800 W/m2. The results show the changes over time and space in electron and ion temperatures, confirming the Crank–Nicholson approach and the reliability of the developed software for modeling energy transfer in low-pressure plasma conditions.Keywords
This publication has 8 references indexed in Scilit:
- FINITE ELEMENT METHOD SCHEMES OF HIGHER ACCURACY FOR SOLVING NON-STATIONARY FOURTH-ORDER EQUATIONSJournal of Mathematics, Mechanics and Computer Science, 2023
- Bandpass approximate Crank-Nicolson implementation for anisotropic gyrotropic plasma open region simulationOptik, 2021
- Convergence of high-precision finite element method schemes for two-temperature plasma equationAIP Conference Proceedings, 2021
- Computationally efficient complex envelope approximate Crank–Nicolson scheme and its open region problem for anisotropic gyrotropic plasmaPhysics of Plasmas, 2020
- Introduction to Plasma Physics and Controlled Fusion2016
- Comparison of some methods for solving the internal wave propagation problem in a weakly stratified fluidMathematical Models and Computer Simulations, 2011
- The Theory of Difference Schemes2001
- A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction typeMathematical Proceedings of the Cambridge Philosophical Society, 1947