【初学者】comsol周期性电流如何设置？例如通电三天断电三天

Question

设置如下，首先用矩形波分段，再用解析函数周期化，然后在需要设置电流的地方调用了函数an1，并引入时间

最终出现重复误差测试失败，可能已达到奇点的报错。请问有能解答相关问题吗？还在学习中可能需要一点详细的回答

Answer 1

该回答引用ChatGPT-4o并结合自身经验回答，望对题主有所帮助/启发；若有帮助，还望采纳🤝🤝🤝

---

根据您的描述和上传的图片信息，您在COMSOL中设置了一个周期性电流，通过矩形波函数分段并周期化电流来实现三天通电三天断电的循环，但遇到了“重复误差测试失败，可能已达到奇点”的报错。这种报错通常在设置周期性边界条件或者间歇性加载时会发生，尤其是在时间步长或边界条件设置不合适的情况下。

问题分析
时间步长设置问题：周期性或分段函数在不同时间段内有显著变化，这可能会导致数值解算时无法稳定收敛，尤其是在每次周期变化的边界点。如果时间步长选择不合适，会影响收敛性。
函数设置错误：如果在定义周期性电流的函数（如an1(t)）时，周期或者函数表达式有错误，也会影响模拟的稳定性。
电流幅值：如果电流幅值设得过大，尤其是周期性切换的过程中，可能会导致数值计算困难，从而产生误差。
解决方案
您可以从以下几个方面调整和检查设置：

1. 调整时间步长
   在研究设置中，将输出时间步长改为更细的步长，以便COMSOL在电流切换的边界时刻获得更精确的解。例如可以改为range(0, 1, 14)，让时间步为1天，甚至可以更小，以确保在通断电切换时有足够的时间点来捕捉变化。
2. 检查周期性函数an1(t)的定义
   检查您定义的解析函数是否正确，例如如果是矩形波的定义，确保其表达式定义符合“三天通电，三天断电”的周期要求。
   另外，确保周期性函数在定义中没有极端的跃变，比如可以考虑在切换时段的边界进行光滑处理，避免突变。
   您可以参考如下的定义方式（假设COMSOL支持这种表达方式）：

an1(t) = if(mod(t,6) < 3, 5e-6, 0)

mod(t,6)表示以6天为周期，if(mod(t,6) < 3, 5e-6, 0)表示前3天电流密度为5e-6，后3天为0，实现周期性变化。
3. 调整求解器设置
进入时间步长设置，尝试将时间步长设置为“严格”或“自由”，让COMSOL自动调整时间步长以适应突变的电流。
也可以在“高级”设置中增加“初始步长的最小值”和“最大值”，比如设置最小时间步长为0.1天，这样COMSOL可以在遇到奇点时自动减小步长以通过求解。
4. 尝试降低电流密度
如果电流密度设置过高，可以尝试将其暂时降低，以确认是否是因为电流过大导致数值计算困难。
5. 检查物理场设置
确认物理场（例如电流、热传导等）是否正确地联立设置，尤其是涉及到的边界条件要与时间依赖设置相符。