基于物理信息的神经网络(Physics-informed Neural Networks;PINNs)Part-1(简单介绍)

【摘要】 基于物理信息的神经网络(Physics-informed Neural Networks,简称PINNs),是一类用于解决有监督学习任务的神经网络,它不仅能够像传统神经网络一样学习到训练数据样本的分布规律,而且能够学习到数学方程描述的物理定律。与纯数据驱动的神经网络学习相比,PINN在训练过程中施加了物理信息约束,因而能用更少的数据样本学习到更具泛化能力的模型。本文主要解析这种神经网络以及相关应用。

1 论文简介

Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations
https://www.sciencedirect.com/science/article/pii/S0021999118307125

2019年,来自布朗大学应用数学的研究团队提出了一种用物理方程作为运算限制的“物理激发的神经网络” (PINNs) 并发表在了计算物理学领域权威杂志《计算物理学期刊》(Journal of Computational Physics) 上。这篇论文一经发表就获得了大量关注。这篇论文因为代码体系的完整性使得开发人员们很容易上手把相关的学习框架应用到不同领域上去。所以在发表不久之后,一系列不同的PINNs也被其他研究者开发出来。甚至可以不夸张的说,PINN是目前AI物理领域论文中最常见到的框架和词汇之一。

在这里插入图片描述

这张图源自?

求解过程

给定一个非线性PDE
在这里插入图片描述
其中 u u u 是要求的解, N N N 是非线性偏微分算子, λ \lambda λ 是需要待定的参数。 用一个具体的例子(Burgers方程)来说明主要的想法和步骤:

简单介绍 Burgers方程:
伯格斯方程(Burgers equation) 是一个模拟冲击波的传播和反射的非线性偏微分方程。
Burgers方程是一种常见的非线性对流偏微分方程,它被广泛应用于流体力学、声学、地震学等领域的建模和仿真。Burgers方程可以写成如下形式:
∂ u ∂ t + u ∂ u ∂ x = ν ∂ 2 u ∂ x 2 \frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}=\nu\frac{\partial^2u}{\partial x^2} tu+uxu=νx22u
其中, u u u是速度场, t t t是时间, x x x是空间位置, ν \nu ν是粘性系数。Burgers方程的左边表示对流项,右边表示扩散项,它们共同描述了速度场的变化规律。
Burgers方程的解析解通常很难求得,因此常常需要使用数值方法进行求解。传统的数值方法包括有限差分法、有限元法等,它们通常需要离散化空间和时间,并且存在精度和稳定性等方面的问题。
参考链接:Burgers方程 - 百度百科

用一个具体的例子(Burgers方程)来说明主要的想法和步骤:
在这里插入图片描述

方程1为x空间上的范围为-1到1,时间上的范围为0到1;
约束条件: u u u在 [0,x] 上的方程为 − s i n ( π x ) -sin(\pi x) sin(πx);这里的0表示t,x表示空间
u u u在 [t,-1] & [t,1]上的方程为0;表示在t时刻,在x为-1或者1时,为0

利用 Neural Networks 来逼近 。 定义损失函数为
在这里插入图片描述
其中:
在这里插入图片描述

我们来分析一下损失函数。 如果神经网络能很好地求解出PDE的解,那么对于来自初边值的任一个点,其值 M S E u MSE_{u} MSEu 趋于零;对于内部的配置点,因为很好地拟合了微分方程, M S E f MSE_{f} MSEf趋于零,也就是说,损失函数的值为 0 时,我们便可以说在训练集上每个点都有神经网络的值趋于真实值 。这样便将问题转化为如何优化损失函数。利用神经网络的反向传播机制L-BFGS便可以求解。

在这里插入图片描述

利用神经网络的反向传播机制和L-BFGS便可以求解。 作者使用了tensorflow1.13 版本来求解。
参考链接:基本模型 PINNs : Physics Informed Neural Networks

相关应用论文

Physics-Informed Neural Network Super Resolution for Advection-Diffusion Models

物理信息神经网络超分辨率方法用于平流扩散模型

https://arxiv.org/abs/2011.02519

在这里插入图片描述

