一种双组分抛物线型系统跨尺度时空扩散的构建方法技术方案

本发明专利技术公开了一种双组分抛物线型系统跨尺度时空扩散的构建方法，包括以下步骤：首先，建立

【技术实现步骤摘要】
一种双组分抛物线型系统跨尺度时空扩散的构建方法


[0001]本专利技术涉及仿真
，具体的说是涉及一种双组分抛物线型系统跨尺度时空扩散的构建方法
。

技术介绍

[0002]自微分方程理论诞生以来，偏微分方程一直发挥着重要的作用，广泛地出现在气象预报
、
海洋流动
、
环境治理
、
量子碰撞
、
航天计算
、
武器研究
、
飞机制造
、
交通涉及
、
水坝建设
、
油田开发
、
生物医学
、
股票期权等邻域，用于描述
、
解释或者预见各种自然现象
、
社会现象和科学工程问题
。
偏微分方程具有三大类型，分别包括抛物型
、
双曲型和椭圆型，对偏微分方程的深入研究，有助于我们更好地理解自然和社会现象中的各种变化规律
。
抛物型偏微分方程是理论中一个重要的分支，它描述了许多自然现象和工程问题中的物理过程，如反应扩散
、
热传导
、
物质运输等
。
近些年，反应扩散系统作为斑图生成系统引起了人们的极大兴趣
。
在数学上，反应扩散系统采用半线性抛物线型偏微分方程的形式，方程的解表现出广泛的行为，包括行波和类波现象的组成，以及其它自组织模式，如点状斑图
、
条形斑图
、
混合状斑图或更复杂的结构
。
这些奇妙的现象统称为“Turing
斑图”。
[0003]因此说，偏微分方程的求解技术是备受关注的研究内容
。
传统的求解技术需要精确地解析表达问题的通解，或者某些特定情形的特解
。
然而，大部分偏微分并不具备解析表达式
。
因此，利用解析表达式计算大量的时空位置的函数值，并不是高效的解决方案
。
随着计算机技术的蓬勃发展，简单便捷地数值求解越来越受到广泛学者的关注
。
相较于最简单的单组分反应扩散方程，双组分系统比单组分系统允许更大范围可能的丰富动力学行为
。
此外，大多数反应扩散系统的位置信息都是一维的
。
当从一个空间维度转向多个空间维度时，是否依然可以用来自一维系统的许多结论陈述，这是一个问题
。
故本专利设计一种双组分抛物型系统跨尺度时空扩散的构建方法
,
建立起抛物型系统在不同空间尺度下的时空动力学行为之间的关联
。

技术实现思路

