Virtual Lab: Plasma Fusion Reactors

Authenticated as Guest

UUID: mit-fusion

EXPLORER

📄 main.py

📄 model_arch.py

📁 datasets

WebGL 3D Simulation

Initializing WebGL Engine...

Interactive Editor

# Magnetic confinement simulation
import numpy as np
import matplotlib.pyplot as plt
import io
import base64

# Simulate plasma temperature and density profiles
r = np.linspace(0, 1, 100) # normalized minor radius
# Parabolic profiles typical of L-mode confinement
temperature = 5.0 * (1 - r**2)**2 # keV
density = 1.0 * (1 - r**2) # x10^20 m^-3

plt.style.use("dark_background")
fig, ax1 = plt.subplots(figsize=(6, 4))

color = "#ef4444" # Red for temperature
ax1.set_xlabel("Normalized Minor Radius (r/a)")
ax1.set_ylabel("Temperature (keV)", color=color)
ax1.plot(r, temperature, color=color, linewidth=2)
ax1.tick_params(axis="y", labelcolor=color)

ax2 = ax1.twinx()  
color = "#3b82f6" # Blue for density
ax2.set_ylabel("Density (10^20 m^-3)", color=color)  
ax2.plot(r, density, color=color, linewidth=2, linestyle="--")
ax2.tick_params(axis="y", labelcolor=color)

plt.title("Tokamak Plasma Radial Profiles")
fig.tight_layout()

buf = io.BytesIO()
plt.savefig(buf, format="png", bbox_inches="tight", transparent=True)
buf.seek(0)
render_image(base64.b64encode(buf.read()).decode("utf-8"))

print("Tokamak Magnetic Confinement Simulation Complete.")

root@phd-lab-vbox:~# python main.py

Waiting for execution...

root@phd-lab-vbox:~#