图1:PINNSR的网络结构。输入的低分辨率图像LR是通过在粗网格上进行模拟生成的,而不是从高分辨率图像HR进行下采样。随机像素(白色)被删除以模拟缺失像素。总损失是像素损失 L p i x L_{pix} Lpix 和物理一致性损失 L p h y s L_{phys} Lphys 的加权和。

电力系统动力学求解

Physics-Informed Neural Networks for Power Systems
PINNs for Power Systems

https://ieeexplore.ieee.org/abstract/document/9282004

在这里插入图片描述

金属表面声波裂缝检测

PHYSICS-INFORMED NEURAL NETWORK FOR ULTRASOUND NONDESTRUCTIVE QUANTIFICATION OF SURFACE BREAKING CRACKS
PINNs for 超声无损评估表面裂纹的数量化

https://link.springer.com/article/10.1007/s10921-020-00705-1

在这里插入图片描述

Reference

物理神经网络(PINN)解读

https://blog.csdn.net/waitingwinter/article/details/103979339

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/26428.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

UFS 2 -UFS架构简介2

UFS 2 -UFS架构简介2 1 UFS架构简介1.1 System Boot and Enumeration1.2 UFS Interconnect (UIC) Layer1.2.1 UFS Physical Layer Signals1.2.2 MIPI UniPro1.2.3 MIPI UniPro Related Attributes 1.3 UFS Transport Protocol (UTP) Layer1.3.1 Architectural Model1.3.1.1 Cli…

图解max{X,Y}和min{X,Y}并求相关概率

