关于模拟：有几个商业求解器可能会有所帮助。

• 用于纳米光子学的一个非常流行的 FDTD Maxwell 方程求解器是Lumerical。它具有光热和液晶（我会说这是您应该考虑的另一个关键字）模拟，应该非常相关。

• 另一个首选的多物理场求解器是 COMSOL。搜索电流体动力学和磁流体动力学会得到很多结果。基于 COMSOL 建模的几篇相关论文：

  1. MG Sailer，BJ Kirby，“与微纳米级流体流动现象有关的电流体动力学系统的计算模拟”，在Proc 中。COMSOL 会议 ，马萨诸塞州波士顿，2008 年。
  2. A. Mulye, S. Potnis，“用于性能优化的电流体动力微型泵建模”，在Proc。COMSOL 会议 ，印度班加罗尔，2013 年。
  3. N. Kleinknecht，“ COMSOL 中的高频磁流体动力学计算”，Proc。COMSOL 会议 ，斯图加特，德国，2011 年。
  4. Y. Yan、S. Smolentsev、M. Abdou，“验证 COMSOL 多物理场在聚变应用中的磁流体动力学 (MHD) 流动，”在Proc。COMSOL 会议 ，马萨诸塞州波士顿，2017 年。

关于数值技术本身：

• 我最近偶然发现了以下论文（这是我写回答您问题的主要动机）：E. Veerapaneni，“三维囊泡电流体动力学的积分方程方法”，J. Comp。物理。, 卷。326，第 278-289 页，2016 年 12 月。摘自摘要：该方法使用半隐式时间步长方案来克服与控制方程相关的数值刚度。
• 几年前我看到的另一篇论文以磁流体动力学的数值研究为特色：M. Kinet、B. Knaepen 和 S. Molokov，“具有导电壁的管道中磁流体动力学流动的不稳定性和过渡”，Phys。牧师莱特。, 卷。103，154501，2009 年 10 月。