湖北智能网站建设找哪家允许发外链的网站

湖北智能网站建设找哪家,允许发外链的网站,网络推广渠道和方式,网站开发人员篡改客户数据多维物理场中张量方程的符号识别Symbolic identification of tensor equationsin multidimensional physical fields陈天逸1#xff0c;杨浩1#xff0c;马文军1#xff0c;张俊1,∗北京航空航天大学航空科学与工程学院#xff0c;北京100191引用格式#xff1a; T. Chen,…多维物理场中张量方程的符号识别Symbolic identification of tensor equationsin multidimensional physical fields陈天逸1杨浩1马文军1张俊1,∗北京航空航天大学航空科学与工程学院北京100191引用格式T. Chen, H. Yang, W. Ma, and J. Zhang, “Symbolic Identification of Tensor Equations in Multidimensional Physical Fields.” Journal of Fluid Mechanics, 1024, A34 (2025) . DOI: https://doi.org/10.1017/jfm.2025.10710.摘要近年来数据驱动方法在从模拟或实验数据中发现控制方程方面展现出巨大潜力。然而现有方法多局限于标量方程能够识别张量关系的研究屈指可数。本文提出一个通用的数据驱动框架用于识别张量方程称为张量方程符号识别方法。该方法的核心理念——采用宿主-质粒结构表征张量方程——源于多维基因表达式编程思想的启发。为提升进化过程的鲁棒性本方法采用了遗传信息保留策略。此外相较于传统进化算法本研究引入了两项关键创新首先通过维度齐次性检查来约束搜索空间排除物理无效表达式其次以张量线性回归技术替代传统线性缩放方法极大提升了数值系数优化的效率。我们通过两个基准案例验证了本方法的有效性在合成数据中精确还原了目标方程同时展现出良好的抗噪能力与灵活的表达性能。进一步地本方法被应用于直接从分子模拟数据中识别本构关系——这些数据完全无需依赖宏观本构模型即可生成。该方法同时适用于可压缩与不可压缩流动工况并成功识别出相应的宏观形式彰显了其在数据驱动的张量方程发现领域的应用前景。一、研究背景近年来从数据中自动发现控制方程已成为AI for Science领域的前沿方向。以符号回归为代表的一系列方法能够基于观测数据自动归纳其内在规律并以显式数学方程的形式加以表达为揭示物理机制、加速工程计算提供了有力支持。然而现有符号回归方法大多局限于标量方程的识别。在处理高维系统时候选函数空间会随系统变量的增多呈指数级增长导致在有限计算成本下难以稳定可靠地提取出具有物理意义的控制方程。将符号回归拓展至张量方程具有重要意义。高维物理系统的数学描述天然依赖张量结构诸多领域的核心物理量本质上均为张量如流体力学中的变形率张量、固体力学中的弹性张量等。张量不仅能刻画物理量的方向依赖性还具备旋转不变性等重要性质因而能够简洁而准确地描述高维系统的物理规律。目前已有的张量符号回归方法往往难以同时兼顾表达能力、物理一致性和搜索效率。为解决以上挑战本研究提出了一种通用的张量方程符号识别Symbolic Identification of Tensor Equations, SITE框架。二、研究方法SITE是一个面向张量方程符号识别的通用框架旨在高效挖掘具有量纲一致性的张量表达式。图1展示了该框架的整体流程主要包括数据预处理、张量与标量表达式的构建以及完整的进化优化过程。首先系统对输入数据进行预处理包括对部分物理量进行梯度计算并构建候选终端库。接着基于候选终端库和一组预定义的候选符号库SITE通过组合不同符号结构生成候选表达式并借助进化流程实现表达式的迭代优化。SITE的核心优化机制为进化算法并通过树结构来实现灵活的表达能力。具体而言框架采用多维基因表达式编程M-GEP中使用的宿主-质粒架构以实现张量与标量的协同表达同时提出基因保留策略确保标量表达式信息在变异、重组等遗传操作中不被丢失。在该架构中宿主个体表示二阶张量质粒个体表示标量二者通过宿主表达式中的特殊p算符进行关联。图1SITE框架示意图在进化算法中SITE引入了两个关键步骤以提升性能。首先在所有表达式被评估前执行量纲齐次性校核dimensional homogeneity check, DHC不满足量纲一致性的表达式将被赋予一个足够大的损失从而在演化过程中被淘汰掉。这一模块还包括一个种子注入策略如图2所示通过向种群中定期注入满足量纲齐次性约束的种子个体以实现在进化前期辅助生成合法表达式、在后期帮助种群跳出局部最优。为提高数值系数的优化效率本文引入了张量线性回归tensor linear regression, TLR技术。在张量符号回归中使用单一全局缩放系数无法有效拟合预测张量表达式的各个分量因此常见的线性放缩技术失效。TLR为同一个体中的多个子表达式分别确定独立的系数并通过张量特征的线性组合来拟合预测张量表达式与目标之间的关系。这个系数向量通过理论分析直接给出并在程序中用于快速求解。这种方法提高了系数最优化的效率并通过为不同基因赋予不同重要性实现了隐式的特征选择。图2种子注入策略示意图三、结果3.1 基准算例验证本文首先通过几个基准方程进行了SITE框架的验证测试。这里以电磁学领域的Maxwell应力张量为例验证SITE识别张量方程的能力并围绕该框架开展系列研究包括对关键模块的消融实验与SOTA方法的性能对比以及针对数据噪声的鲁棒性测试。我们按照特定的函数形式在空间中设置了电磁场并生成了稀疏的电磁场数据。Maxwell应力张量的正确形式如下TLR的消融实验结果如表1所示无论是否使用随机数值常数random numerical constants, RNCTLR都能提供足够快的系数优化效率且识别出正确的方程。与M-GEP方法的对比实验结果如表2所示可以看到无论是识别出正确方程的成功率还是算法的搜索效率SITE都优于M-GEP。最后我们在原始数据基础上添加了不同程度的高斯噪声并进行测试表3的结果表明SITE对于数值噪声具有良好的鲁棒性。3.2基于分子模拟数据的方程识别我们进一步将SITE方法应用于通过分子模拟产生的数据未预先嵌入任何宏观控制方程如NS方程或预定义本构关系。本文采用的分子模拟方法是直接模拟蒙特卡洛direct simulation Monte Carlo, DSMC方法其基于分子动理论模拟分子的运动和碰撞过程通过统计平均获得宏观物理量是更底层的数据生成方式。本研究的目标正是验证SITE能够从不施加任何宏观本构假设的分子模拟数据中挖掘出宏观张量本构关系同时评估该框架在不同流动条件下识别对应宏观方程的能力。我们采用DSMC方法模拟了全局Knudsen数为0.005的方腔流动顶盖速度分别为50m/s不可压缩情况和337m/s可压缩情况模拟气体为氩气。流场的宏观量通过对分子信息的统计平均得到另外我们训练了一个全连接神经网络并利用自动微分技术计算宏观量的梯度。我们选择400×400个空间点中靠近中心的200组数据点作为训练数据以排除边界效应的影响。我们知道连续介质假设下给出的牛顿本构方程为SITE在两种流场条件下识别得到的方程如表4所示。在低速条件下SITE发现的本构关系仅包含前述本构关系中的剪切应力项与静压项这与低速流动时可压缩效应较弱的物理预期一致而在高速条件下SITE发现的本构关系中出现了额外的可压缩项这一结果表明SITE框架能够成功分辨出与流动状态相关的物理贡献展现了其根据不同流动特性自动识别相应本构关系的能力。四、总结本研究提出了张量方程的符号识别SITE框架。该框架通过宿主-质粒架构实现张量与标量的协同运算结合量纲齐次性校核与种子注入策略确保了所得方程的物理有效性。此外通过引入张量线性回归技术提升了系数的优化效率与精度。我们在基准算例中验证了SITE框架在噪声场景下的高准确率与鲁棒性且成功将其应用于DSMC方法模拟产生的数据实现了不可压缩与可压缩工况下本构关系的自动识别。SITE框架有望为高维物理系统的可解释性建模提供一个有潜力的工具。代码与数据本研究相关代码与数据已公开于https://github.com/BUAA-MARS-group/SITE原文下载https://doi.org/10.1017/jfm.2025.10710注文章由原作者投稿分享向本公众号授权发布。