【技术实现步骤摘要】
本专利技术属于数值模拟
，涉及一种基于离散分数阶差分的反常扩散模拟方 法。
技术介绍
扩散方程是常用来描述自然界传热、传质现象。对于多孔介质如织物、煤岩体、土 壤等而言，流体常体现出反常扩散的特点。分数阶导数由于具有遗传特性，和反常扩散特征 吻合。最近分数阶导数被用于修正扩散方程，不少研究学者提出了时间分数阶的扩散方程[1_4]，用来描述流体在多孔介质中扩散的复杂性，避免了对于多孔介质复杂结构的刻画。 然而，分数阶导数算子的遗传特性导致了反常扩散的长历程或长时间计算变得较 为困难。寻找高精度的数值模拟方法或新的分数阶导数建模方法成为解决的方法。最近， Atici、T〇rreS[5,6]等人提出了时标上的分数阶差分。本专利技术基于该离散分数阶差分方法， 首次提出了一种新的离散分数阶扩散方程。和利用经典的分数阶扩散方程相比，该扩散模 型在保持了记忆效应的同时，表达式和数值模拟更为简洁，便于工程应用。
技术实现思路
本专利技术的目的在于提供，解决了 经典分数阶导数算子的遗传特性导致了反常扩散的长历程或长时间计算困难的问题。 本专利技术所采用的技术方案是按照以下步骤进行： 步骤1 :采用描述扩散现象的一维经典扩散方程： Ut (X，t) = Kuxx (X，t) (1) K为扩散系数，u(x，t)表示随时间和空间浓度函数，Ut (x，t)表示浓度对时间的偏 导数，uxx(x，t)表示浓度对空间的二阶偏导数，a为初始时刻，L为介质长度，方程（1)的初 边界条件为， u (x, a) = f (...
【技术保护点】
一种基于离散分数阶差分的反常扩散模拟方法，其特征在于按照以下步骤进行：步骤1：采用描述扩散现象的一维经典扩散方程：ut(x,t)＝Kuxx(x,t) (1)K为扩散系数，u(x,t)表示随时间和空间浓度函数，ut(x,t)表示浓度对时间的偏导数，uxx(x,t)表示浓度对空间的二阶偏导数，a为初始时刻，L为介质长度，方程(1)的初边界条件为，u(x,a)＝f(x),u(0,t)＝φ(t),u(L,t)＝ψ(t)f(x)表示初始时刻浓度空间各点的浓度分布，φ(t)和ψ(t)表示介质二端浓度随时间的变化；步骤2：利用离散分数阶差分：aCu(t):=hΓ(1-v)Σs=ahth-(1-v)(t-σ(sh))h(-v)Δu(sh),t∈(hN)a+(1-v)h---(2)]]>h为时间方向上的差分步长，v为分数阶差分阶数，介于0到1，为分数阶差分算子，(hN)a+(1‑ν)h为离散集合{a+(1‑ν)h,a+(2‑ν)h,...}，s,t为其中元素，u(t)为定义在该离散集合上的函数，函数σ和分别定义为：σ(sh)＝(s+1)h...
【技术特征摘要】
1. 一种基于离散分数阶差分的反常扩散模拟方法，其特征在于按照以下步骤进行： 步骤1 :采用描述扩散现象的一维经典扩散方程： Ut (X，t) = Kuxx (X，t) (1) K为扩散系数，u(x，t)表示随时间和空间浓度函数，Ut (X，t)表示浓度对时间的偏导数， uxx(x，t)表示浓度对空间的二阶偏导数，a为初始时刻，L为介质长度，方程（1)的初边界条 件为， u (x, a) = f(x)，u (0, t) = Φ (t)，u (L，t) = Ψ (t) f(x)表示初始时刻浓度空间各点的浓度分布，Φ (t)和Ψ (t)表示介质二端浓度随时 间的变化； 步骤2:利用离散分数阶差分：h为时间方向上的差分步长，V为分数阶差分阶数，介于0到1，f <为分数阶差分算子， (hN)a+(1_v)h为离散集合{a+(l_ v)h, a+(2_ v)h, . . . }, s, t为其中元素,u(t)为定义在该离 散集合上的函数，函数σ和C1:分别定义为： 〇 (sh) = (s+l)h步骤3 :利用离散化连续方程方法，利用方程（2)将方程（1)离散化为：方程（3)初边界条...
【专利技术属性】
技术研发人员：吴国成，曾生达，杜米特鲁·伯莱亚努，黄蓝蓝，
申请(专利权)人：内江师范学院，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。