智能冰袋降温优化系统

 

一、实际应用场景描述

 

某三甲医院急诊科日均接诊外伤患者80余人次，其中约60%需要使用冰袋进行冷敷治疗。护理部门发现现有冰袋使用流程存在严重问题：1）护士凭经验估算冰袋使用时长，导致30%的患者出现局部冻伤；2）冰袋反复冷冻后降温效果衰减30-40%，但缺乏量化评估手段；3）急救转运过程中冰袋保温措施不当，到达医院时温度已回升至5-8℃，延误最佳处置时机；4）科室每月因冰袋采购和管理不当造成的浪费约1200元。急诊科主任需要一个科学的系统，能够根据冰袋材质、环境温度、包裹方式等因素，精确计算融化吸热过程和降温效果，优化冰袋使用时长和保存策略。

 

二、引入痛点

 

1. 经验依赖性强：传统使用基于"半小时更换"的经验法则，忽视个体差异和环境变量

2. 降温效果不可控：无法预测特定条件下冰袋能维持多长时间的有效低温区间（2-15℃）

3. 材质衰减难追踪：重复使用导致冰袋隔热层和相变材料性能退化，但缺乏量化指标

4. 转运保温盲区：救护车、救护直升机等转运环境中，缺乏动态的保温策略指导

5. 安全风险预警缺失：无法提前预警冻伤风险，依赖事后观察

6. 成本效益失衡：过度保守的使用策略导致冰袋消耗量增加，成本上升

 

三、核心逻辑讲解

 

相变传热基础理论

 

傅里叶导热定律：

 

q = -k \cdot A \cdot \frac{dT}{dx}

 

其中：

 

- q ：热流量 (W)

- k ：材料导热系数 (W/(m·K))

- A ：传热面积 (m²)

- \frac{dT}{dx} ：温度梯度 (K/m)

 

牛顿冷却定律：

 

q = h \cdot A \cdot (T_{surface} - T_{ambient})

 

其中：

 

- h ：对流换热系数 (W/(m²·K))

- T_{surface} ：表面温度 (K)

- T_{ambient} ：环境温度 (K)

 

相变潜热计算：

 

Q = m \cdot L_f

 

其中：

 

- Q ：相变吸热量 (J)

- m ：相变材料质量 (kg)

- L_f ：熔化潜热 (J/kg)，水为334 kJ/kg

 

冰袋融化动力学模型

 

基于能量平衡的微分方程：

 

\rho \cdot c_p \cdot V \cdot \frac{dT}{dt} = q_{in} - q_{out} - q_{phase}

 

展开为：

 

\rho \cdot c_p \cdot V \cdot \frac{dT}{dt} = h_{conv} \cdot A \cdot (T_{env} - T) + \frac{k_{ins} \cdot A}{d_{ins}} \cdot (T_{core} - T) - \delta_{melt} \cdot L_f \cdot \frac{dm}{dt}

 

其中：

 

- \rho ：冰的密度 (917 kg/m³)

- c_p ：比热容 (冰：2.09 kJ/(kg·K)，水：4.18 kJ/(kg·K))

- V ：冰袋体积 (m³)

- \delta_{melt} ：熔化指示函数（固相为1，液相为0）

- \frac{dm}{dt} ：熔化速率 (kg/s)

 

熔化速率模型（Stefan问题简化）

 

基于界面能量平衡的熔化速率：

 

\frac{dm}{dt} = \frac{k_{ice} \cdot A \cdot (T_{fusion} - T_{surface})}{L_f \cdot d_{effective}}

 

其中：

 

- k_{ice} ：冰的导热系数 (2.18 W/(m·K))

- T_{fusion} ：冰点 (273.15 K)

- d_{effective} ：有效传热厚度 (m)

 

有效降温时间计算

 

定义有效降温区间为 [T_{min}, T_{max}] ，计算冰袋温度降至 T_{max} 后，在 T_{min} 以上维持的时间：

 

t_{effective} = t_{cooling} + t_{maintaining}

 

其中：

 

- t_{cooling} ：从初始温度冷却到 T_{max} 的时间

- t_{maintaining} ：在 [T_{min}, T_{max}] 区间内的维持时间

 

保温效果评估

 

热阻网络模型：

 

R_{total} = R_{conv,out} + R_{ins} + R_{conv,in}

 

热损失率：

 

\dot{Q}_{loss} = \frac{T_{inside} - T_{outside}}{R_{total}}

 

