[1] import astropy.units as u import matplotlib.pyplot as plt import numpy as np from astropy.constants import h, m_e, k_B from astropy.constants import u as atomicu def Saha (T,Z_i,Z_ii,E_ionize,P_e): N_iitoi = N_iitoi = N_iitoi.decompose() return (N_iitoi) 2*k_B*T*Z_ii/(P_e*Z_i)*(2*np.pi*m_e*k_B*T/h**2)**(3/2)*np.exp(-E_ionize/(k_B*T)) [3] def Boltzmann_H(T, low, high): g_low = 2*low**2 g_high = 2*high**2 E_low = -13.6/low**2*u.eV E_high = -13.6/high**2*u.ev exponent = −1* (E_high-E_low)/(k_B*T) exponent = exponent.decompose() N_htol = g_high/g_low*np.exp(exponent) N_htol = N_htol.decompose() return (N_htol)

Computer Networking: A Top-Down Approach (7th Edition)
7th Edition
ISBN:9780133594140
Author:James Kurose, Keith Ross
Publisher:James Kurose, Keith Ross
Chapter1: Computer Networks And The Internet
Section: Chapter Questions
Problem R1RQ
icon
Related questions
Question

Answer in Python. Given this code to start with (definitions for the Saha and Boltzmann Equations), I must: use the Saha equation (defined in the second code block) to find the temperature of the partial ionization zone for a certain gas.

(Hint: The Saha function will return answers for the Saha equation for gases of arbitrary composition) 

[1] import astropy.units as u
import matplotlib.pyplot as plt
import numpy as np
from astropy.constants import h, m_e, k_B
from astropy.constants import u as atomicu
def Saha (T,Z_i,Z_ii,E_ionize,P_e):
N_iitoi = 2*k_B*T*Z_ii/(P_e*Z_i)*(2*np.pi*m_e*k_B*T/h**2)**(3/2)*np.exp(-E_ionize/(k_B*T))
N_iitoi = N_iitoi.decompose()
return (N_iitoi)
[3] def Boltzmann_H(T,low, high):
g_low
2*low**2
g_high 2*high**2
E low =
=
-13.6/low**2*u.eV
E_high = -13.6/high**2*u.eV
exponent = -1* (E_high-E_low)/(k_B*T)
exponent = exponent.decompose()
N htol g_high/g_low*np.exp(exponent)
N_htol N_htol.decompose()
return (N_htol)
Transcribed Image Text:[1] import astropy.units as u import matplotlib.pyplot as plt import numpy as np from astropy.constants import h, m_e, k_B from astropy.constants import u as atomicu def Saha (T,Z_i,Z_ii,E_ionize,P_e): N_iitoi = 2*k_B*T*Z_ii/(P_e*Z_i)*(2*np.pi*m_e*k_B*T/h**2)**(3/2)*np.exp(-E_ionize/(k_B*T)) N_iitoi = N_iitoi.decompose() return (N_iitoi) [3] def Boltzmann_H(T,low, high): g_low 2*low**2 g_high 2*high**2 E low = = -13.6/low**2*u.eV E_high = -13.6/high**2*u.eV exponent = -1* (E_high-E_low)/(k_B*T) exponent = exponent.decompose() N htol g_high/g_low*np.exp(exponent) N_htol N_htol.decompose() return (N_htol)
Expert Solution
trending now

Trending now

This is a popular solution!

steps

Step by step

Solved in 2 steps

Blurred answer
Similar questions
  • SEE MORE QUESTIONS
Recommended textbooks for you
Computer Networking: A Top-Down Approach (7th Edi…
Computer Networking: A Top-Down Approach (7th Edi…
Computer Engineering
ISBN:
9780133594140
Author:
James Kurose, Keith Ross
Publisher:
PEARSON
Computer Organization and Design MIPS Edition, Fi…
Computer Organization and Design MIPS Edition, Fi…
Computer Engineering
ISBN:
9780124077263
Author:
David A. Patterson, John L. Hennessy
Publisher:
Elsevier Science
Network+ Guide to Networks (MindTap Course List)
Network+ Guide to Networks (MindTap Course List)
Computer Engineering
ISBN:
9781337569330
Author:
Jill West, Tamara Dean, Jean Andrews
Publisher:
Cengage Learning
Concepts of Database Management
Concepts of Database Management
Computer Engineering
ISBN:
9781337093422
Author:
Joy L. Starks, Philip J. Pratt, Mary Z. Last
Publisher:
Cengage Learning
Prelude to Programming
Prelude to Programming
Computer Engineering
ISBN:
9780133750423
Author:
VENIT, Stewart
Publisher:
Pearson Education
Sc Business Data Communications and Networking, T…
Sc Business Data Communications and Networking, T…
Computer Engineering
ISBN:
9781119368830
Author:
FITZGERALD
Publisher:
WILEY