HelioCoreNode 类是一个用于模拟非线性动力学系统（如临界失稳边界附近的演化）的核心组件。其实现基于一组耦合的常微分方程，并通过欧拉法进行数值积分。

核心实现解析

1. 类定义与初始化

class HelioCoreNode:
    def __init__(self, rho: float, alpha: float, rho_c: float, alpha_c: float, eta: float):
        # G层本体状态
        self.rho = rho
        self.alpha = alpha # 临界失稳边界 self.rho_c = rho_c self.alpha_c = alpha_c
        # B层归藏阻尼（裂隙吞吐阻力）
        self.eta = eta

• 状态变量：rho（实存深度）和 alpha（张力强度）是系统的两个核心状态。
• 临界参数：rho_c 和 alpha_c 定义了系统的临界点或失稳边界。
• 阻尼系数：eta 控制系统的耗散或回归速率。

2. 动力学方程与数值积分 (step 方法)

def step(self, dt: float):
    """LD-002 非线性动力学步进（欧拉数值积分）
    dρ/dt = -η*(ρ - ρ_c) - alpha*ρ  本体回归临界区间 + 张力消耗实存
    dα/dt = (rho_c - rho) - eta*alpha  本体偏差生成张力 + 归藏阻尼缓释张力
    """
    drho_dt = -self.eta * (self.rho - self.rho_c) - self.alpha * self.rho
    dalpha_dt = (self.rho_c - self.rho) - self.eta * self.alpha self.rho += drho_dt * dt self.alpha += dalpha_dt * dt

• drho_dt：由两项组成。-eta*(rho - rho_c) 表示系统向临界点 rho_c 的线性回归趋势；-alpha*rho 表示张力 alpha 对实存深度 rho 的非线性消耗。
• dalpha_dt：同样由两项组成。(rho_c - rho) 表示当 rho 偏离临界点时会产生张力；-eta*alpha 表示张力的线性耗散。
• 积分方法：采用显式欧拉法，使用当前时间步的状态计算导数，并乘以时间步长 dt 来更新状态。

3. 扰动与观测方法

def perturb(self, delta_rho: float, delta_alpha: float):
    """瞬时外部扰动，放大G-S裂隙"""
    self.rho += delta_rho
    self.alpha += delta_alpha

def observe(self, t: float) -> Dict:
    """快照观测，输出三层核心指标"""
    # S层衍生指标：符号演化率γ
    gamma = np.abs(self.alpha / (self.rho + 1e-8))
    # B层归藏指标：距临界边界拓扑距离
    dist_critical = np.abs(self.rho - self.rho_c) + np.abs(self.alpha - self.alpha_c)
    return {
        "t": t,
        "rho": self.rho,
        "alpha": self.alpha,
        "gamma": gamma,
        "dist_to_critical": dist_critical
    }

• perturb：对状态变量 rho 和 alpha 施加瞬时、加性的扰动。
• observe：返回包含原始状态和衍生指标的快照。
  • gamma：定义为 |alpha / rho|，可解释为张力与实存的比率，其临界值通常为1。
  • dist_to_critical：使用曼哈顿距离度量当前状态与临界点的距离。

典型使用模式HelioCoreNode 通常被封装在 Experiment 类中进行系统性模拟和可视化。标准工作流程如下：

# 1. 配置并创建实验
exp = Experiment(
    name="稳定性测试",
    rho_0=0.9, # 初始实存深度 alpha_0=0.48,     # 初始张力强度 delta_rho=0.05,   # 对rho的初始扰动
    delta_alpha=0.02, # 对alpha的初始扰动
    T=50.0,           # 模拟总时长
    dt=0.02,          # 积分步长    rho_c=1.0,        # 临界实存深度 alpha_c=0.5,      # 临界张力强度 eta=0.3           # 阻尼系数
)

# 2. 运行模拟
times, history = exp.run()
# `history` 是包含每个时间步 `observe()` 输出字典的列表

# 3. 可视化结果
fig = exp.plot(times, history)
plt.show()

关键参数影响分析

应用场景与扩展

1. 稳定性分析：通过设置不同的初始条件 (rho_0, alpha_0) 和扰动，观察系统是回归原状态、收敛到新平衡点还是失稳发散。
2. 参数扫描研究：系统性地改变 eta、rho_c 等参数，研究它们对系统稳定边界和演化轨迹的影响。
3. 模型扩展：
   • 更复杂的动力学：可在 step 方法中修改微分方程，例如引入非线性项、随机噪声或时变参数。
   • 高级积分器：将欧拉法替换为龙格-库塔法等精度更高的数值积分方法。
   • 网络耦合：创建多个 HelioCoreNode 实例，并通过耦合项连接它们的 step 方法，以模拟相互作用的节点网络。