MOOSE框架中SolutionIC初始条件的子域限制功能增强

背景介绍

在基于MOOSE框架的多物理场仿真中，初始条件的设置对于模拟结果的准确性至关重要。SolutionIC是MOOSE中一个常用的初始条件类型，它允许用户从已有的Exodus格式结果文件中读取场变量值作为新模拟的初始条件。然而，在实际应用中，用户经常需要将源文件中的场变量值仅应用到目标网格的特定子域上，或者源文件和目标网格具有不同的子域划分。

原有功能限制

在MOOSE框架的早期版本中，SolutionIC通过block参数同时控制两个功能：

1. 目标网格中应用初始条件的子域限制
2. 源Exodus文件中读取数据的子域限制

这种设计在源文件和目标网格具有相同子域划分时工作良好，但在以下场景会出现问题：

• 源文件没有子域划分而目标网格有子域划分
• 源文件和目标网格的子域划分不同
• 用户希望从源文件的特定区域读取数据但应用到目标网格的不同区域

功能增强方案

为解决上述限制，MOOSE框架对SolutionIC进行了功能增强，新增了专门用于控制源文件子域限制的参数。具体改进包括：

1. 引入新的参数source_block，专门用于指定源Exodus文件中读取数据的子域限制
2. 保留原有block参数，仅用于控制目标网格中应用初始条件的子域限制
3. 当source_block未指定时，保持向后兼容性，使用block参数同时控制源和目标子域限制

应用示例

假设用户需要将一个无子域划分的2D模拟结果应用到具有子域划分的新网格中，且只希望在新网格的特定子域上应用初始条件：

[Mesh]
  [2d]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 10
    ny = 10
    xmax = 1
    ymax = 1
  []
  [subdomain]
    type = ParsedSubdomainMeshGenerator
    block_id = 1
    combinatorial_geometry = 'y < 0.5'
    input = 2d
  []
[]

[ICs]
  [solution_ic]
    type = SolutionIC
    from_variable = 'temperature'
    solution_uo = source_solution
    variable = temperature
    block = 1  # 仅应用到目标网格的子域1
    source_block = 'all'  # 从源文件的所有区域读取数据
  []
[]

技术实现细节

在实现层面，这一增强涉及以下关键技术点：

1. 参数解析系统：扩展了SolutionIC的参数解析能力，新增source_block参数处理逻辑
2. 数据映射机制：改进了场变量数据从源网格到目标网格的映射算法，支持不同的子域限制条件
3. 错误处理：增强了子域检查逻辑，提供更清晰的错误提示信息
4. 兼容性保障：确保旧输入文件仍能正常工作，同时支持新功能

应用场景

这一功能增强特别适用于以下仿真场景：

1. 多尺度模拟：将粗网格的整体结果作为细网格局部区域的初始条件
2. 多物理场耦合：将一个物理场的全域结果作为另一个物理场局部区域的初始条件
3. 参数化研究：基于基准案例结果初始化参数化研究的特定区域
4. 缺陷模拟：在全域健康模型基础上初始化局部缺陷区域的场变量

总结

MOOSE框架对SolutionIC初始条件的子域限制功能增强，显著提高了初始条件设置的灵活性和精确性。通过分离源文件和目标网格的子域控制，用户可以更精确地控制初始条件的应用范围，为复杂多物理场仿真提供了更强大的支持。这一改进不仅解决了现有用户的实际需求，也为MOOSE框架在更广泛领域的应用奠定了基础。