图解max{X,Y}和min{X,Y}并求相关概率 对max{X,Y}或min{X,Y}进行分解再求解 P ( m a x { X , Y } ≥ c ) P [ ( X ≥ c ) ∪ ( Y ≥ c ) ] P ( m a x { X , Y } ≤ c ) P [ ( X ≤ c ) ∩ ( Y ≤ c ) ] P ( m i n { X , Y } ≥ c ) P [ ( X ≥ c ) ∩ ( Y ≥ c ) ] P ( m i…

k8s功能优势应用场景介绍(一)

一,K8S功能: 1、数据卷 pod中容器之间共享数据,可以使用数据卷 2、应用程序健康检查 容器内服务可能进程阻塞无法处理请求,可以设置监控检查策略保证应用健壮性 3、复制应用程序实例 控制器维护着pod副本数量,保证一个pod或一组同…

C++11 auto类型推导

1.类型推导 C11引入了auto 和 decltype 关键字实现类型推导,通过这两个关键字不仅能方便地获取复杂的类型,而且还能简化书写,提高编码效率。 auto 类型推导的语法和规则 在之前的 C 版本中,auto 关键字用来指明变量的存储类型…

Allure测试报告定制全攻略,优化你的Web自动化测试框架!

目录 前言: 1. Allure测试报告简介 2. Web自动化测试框架简介 3. 封装Web自动化框架 3.1 安装Selenium 3.2 封装Selenium 3.3 定制Allure测试报告 3.3.1 适配翻译插件 3.3.2 定制测试报告样式 4. 示例代码 5. 总结 前言: 随着现在Web应用的普…

SciencePub学术 | 计算机科学类重点SCIEI征稿中

SciencePub学术刊源推荐: 计算机科学类SCI&EI征稿中!录用率高,自引率低,进展顺利。信息如下,录满为止: 一、期刊概况: 【期刊简介】IF:4.0-4.5↑, JCR 2区,中科院3区…

SpringAOP简介及实现(包含切面、切点、连接点和通知)

目录 1.什么是AOP、SpringAOP? 2.AOP的组成 3.SpringAOP的实现 4.切点的表达式 1.什么是AOP、SpringAOP? 在学习SpringAOP之前,我们得先了解一下什么是AOP。AOP是一种面向切面编程的思想。那什么是切面呢?它其实是对某一类事情…

【HR专用】Vue+SpringBoot,实现人才招聘库的开发(后端部分)

人才招聘库是企业用于储存和管理潜在候选人信息的数据库。通常情况下,这些候选人可能已经应聘过公司的职位,也可能是通过其他途径获取的,例如社交网络、招聘网站等。 对于一个中小公司来说,人力资源部绝对是一个重要部门&#xff…

测试类型(单元、集成、系统或手动测试)

测试类型(单元、集成、系统或手动测试) 单元测试 单元是系统的单个组件,例如类或单个方法。孤立地测试单元称为单元测试。 优点:速度快/易控/易写 缺点:缺乏现实性/无法捕获所有错误(例如与其他组件或服务的交互) 单元…

Arthas-JVM相关命令使用

tip:作为程序员一定学习编程之道,一定要对代码的编写有追求,不能实现就完事了。我们应该让自己写的代码更加优雅,即使这会费时费力。 开头: 我们先说下生产使用频率较高的有哪些:dashboard、heapdump、jvm…

【连续介质力学】二阶张量的图像表示

二阶张量在特定方向的投影 法向和切向分量 二阶张量T投影到 n ^ \hat n n^方向的结果是 t ⃗ ( n ^ ) T ⋅ n ^ \vec t^{(\hat n)}T \cdot \hat n t (n^)T⋅n^,其中 t ⃗ ( n ^ ) \vec t^{(\hat n)} t (n^)可以分解成: t ⃗ ( n ^ ) T ⃗ N T ⃗ S…

2023年上半年系统规划与管理师上午真题及答案解析

1.香农用概率来定量描述信息的公式如下,其中H(x)表示X的( ),Pi是( )出现第i种状态的( )。 A.信息熵 事件 概率 B.总熵 单位 概率 C.信息熵 单位 概率 D.总熵 单位 度量 2.信息传输模型中,( )负责信息的向外传播,( )负责…

VSLAM视觉里程计总结

相机模型是理解视觉里程计之前的基础。视觉里程计(VIO)主要分为特征法和直接法。如果说特征点法关注的是像素的位置差,那么,直接法关注的则是像素的颜色差。特征点法通常会把图像抽象成特征点的集合,然后去缩小特征点之…

Android 应用快捷ShortcutManager与ShortcutManagerCompat详解与实战(二)

一、介绍 之前我已通过一篇文章介绍了应用快捷的接入与Demo。如果还未看过上一篇的文章可以先了解入门。 传送门:Android 应用快捷(shortcut)功能的详解(一)_蜗牛、Z的博客-CSDN博客 有创建自然就会有管理,否则一个完美的方案不应该这么被推荐出来。如何…

Zemax Lumerical | 二维光栅出瞳扩展系统优化(下)

简介 本文提出并演示了一种以二维光栅耦出的光瞳扩展(EPE)系统优化和公差分析的仿真方法。 在这个工作流程中,我们将使用3个软件进行不同的工作 ,以实现优化系统的大目标。首先,我们使用 Lumerical 构建光栅模型并使用…

数据库实验报告--安全性实验

一、 实验目的 (1)理解SQL Server验证用户身份的过程,掌握设置身份验证模式的方法。 (2)理解登录账号的概念,掌握混合认证模式下登录账号的建立与取消方法。 (3)掌握混合认证模式…

Centos7中mysql安装配置

前提:先关闭防火墙或开启tcp的3306端口 1、查看服务器上是否有现成的安装包 yum list mysql* 2、去mysql官网的yum资源库找到对应的rpm文件的下载链接 确定系统版本 cat /etc/redhat-release 到mysql官网复制对应版本的资源下载链接 MySQL :: Download MySQL Yum…

git Husky

虽然我们已经要求项目使用eslint了,但是不能保证组员提交代码之前都将eslint中的问题解决掉了: 也就是我们希望保证代码仓库中的代码都是符合eslint规范的; 那么我们需要在组员执行 git commit 命令的时候对其进行校验,如果不符合…

说说验证码功能的实现

前言 大家好,我是 god23bin,今天说说验证码功能的实现,相信大家都经常接触到验证码的,毕竟平时上网也能遇到各种验证码,需要我们输入验证码进行验证我们是人类,而不是机器人。 验证码有多种类型&#xff…

项目中使用es(一):使用springboot操作elasticsearch

使用springboot操作es 写在前面搭建项目环境和选择合适版本具体的代码实现(1)继承ProductInfoRepository具体的代码实现(2)使用ElasticsearchRestTemplate操作问题总结最后放个demo 写在前面 对于elasticsearch的搭建&#xff0c…