温度回升时间常数：

 

\tau = \frac{m \cdot c_{p,eff}}{1/R_{total}}

 

优化目标函数

 

多目标优化：

 

\min \left\{ -\text{CoolingEffect}, \ \text{TimeToThreshold}, \ \text{FrostbiteRisk} \right\}

 

约束条件：

 

- 温度安全： T_{skin} \geq 2°C 或 T_{contact} \geq 0°C 

- 时间约束： t_{use} \leq t_{max\_safe} 

- 材料约束： m_{ice} \geq m_{min} ，保温层完整性

 

四、代码模块化实现

 

项目结构

 

ice_cooling_optimizer/

├── main.py # 主程序入口

├── config.py # 配置文件

├── core/ # 核心计算模块

│ ├── __init__.py

│ ├── phase_change_model.py # 相变传热模型

│ ├── heat_transfer.py # 热传递计算

│ ├── ice_melting.py # 冰融化动力学

│ └── cooling_effect.py # 降温效果评估

├── materials/ # 材料数据库

│ ├── __init__.py

│ └── material_database.py

├── optimization/ # 优化模块

│ ├── __init__.py

│ └── usage_optimizer.py

├── simulation/ # 仿真模块

│ ├── __init__.py

│ └── transport_simulator.py

├── sensors/ # 传感器接口

│ ├── __init__.py

│ └── temperature_sensor.py

├── visualization/ # 可视化模块

│ ├── __init__.py

│ └── cooling_visualizer.py

├── utils/ # 工具函数

│ ├── __init__.py

│ ├── unit_converter.py

│ └── validators.py

├── tests/ # 单元测试

├── requirements.txt

└── README.md

 

核心代码实现

 

1. config.py - 配置文件

 

"""

配置模块：存储系统全局参数和材料物理属性

基于智能与分子化学工程课程中的相变传热理论

"""

from dataclasses import dataclass, field

from typing import Dict, List, Optional, Tuple

import json

import math

 

# 物理常数

PHYSICAL_CONSTANTS = {

    "water_density": 1000.0, # 水密度 (kg/m³)

    "ice_density": 917.0, # 冰密度 (kg/m³)

    "latent_heat_fusion": 334000.0, # 冰熔化潜热 (J/kg)

    "ice_specific_heat": 2090.0, # 冰比热容 (J/(kg·K))

    "water_specific_heat": 4180.0, # 水比热容 (J/(kg·K))

    "absolute_zero_celsius": -273.15, # 绝对零度 (°C)

    "ice_melting_point": 0.0, # 冰熔点 (°C)

    "gas_constant": 8.314, # 气体常数 (J/(mol·K))

    "stefan_boltzmann": 5.67e-8, # Stefan-Boltzmann常数 (W/(m²·K⁴))

}

 

@dataclass

class MaterialProperties:

    """材料物理属性数据类"""

    # 基本属性

    name: str

    density: float # 密度 (kg/m³)

    thermal_conductivity: float # 导热系数 (W/(m·K))

    specific_heat: float # 比热容 (J/(kg·K))

    latent_heat: Optional[float] = None # 相变潜热 (J/kg)

    melting_point: Optional[float] = None # 熔点 (°C)

    

    # 热辐射属性

    emissivity: float = 0.9 # 发射率 (0-1)

    

    # 机械属性

    thickness: float = 0.01 # 厚度 (m)

    

    # 老化/衰减参数

    degradation_rate: float = 0.02 # 性能衰减率 (%/cycle)

    max_cycles: int = 100 # 最大使用循环次数

    

    # 成本参数

    cost_per_unit: float = 1.0 # 单位成本

    unit: str = "m²" # 成本计算单位

 

@dataclass