[0004]本专利技术所要解决的技术问题在于充实延伸了偏微分方程建模与仿真的研究领域，提出一种双组分抛物线型系统跨尺度时空扩散的构建方法，探究双组分抛物型系统更加丰富的动力学行为，建立抛物型系统在不同空间尺度下的时空动力学行为之间的关联
。
[0005]本专利技术是通过以下技术方案来实现的，一种双组分抛物线型系统跨尺度时空扩散的构建方法，所述方法包括如下步骤：
[0006]步骤1：建立一类经典的双组分抛物线型
Gray
‑
Scott
系统，并针对不同的空间位置信息，对系统进行等距时空网格划分；
[0007]步骤2：利用
Taylor
展开定理，对步骤1中所述的系统在不同空间尺度下的偏微分方程进行离散化处理，将依赖续变量的微分方程定解问题转化为相应的依赖离散变量的代数问题；
[0008]步骤3：根据系统在不同空间尺度下的偏微分方程的初始条件和边界条件，为相应的等距时空网格点附上初值；
[0009]步骤4：选取合适的时空网格比和扩散系数，利用
MATLAB
仿真软件，分别对不同空间尺度影响下的
Gray
‑
Scott
系统进行数值分析与算法模拟，构建抛物性系统在不同空间尺度下时空动力学行为之间的关联
。
[0010]优选的，步骤1所述的一类经典的双组分抛物线型
Gray
‑
Scott
系统，其数学表达式为：
[0011][0012][0013]其中，表示拉普拉斯算子，
u
代表基底物，
v
代表催化剂，
u
和
v
代表两种化学物质，
D
u
和
D
v
分别表示两种化学物质扩散系数，
F
代表物质
u
的进给率，
K
代表物质
v
的去除率
。
[0014]优选的，步骤1所述的针对不同的空间位置信息，对系统进行等距时空网格划分，具体为：对于一维空间求解区域
Ω
∈R
，数值分析选取等距线段，对于二维空间求解区域
Ω
∈R2，数值分析选取规则矩形，对于三维空间求解区域
Ω
∈R3，数值分析选取规则立方体，且不同尺度的空间网格步长划分一致
。
[0015]优选的，步骤2所述的将依赖续变量的微分方程定解问题转化为相应的依赖离散变量的代数问题，具体为：利用
Newton
差商理论和
Taylor
展开公式，对系统在不同空间尺度下的偏微分方程进行离散化处理，时间尺度采用一阶向前差商，空间尺度采用二阶中心差商
。
[0016]优选的，步骤3所述的根据系统在不同空间尺度下的偏微分方程的初始条件和边界条件，为相应的等距时空网格点附上初值，具体为：对于一维空间扩散的初始条件，算法仿真模拟选取的初值为正弦函数，对于二维及三维空间扩散的初始条件，算法仿真模拟选取的初值为局部网格常数赋值；以上不同空间尺度下的系统的边界条件，均选取与外界零通量的
Neumann
边界
。
[0017]优选的，步骤4所述的选取合适的时空网格比和扩散系数，利用
MATLAB
仿真软件，分别对不同空间尺度影响下的
Gray
‑
Scott
系统进行数值分析与算法模拟，具体为：利用离散最大模原理，即差分格式有凸组合系数结构
。
[0018]进一步地，依据
Gray
‑
Scott
系统的空间维度对计算区域进行划分，当系统的空间尺度为一维时，其数学表达式为：
[0019][0020][0021]x
＝
Ω
x
∈[0,L],
[0022]对求解区域
Ω
x
进行以下离散化处理：
[0023][0024]其中，
h
称为空间步长，
τ
称为时间步长，
L
为空间轴长度，
T本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种双组分抛物线型系统跨尺度时空扩散的构建方法，其特征在于：所述方法包括如下步骤：步骤1：建立一类经典的双组分抛物线型
Gray
‑
Scott
系统，并针对不同的空间位置信息，对系统进行等距时空网格划分；步骤2：利用
Taylor
展开定理，对步骤1中所述的系统在不同空间尺度下的偏微分方程进行离散化处理，将依赖续变量的微分方程定解问题转化为相应的依赖离散变量的代数问题；步骤3：根据系统在不同空间尺度下的偏微分方程的初始条件和边界条件，为相应的等距时空网格点附上初值；步骤4：选取合适的时空网格比和扩散系数，利用
MATLAB
仿真软件，分别对不同空间尺度影响下的
Gray
‑
Scott
系统进行数值分析与算法模拟，构建抛物性系统在不同空间尺度下时空动力学行为之间的关联
。2.
根据权利要求1所述的双组分抛物线型系统跨尺度时空扩散的构建方法，其特征在于：步骤1所述的一类经典的双组分抛物线型
Gray
‑
Scott
系统，其数学表达式为：系统，其数学表达式为：其中，表示拉普拉斯算子，
u
代表基底物，
v
代表催化剂，
u
和
v
代表两种化学物质，
D
u
和
D
v
分别表示两种化学物质扩散系数，
F
代表物质
u
的进给率，
K
代表物质
v
的去除率
。3.
根据权利要求2所述的双组分抛物线型系统跨尺度时空扩散的构建方法，其特征在于：步骤1所述的针对不同的空间位置信息，对系统进行等距时空网格划分，具体为：对于一维空间求解区域
Ω
∈R
，数值分析选取等距线段，对于二维空间求解区域
Ω
∈R2，数值分析选取规则矩形，对于三维空间求解区域
Ω
∈R3，数值分析选取规则立方体，且不同尺度的空间网格步长划分一致
。4.
根据权利要求2所述的双组分抛物线型系统跨尺度时空扩散的构建方法，其特征在于：步骤2所述的将依赖续变量的微分方程定解问题转化为相应的依赖离散变量的代数问题，具体为：利用
Newton
差商理论和
Taylor
展开公式，对系统在不同空间尺度下的偏微分方程进行离散化处理，时间尺度采用一阶向前差商，空间尺度采用二阶中心差商
。5.
根据权利要求2所述的双组分抛物线型系统跨尺度时空扩散的构建方法，其特征在于：步骤3所述的根据系统在不同空间尺度下的偏微分方程的初始条件和边界条件，为相应的等距时空网格点附上初值，具体为：对于一维空间扩散的初始条件，算法仿真模拟选取的初值为正弦函数，对于二维及三维空间扩散的初始条件，算法仿真模拟选取的初值为局部网格常数赋值；以上不同空间尺度下的系统的边界条件，均选取与外界零通量的
Neumann
边界
。6.
根据权利要求2所述的双组分抛物线型系统跨尺度时空扩散的构建方法，其特征在于：步骤4所述的选取合适的时空网格比和扩散系数，利用
MATLAB
仿真软件，分别对不同空间尺度影响下的
Gray
‑
Scott
系统进行数值分析与算法模拟，具体为：利用离散最大模原理，即差分格式有凸组合系数结构
。
7.
根据权利要求2所述的双组分抛物线型系统跨尺度时空扩散的构建方法，其特征在于：当系统的空间尺度为一维时，其数学表达式为：于：当系统的空间尺度为一维时，其数学表达式为：
x
＝
Ω
x
∈[0,L],
对求解区域
Ω
x
进行以下离散化处理：其中，
h
称为空间步长，
τ
称为时间步长，
L
为空间轴长度，
T
为时间轴长度
,L
以
h
为步长等分为
M
份，
T
以
τ
为步长等分为
N
份；将真解
u(x,t)
限制在
(M+1)
×
(N+1)
离散网格上，相应的离散数据集合：取网格的空间步长为
h
，时间步长为
τ
，在
(x
i
,t
l
)
处，利用
Newton
差商理论和
Taylor
展开公式，时间导数离散为一阶向前差商，空间导数离散为二阶中心差商，可得：公式，时间导数离散为一阶向前差商，空间导数离散为二阶中心差商，可得：其中
i
＝
2,3
…
M,l
＝
1,2N,
略去无穷小量，进一步得到位置...

【专利技术属性】
技术研发人员：肖敏，陆云翔，
申请(专利权)人：南京邮电大学，
类型：发明
国别省市：