Add physics visualisation panel

This commit is contained in:
2025-12-02 03:09:58 +00:00
parent f210e00743
commit 29bf371834

View File

@@ -4,8 +4,119 @@ This module provides an interactive dashboard for visualising the physics
engine directly, demonstrating thermal-electrical coupling in real-time.
"""
from collections import deque
from dataclasses import dataclass, field
import streamlit as st
from py_dvt_ate.simulation.physics.engine import PhysicsEngine
# History buffer size for charts
HISTORY_SIZE = 500
@dataclass
class SimulationHistory:
"""Stores time series data for visualisation."""
time: deque[float] = field(default_factory=lambda: deque(maxlen=HISTORY_SIZE))
chamber_temp: deque[float] = field(
default_factory=lambda: deque(maxlen=HISTORY_SIZE)
)
case_temp: deque[float] = field(default_factory=lambda: deque(maxlen=HISTORY_SIZE))
junction_temp: deque[float] = field(
default_factory=lambda: deque(maxlen=HISTORY_SIZE)
)
output_voltage: deque[float] = field(
default_factory=lambda: deque(maxlen=HISTORY_SIZE)
)
power_dissipation: deque[float] = field(
default_factory=lambda: deque(maxlen=HISTORY_SIZE)
)
def init_session_state() -> None:
"""Initialise Streamlit session state."""
if "engine" not in st.session_state:
st.session_state.engine = PhysicsEngine(update_rate_hz=100.0)
if "history" not in st.session_state:
st.session_state.history = SimulationHistory()
if "running" not in st.session_state:
st.session_state.running = False
def step_simulation(steps: int = 10) -> None:
"""Advance the simulation by the given number of steps."""
engine: PhysicsEngine = st.session_state.engine
history: SimulationHistory = st.session_state.history
for _ in range(steps):
engine.step()
# Record current state in history
thermal = engine.get_thermal_state()
electrical = engine.get_electrical_state()
history.time.append(thermal.timestamp)
history.chamber_temp.append(thermal.chamber_temperature)
history.case_temp.append(thermal.case_temperature)
history.junction_temp.append(thermal.junction_temperature)
history.output_voltage.append(electrical.output_voltage)
history.power_dissipation.append(electrical.power_dissipation)
def display_thermal_chart() -> None:
"""Display temperature chart."""
history: SimulationHistory = st.session_state.history
if len(history.time) < 2:
st.info("Start the simulation to see temperature data")
return
chart_data = {
"Time (s)": list(history.time),
"Chamber": list(history.chamber_temp),
"Case": list(history.case_temp),
"Junction": list(history.junction_temp),
}
st.line_chart(
chart_data,
x="Time (s)",
y=["Chamber", "Case", "Junction"],
color=["#1f77b4", "#ff7f0e", "#d62728"],
)
def display_current_state() -> None:
"""Display current simulation state metrics."""
engine: PhysicsEngine = st.session_state.engine
thermal = engine.get_thermal_state()
electrical = engine.get_electrical_state()
col1, col2, col3, col4 = st.columns(4)
with col1:
st.metric("Chamber Temp", f"{thermal.chamber_temperature:.2f} °C")
with col2:
st.metric("Case Temp", f"{thermal.case_temperature:.2f} °C")
with col3:
st.metric("Junction Temp", f"{thermal.junction_temperature:.2f} °C")
with col4:
st.metric("Output Voltage", f"{electrical.output_voltage:.4f} V")
col5, col6, col7, col8 = st.columns(4)
with col5:
st.metric("Input Voltage", f"{electrical.input_voltage:.2f} V")
with col6:
st.metric("Load Current", f"{electrical.load_current * 1000:.1f} mA")
with col7:
st.metric("Power Diss.", f"{electrical.power_dissipation * 1000:.2f} mW")
with col8:
st.metric("Sim Time", f"{engine.simulation_time:.2f} s")
def main() -> None:
"""Main entry point for the Streamlit dashboard."""
@@ -23,6 +134,21 @@ def main() -> None:
"""
)
init_session_state()
# Current state display
st.subheader("Current State")
display_current_state()
# Temperature chart
st.subheader("Temperature History")
display_thermal_chart()
# Auto-refresh when running
if st.session_state.running:
step_simulation(steps=10)
st.rerun()
if __name__ == "__main__":
main()