class IcePackConfiguration:

    """冰袋配置数据类"""

    # 尺寸 (m)

    length: float

    width: float

    height: float

    

    # 内部相变材料

    pcm_mass: float # 相变材料质量 (kg)

    pcm_type: str = "pure_water" # 相变材料类型

    

    # 保温层

    insulation_material: str = "foam_pe" # 保温材料

    insulation_thickness: float = 0.02 # 保温层厚度 (m)

    

    # 外部包装

    outer_layer_material: str = "hdpe" # 外壳材料

    surface_area: Optional[float] = None # 表面积 (m²)，自动计算

    

    # 使用参数

    initial_temperature: float = -18.0 # 初始温度 (°C，冷冻室温度)

    target_min_temp: float = 2.0 # 目标最低温度 (°C)

    target_max_temp: float = 15.0 # 有效降温上限 (°C)

    

    def get_volume(self) -> float:

        """计算冰袋体积 (m³)"""

        return self.length * self.width * self.height

    

    def get_surface_area(self) -> float:

        """计算外表面积 (m²) - 长方体近似"""

        l, w, h = self.length, self.width, self.height

        return 2 * (l*w + w*h + h*l)

    

    def get_heat_transfer_area(self) -> float:

        """计算有效传热面积 (m²) - 扣除保温层后"""

        if self.surface_area is None:

            self.surface_area = self.get_surface_area()

        return self.surface_area

 

@dataclass

class EnvironmentalConditions:

    """环境条件数据类"""

    ambient_temperature: float # 环境温度 (°C)

    relative_humidity: float = 50.0 # 相对湿度 (%)

    air_velocity: float = 0.1 # 空气流速 (m/s)

    external_radiation_temp: float = None # 外部辐射温度 (°C)，默认等于环境温度

    pressure: float = 101325.0 # 气压 (Pa)

    

    def get_radiation_temperature(self) -> float:

        """获取辐射温度"""

        if self.external_radiation_temp is None:

            return self.ambient_temperature

        return self.external_radiation_temp

 

@dataclass

class PatientParameters:

    """患者参数数据类"""

    skin_temperature: float = 33.0 # 皮肤温度 (°C)

    contact_area: float = 0.02 # 接触面积 (m²)

    blood_flow_rate: float = 0.001 # 局部血流速率 (kg/s)

    metabolic_heat: float = 50.0 # 代谢产热 (W)

    sensitivity_factor: float = 1.0 # 敏感因子（儿童、老人较高）

 

# 材料数据库 - 常见冰袋材料和保温材料

MATERIAL_DATABASE = {

    # 相变材料 (PCM)

    "pure_water": MaterialProperties(

        name="纯水",

        density=1000.0,

        thermal_conductivity=0.6,

        specific_heat=4180.0,

        latent_heat=334000.0,

        melting_point=0.0

    ),

    "salt_water_10pct": MaterialProperties(

        name="10%盐水",

        density=1070.0,

        thermal_conductivity=0.58,

        specific_heat=3900.0,

        latent_heat=310000.0,

        melting_point=-6.0

    ),

    "gel_water_base": MaterialProperties(

        name="水基凝胶",

        density=950.0,

        thermal_conductivity=0.45,

        specific_heat=3500.0,

        latent_heat=280000.0,

        melting_point=-2.0

    ),

    "urea_water": MaterialProperties(

        name="尿素水溶液",

        density=1050.0,

        thermal_conductivity=0.52,

        specific_heat=3700.0,

        latent_heat=295000.0,

        melting_point=-11.0

    ),

    

    # 保温材料

    "foam_pe": MaterialProperties(

        name="聚乙烯泡沫",

        density=30.0,

        thermal_conductivity=0.03,

        specific_heat=1300.0,

        thickness=0.02,

        emissivity=0.9,

        degradation_rate=0.01,

        max_cycles=200,

        cost_per_unit=0.5,

        unit="m²"

    ),

    "vacuum_insulation": MaterialProperties(

        name="真空隔热板(VIP)",

        density=200.0,

        thermal_conductivity=0.004,

        specific_heat=1000.0,

        thickness=0.008,

        emissivity=0.02,

        degradation_rate=0.005,

        max_cycles=500,

        cost_per_unit=15.0,

        unit="m²"

    ),

    "aerogel": MaterialProperties(

        name="气凝胶",

        density=100.0,

        thermal_conductivity=0.014,

        specific_heat=1100.0,

        thickness=0.005,

        emissivity=0.3,

        degradation_rate=0.008,

        max_cycles=300,

        cost_per_unit=25.0,

        unit="m²"

    ),

    "wool_felt": MaterialProperties(

        name="羊毛毡",

        density=150.0,

        thermal_conductivity=0.04,

        specific_heat=1400.0,

        thickness=0.01,

        emissivity=0.95,

        degradation_rate=0.03,

        max_cycles=50,

        cost_per_unit=2.0,

        unit="m²"

    ),

    

    # 外壳材料

    "hdpe": MaterialProperties(

        name="高密度聚乙烯",

        density=950.0,

        thermal_conductivity=0.5,

        specific_heat=1900.0,

        emissivity=0.9

    ),

    "pp": MaterialProperties(

        name="聚丙烯",

        density=900.0,

        thermal_conductivity=0.14,

        specific_heat=1900.0,

        emissivity=0.9

    ),

    "silicone": MaterialProperties(

        name="硅胶",

        density=1200.0,

        thermal_conductivity=0.2,

        specific_heat=1700.0,

        emissivity=0.92

    ),

}

 

