[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: What is the difference between a host and an end system? List several different types of end...
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
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