# 预设冰袋配置

ICE_PACK_PRESETS = {

    "standard_ice_pack": IcePackConfiguration(

        length=0.25,

        width=0.18,

        height=0.03,

        pcm_mass=0.3,

        pcm_type="pure_water",

        insulation_material="foam_pe",

        insulation_thickness=0.02,

        initial_temperature=-18.0

    ),

    "large_clinical_pack": IcePackConfiguration(

        length=0.35,

        width=0.25,

        height=0.04,

        pcm_mass=0.8,

        pcm_type="gel_water_base",

        insulation_material="foam_pe",

        insulation_thickness=0.025,

        initial_temperature=-18.0

    ),

    "compact_portable": IcePackConfiguration(

        length=0.15,

        width=0.12,

        height=0.02,

        pcm_mass=0.1,

        pcm_type="salt_water_10pct",

        insulation_material="wool_felt",

        insulation_thickness=0.015,

        initial_temperature=-18.0

    ),

    "premium_vip_pack": IcePackConfiguration(

        length=0.3,

        width=0.22,

        height=0.035,

        pcm_mass=0.5,

        pcm_type="pure_water",

        insulation_material="vacuum_insulation",

        insulation_thickness=0.008,

        initial_temperature=-18.0

    ),

}

 

# 系统配置

SYSTEM_CONFIG = {

    "default_time_step": 60.0, # 默认时间步长 (s)

    "default_simulation_duration": 7200.0, # 默认仿真时长 (s，2小时)

    "numerical_tolerance": 1e-6, # 数值计算容差

    "max_iterations": 10000, # 最大迭代次数

    "safety_margin": 0.1, # 安全裕度

    "frostbite_threshold": 2.0, # 冻伤风险温度阈值 (°C)

    "freeze_damage_threshold": -5.0, # 组织冻结损伤温度 (°C)

    "output_precision": 3, # 输出精度

}

 

def load_material_database(db_path: str = "materials/material_database.json") -> Dict[str, MaterialProperties]:

    """从JSON文件加载材料数据库"""

    try:

        with open(db_path, 'r', encoding='utf-8') as f:

            raw_data = json.load(f)

            materials = {}

            for name, props in raw_data.items():

                materials[name] = MaterialProperties(**props)

            return materials

    except FileNotFoundError:

        print(f"警告：未找到材料数据库文件 {db_path}，使用内置数据库")

        return MATERIAL_DATABASE.copy()

 

def save_material_database(materials: Dict[str, MaterialProperties], db_path: str):

    """保存材料数据库到JSON文件"""

    data = {}

    for name, props in materials.items():

        data[name] = {

            'name': props.name,

            'density': props.density,

            'thermal_conductivity': props.thermal_conductivity,

            'specific_heat': props.specific_heat,

            'latent_heat': props.latent_heat,

            'melting_point': props.melting_point,

            'emissivity': props.emissivity,

            'thickness': props.thickness,

            'degradation_rate': props.degradation_rate,

            'max_cycles': props.max_cycles,

            'cost_per_unit': props.cost_per_unit,

            'unit': props.unit

        }

    with open(db_path, 'w', encoding='utf-8') as f:

        json.dump(data, f, indent=2, ensure_ascii=False)

 

def create_custom_ice_pack(

    length: float,

    width: float,

    height: float,

    pcm_mass: float,

    pcm_type: str = "pure_water",

    insulation_material: str = "foam_pe",

    insulation_thickness: float = 0.02,

    **kwargs

) -> IcePackConfiguration:

    """创建自定义冰袋配置"""

    return IcePackConfiguration(

        length=length,

        width=width,

        height=height,

        pcm_mass=pcm_mass,

        pcm_type=pcm_type,

        insulation_material=insulation_material,

        insulation_thickness=insulation_thickness,

        **kwargs

    )

 

# 单位转换工具

class UnitConverter:

    """单位转换工具类"""

    

    # 长度转换

    LENGTH_CONVERSIONS = {

        ("m", "cm"): 100,

        ("m", "mm"): 1000,

        ("m", "inch"): 39.37,

        ("cm", "m"): 0.01,

        ("mm", "m"): 0.001,

        ("inch", "m"): 0.0254,

    }

    

    # 温度转换 (相对于摄氏度的偏移)

    TEMPERATURE_OFFSETS = {

        ("C", "K"): 273.15,

        ("K", "C"): -273.15,

        ("C", "F"): lambda c: c * 9/5 + 32,

        ("F", "C"): lambda f: (f - 32) * 5/9,

    }

    

    # 功率转换

    POWER_CONVERSIONS = {

        ("W", "kW"): 0.001,

        ("kW", "W"): 1000,

        ("W", "BTU/hr"): 3.412,

        ("BTU/hr", "W"): 0.293,

    }

    

    # 能量转换

    ENERGY_CONVERSIONS = {

        ("J", "kJ"): 0.001,

        ("kJ", "J"): 1000,

        ("J", "kcal"): 0.000239,

        ("kcal", "J"): 4184,

    }

    

    @classmethod

    def convert_length(cls, value: float, from_unit: str, to_unit: str) -> float:

        """长度单位转换"""

        if from_unit == to_unit:

            return value

        key = (from_unit, to_unit)

        if key in cls.LENGTH_CONVERSIONS:

            return value * cls.LENGTH_CONVERSIONS[key]

        raise ValueError(f"不支持的长度单位转换: {from_unit} -> {to_unit}")

    

    @classmethod

    def convert_temperature(cls, value: float, from_unit: str, to_unit: str) -> float:

        """温度单位转换"""

        if from_unit == to_unit:

            return value

        key = (from_unit, to_unit)

        if key in cls.TEMPERATURE_OFFSETS:

            converter = cls.TEMPERATURE_OFFSETS[key]

            if callable(converter):

                return converter(value)

            return value + converter

        raise ValueError(f"不支持的温度单位转换: {from_unit} -> {to_unit}")

    

    @classmethod

    def convert_power(cls, value: float, from_unit: str, to_unit: str) -> float:

        """功率单位转换"""

        if from_unit == to_unit:

            return value

        key = (from_unit, to_unit)

        if key in cls.POWER_CONVERSIONS:

            return value * cls.POWER_CONVERSIONS[key]

        raise ValueError(f"不支持的功率单位转换: {from_unit} -> {to_unit}")

    

    @classmethod

    def convert_energy(cls, value: float, from_unit: str, to_unit: str) -> float:

        """能量单位转换"""

        if from_unit == to_unit:

            return value

        key = (from_unit, to_unit)

        if key in cls.ENERGY_CONVERSIONS:

            return value * cls.ENERGY_CONVERSIONS[key]

        raise ValueError(f"不支持的能量单位转换: {from_unit} -> {to_unit}")

 

2. core/phase_change_model.py - 相变传热模型

 

"""

相变传热模型模块：实现冰融化过程的数学建模

基于Stefan问题和能量守恒定律

这是智能与分子化学工程课程的核心应用

"""

import math

import numpy as np

from dataclasses import dataclass, field

from typing import Tuple, Optional, Callable, List

from enum import Enum

import logging

 

from config import PHYSICAL_CONSTANTS, MaterialProperties, IcePackConfiguration, EnvironmentalConditions

from utils.validators import validate_positive, validate_temperature_range

 

logger = logging.getLogger(__name__)

 

class PhaseState(Enum):

    """相态枚举"""

    SOLID = "solid" # 固态（冰）

    LIQUID = "liquid" # 液态（水）

    MIXED = "mixed" # 固液共存（相变中）

    SUPERCOOLED = "supercooled" # 过冷态

 

@dataclass

class PhaseChangeState:

    """相变状态数据类"""

    temperature: float # 当前温度 (°C)

    phase: PhaseState # 当前相态

    solid_fraction: float # 固体质量分数 (0-1)

    liquid_fraction: float # 液体质量分数 (0-1)

    interface_position: float = 0.0 # 固液界面位置 (m)，用于Stefan问题

    heat_flux_in: float = 0.0 # 输入热流 (W)

    heat_flux_out: float = 0.0 # 输出热流 (W)

    latent_heat_absorbed: float = 0.0 # 已吸收的潜热 (J)

    sensible_heat_change: float = 0.0 # 显热变化 (J)

 

class PhaseChangeModel:

    """

    相变传热模型核心类

    实现冰袋融化的完整物理建模

    

    理论基础：

    1. 能量守恒定律

    2. Stefan问题（移动边界问题）

    3. 傅里叶导热定律

    4. 牛顿冷却定律

    """

    

    def __init__(

        self,

        ice_pack_config: IcePackConfiguration,

        material_db: Dict[str, MaterialProperties]

    ):

        """

        初始化相变模型

        

        Args:

            ice_pack_config: 冰袋配置

            material_db: 材料数据库

        """

        self.config = ice_pack_config

        self.material_db = material_db

        

        # 获取材料属性

        self.pcm = material_db.get(ice_pack_config.pcm_type, material_db["pure_water"])

        self.insulation = material_db.get(ice_pack_config.insulation_material, material_db["foam_pe"])

        self.outer_layer = material_db.get(ice_pack_config.outer_layer_material, material_db["hdpe"])

        

        # 计算几何参数

        self.volume = ice_pack_config.get_volume()

        self.surface_area = ice_pack_config.get_surface_area()

        self.heat_transfer_area = ice_pack_config.get_heat_transfer_area()

        

        # 初始化状态

        self.state = PhaseChangeState(

            temperature=ice_pack_config.initial_temperature,

            phase=PhaseState.SOLID,

            solid_fraction=1.0,

            liquid_fraction=0.0

        )

        

        # 历史记录

        self.temperature_history: List[Tuple[float, float]] = [] # [(time, temp), ...]

        self.phase_history: List[Tuple[float, PhaseState]] = [] # [(time, phase), ...]

        self.energy_history: List[Tuple[float, float]] = [] # [(time, cumulative_energy), ...]

        

        # 数值求解参数

        self.time_step = 1.0 # 时间步长 (s)

        self.convergence_tolerance = 1e-6

        self.max_newton_iterations = 100

        

        logger.info(f"相变模型初始化完成: PCM={self.pcm.name}, 体积={self.volume:.6f}m³")

    

    def calculate_thermal_resistances(

        self,

        env_conditions: EnvironmentalConditions

    ) -> Dict[str, float]:

        """

        计算热阻网络

        

        热阻网络模型：

        Q = ΔT / (R_conv_out + R_ins + R_conv_in)

        

        Args:

            env_conditions: 环境条件

            

        Returns:

            各热阻值字典 (K/W)

        """

        # 外部对流热阻

        # h_out = 10.45 - v + 10*v^(1/2) (风速相关经验公式)

        air_velocity = env_conditions.air_velocity

        h_out = 10.45 - air_velocity + 10 * math.sqrt(air_velocity) if air_velocity > 0.1 else 10.0

        R_conv_out = 1.0 / (h_out * self.heat_transfer_area)

        

        # 保温层导热热阻

        # R = d / (k * A)

        R_ins = self.insulation.thickness / (self.insulation.thermal_conductivity * self.heat_transfer_area)

        

        # 内部对流热阻（冰水混合物到内壁）

        h_in = 500.0 # 强制对流系数（接触传热）

        R_conv_in = 1.0 / (h_in * self.heat_transfer_area * 0.1) # 假设10%面积接触

        

        # 辐射热阻（可选）

        radiation_temp = env_conditions.get_radiation_temperature()

        delta_T_rad = self.state.temperature - radiation_temp

        if delta_T_rad > 0:

            sigma = PHYSICAL_CONSTANTS["stefan_boltzmann"]

            epsilon = self.outer_layer.emissivity

            R_rad = 1.0 / (epsilon * sigma * self.surface_area * (self.state.temperature + 273.15 + radiation_temp + 273.15) * (self.state.temperature + 273.15)**2 + (radiation_temp + 273.15)**2)

        else:

            R_rad = float('inf') # 无辐射换热

        

        return {

            "R_conv_out": R_conv_out,

            "R_insulation": R_ins,

            "R_conv_in": R_conv_in,

            "R_radiation": R_rad,

            "R_total": R_conv_out + R_ins + R_conv_in + R_rad,

            "h_external": h_out

        }

    

    def calculate_heat_flux(

        self,

        env_conditions: EnvironmentalConditions,

        resistances: Dict[str, float]

    ) -> Tuple[float, Dict[str, float]]:

        """

        计算净热流

利用AI解决实际问题，如果你觉得这个工具好用，欢迎关注长安牧笛！