二维玻尔兹曼方程的不可压缩极限

合集下载

流体力学习题课 (4)

流体力学习题课 (4)

(3)
∂u =0 ∂x
∂v =0 ∂y
∂w =0 ∂z
满足不可压缩流体连续性方程
习题四
4. 二维、定常不可压缩流动,x 方向的速度分量为 二维、定常不可压缩流动,
u = e cosh y +1
求 y 方向的速度分量 v 。 已知 y = 0 时 v = 0。 。 [解] 不可压缩流体的 解 连续性方程: 连续性方程:
条件
v = byz
kyzt − kxzt2 + kz(xt2 − yt) = 0
0≡0
无条件满足
习题四
6. 假定流管形状不随时间变化,设A为流管的横断面积,且在 断 假定流管形状不随时间变化, 为流管的横断面积, 为流管的横断面积 且在A断 面上的流动物理量是均匀的。试证明连续方程具有下述形式: 面上的流动物理量是均匀的。试证明连续方程具有下述形式:
u = 2x2 + y v = 2y2 + z w = −4(x + y)z + xy 2xyz y (x 2 − y 2 )z u=− 2 2 2 v = 2 w= 2 2 2 (x + y ) x + y2 (x + y )
u = yzt
v = xzt
w= xyt
满足不可压缩流 体连续性方程
[解] 考察是否满足不可压缩流体的连续性方程: 解 考察是否满足不可压缩流体的连续性方程: (1)
−x
∂u ∂v + =0 ∂x ∂y − e−x cosh y +
已知
∂v −x = e cosh y v = e−x sinh y + vc (x) ∂y
∂v =0 ∂y
y =0 v =0

二维不可压缩navier-stokes-landau-lifshitz方程组的整体强解

二维不可压缩navier-stokes-landau-lifshitz方程组的整体强解

二维不可压缩navier-stokes-landau-lifshitz方程组的整体强解Navier-Stokes-Landau-Lifshitz方程组是用来描述二维不可压缩流体运动的普遍方程。

它是对流体在动力、热力、物质等多种物理效应作用下进行多尺度模拟的有效工具,是研究复杂流体问题的关键。

1.Naviar-Stokes-Landau-Lifshitz方程组的基本形式Navier-Stokes-Landau-Lifshitz方程组的基本形式包含了方程的结构定义、空间变量、时间变量和输入变量,组成如下:(1)流场方程:对密度,速度和压力的描述;(2)能量方程:描述传热过程;(3)物质守恒方程:将粒子的变量连同流量和能量一起涵盖;(4)边界条件:将流体运动于受定义空间或介质内,并约束方程组的解。

2. Navier-Stokes-Landau-Lifshitz方程组的解Navier-Stokes-Landau-Lifshitz方程组的整体强解是对流体运动过程的完整描述,它包括流体的结构、动力学、热力学等理论模型。

该方程组可以使用多种空间和时间分解技术解决,比如:(1)特征定积分分解技术:特征定积分计算方法,通过积分可以得到流体的历史数据,从而求出解;(2)Galerkin有限元分解技术:通过Galerkin有限元分解方法可以得到很小的解,这也是一种将初值和边界条件一起求解的方法;(3)局部分解技术:通过局部分解计算方法可以得到相对准确的解,这是计算复杂性处理较低的一种解法;(4)关联循环求解器:通过关联循环求解器就可以得到Navier-Stokes-Landau-Lifshitz方程组的解,同时也可以求解空间多尺度的复杂流体问题。

3. Navier-Stokes-Landau-Lifshitz方程组的应用Navier-Stokes-Landau-Lifshitz方程组的本质就是对二维不可压缩流体运动过程作出数学描述,应用适用于众多流体科学领域,主要包括:(1)宇宙飞行器:在设计宇宙飞行器时,舱壁的低速、高压、非定常流体流动一直是设计中重要的一环;(2)津浦发电厂:津浦发电厂是典型的斜坡发电厂,通过模拟流体在津浦斜坡水轮发电机间的流动,可以有效提升发电效率;(3)空心叶轮压气机:空心叶轮压气机的设计要求考虑到流体的动力学特性,流场方程则可以作为压气机设计的主要辅助工具;(4)船舶航行模拟:船舶在水域的航行特别是汽轮船的航行模拟,都可以使用Navier-Stokes-Landau-Lifshitz方程组进行研究。

高等流体力学讲义二维势流

高等流体力学讲义二维势流
u = 0 Φ = 0
在不可压缩流体条件下Φ满足拉普拉斯方程
势流基本方程组
2Φ = 0 Φ + p + 1 Φ Φ + gz = f(t) t ρ 2
边界条件
在静止固壁上 ,
Φ = 0 n
无穷远处, r , u u
势流方程组与一般理想不可压缩流动方程组相比在数学上有了较大旳简化:
•后者有四个方程,而前者只有两个方程。
ln
z
-
z0
点汇
以-m 替代 m 就得到点汇旳复位势,
F(z) -m ln z 2π

F( z )
-m 2π
ln
z
-
z0
4.4 点源(汇)和点涡
点涡:势函数 流函数
F(z) ic ln z ic ln(Reiθ )
cθ ic ln R
Φ = c θ Ψ = - c ln R 等势线 c , 从圆点出发旳射线族; 流线 R=c, 同心圆族。
点源: 速度场
4.4 点源(汇)和点涡
W(z) =
dF dz
=
c z
=
c R
e-iθ
=
uR
-i

e-iθ
uR
=
c R
uθ = 0
可看作在原点有一点源释放流体向四面均匀流出,速度只有R方向分量,离 开原点愈远速度愈小。根据连续方程,经过每个同心圆旳流体流量相等。
原点是奇点,速度无穷大 R 0, uR
F(z)=Φ+ iψ
z= x + i y F(z) 旳实数部分是速度势函数Φ,虚数部分是流函数Ψ。 Φ,Ψ 满足柯西-黎曼条件,根据复变函数理论,F(Z) 是解析函数。

具有minkowski型分形边界的二维波动方程

具有minkowski型分形边界的二维波动方程

具有minkowski型分形边界的二维波动方程分形是一种具有自相似性的几何形状,其边界具有非整数维度的特点。

Minkowski型分形是一种特殊的分形,其边界具有连续的曲率。

在物理学中,波动方程是描述波动现象的重要方程之一。

本文将探讨具有Minkowski型分形边界的二维波动方程。

首先,我们来回顾一下二维波动方程的一般形式:∂²u/∂t² = c²(∂²u/∂x² + ∂²u/∂y²)其中,u是波函数,t是时间,x和y是空间坐标,c是波速。

现在考虑一个具有Minkowski型分形边界的二维区域。

我们可以将该区域划分为许多小区域,每个小区域都具有与整体相似的形状。

假设每个小区域的边界长度为L,那么整个区域的边界长度可以表示为N×L,其中N是小区域的数量。

在这个具有Minkowski型分形边界的区域中,波函数u满足以下边界条件:u(x,y) = 0, 当(x,y)在边界上时这意味着波函数在边界上必须为零,即波不能穿过边界。

为了求解这个具有Minkowski型分形边界的二维波动方程,我们可以采用分离变量法。

假设波函数可以表示为u(x,y,t) = X(x)Y(y)T(t),将其代入波动方程中,可以得到以下三个方程:∂²X/∂x² + k²X = 0∂²Y/∂y² + k²Y = 0∂²T/∂t² + c²k²T = 0其中,k是一个常数。

解这三个方程可以得到波函数的形式。

由于边界条件的限制,我们需要选择满足边界条件的解。

在具有Minkowski型分形边界的情况下,边界的形状是连续变化的,因此我们需要考虑边界的分形特性。

具体来说,我们可以使用分形几何的方法来描述边界的形状。

通过计算边界的分形维度,我们可以得到边界的特征参数。

然后,我们可以将这些参数代入波函数的解中,得到满足边界条件的波函数。

二维定常不可压缩N-S方程无量纲分析

二维定常不可压缩N-S方程无量纲分析

二维定常不可压缩N-S 方程无量纲分析一、引言计算流体力学的控制方程通常认为是N-S(Navier-Strokes)方程组,包含了能量方程、动量方程、连续性方程等方程组的总称。

当考虑流体的黏性时,作用在流体质点上的力除了质量力、法向应力(垂直于作用面的压力)外,还有与作用面相切的切向力,N-S 方程建立了流体微团的动量变化率与作用在微团上的惯性力,压力以及粘性剪切力之间的关系,反映了黏性流体运动的基本规律,对计算流体力学有着十分重要的意义。

本文旨在对二维定常不可压缩N-S 方程进行无量纲化,方便简化计算和分析相似实验。

量纲分析就是对有量纲的物理方程进行参数的组合,实现参数和方程的无量纲化,将方程无量纲化有以下几点好处:(1)方程形式可以得到简化并且可能减少方程个数,进而提高实际计算速度;(2)通过无量纲化尽可能的减少方程中的常数运算,将这些常数转化为某个特征参数,这样可以降低计算难度;(3)防止方程中的物理参数在数量级上造成差异,从而降低精度损失;(4)将方程中的物理量无量纲化后容易实现计算中的相似模拟。

流体力学中的相似通常可以分为几何相似、运动相似和动力相似。

流动相似的概念来源于几何相似的概念,两个流动如果相似,例如模型流动与实际流动相似,则其流场中相应点上各同类物理量将具有各自固定的比例关系,也即可将模型实验的成果应用于实际流动中。

相似原理指出,两个流动若相似必满足一定条件,即满足几何相似、运动相似、动力相似,这些条件还应包括边界条件和初始条件相似。

根据相似原理,两个流动现象只要同时满足上面的相似条件,它们之间就存在相似关系,其对应物理量都成一定的比例关系。

在应用中,首先需要分析所要研究的流体,找出影响流动问题的作用力,我们只需要满足一个主要作用力相似,而不必计较其它作用力是否达到相似。

例如对于一些流动现象,只要流动的雷诺数不是很大,一般其相似条件都依赖于雷诺数。

雷诺数是用来判断流体流动特性的无量量,共有18个应力分量X 轴的运动微分方程:(2.1)最后导出沿x 轴的(2.2) (2.3)(2.4)纲量,对于封闭环境内的流动, 当雷诺数小于 2300时的流动为层流, 能用N-S方程表示;当雷诺数大于4000时的流动为湍流,不能用 N-S 方程表示。

二维不可压缩定常流动(平板)边界层方程

二维不可压缩定常流动(平板)边界层方程

二维不可压缩定常流动(平板)边界层方程下载温馨提示:该文档是我店铺精心编制而成,希望大家下载以后,能够帮助大家解决实际的问题。

文档下载后可定制随意修改,请根据实际需要进行相应的调整和使用,谢谢!并且,本店铺为大家提供各种各样类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,如想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by the editor. I hope that after you download them, they can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you!In addition, our shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!二维不可压缩定常流动(平板)边界层方程一、引言在流体力学领域,边界层方程是研究流体在边界附近流动的基本方程之一。

格子Boltzmann方法模拟二维轴对称狭窄血管内的脉动流

格子Boltzmann方法模拟二维轴对称狭窄血管内的脉动流

201020446(2) 北京师范大学学报(自然科学版)Journal of Beijing Normal University (Natural Science ) 139 格子Boltzmann 方法模拟二维轴对称狭窄血管内的脉动流3张立换 康秀英 吉驭嫔(北京师范大学物理学系,100875,北京)摘要 将格子Boltzmann 方法应用到二维轴对称余弦狭窄血管模型,模拟比较加入脉动后流场速度、压强和剪切应力分布,并详细分析了不同狭窄模型、Reynolds 数和Womersley 数对血液流动规律的影响,从而为研究血管壁病变和动脉硬化形成机制提供了有用的理论参考.关键词 格子Boltzmann 方法;Reynolds 数;Womersley 数;脉动流;动脉狭窄3北京师范大学青年科学基金资助项目通信作者收稿日期:2009205219 格子Boltzmann 方法(lattice Boltzmann met hod ,简称LBM )是20世纪80年代迅速发展起来的一种新的流体动力学数值模拟方法[122].与以宏观连续方程的离散化为基础的传统数值方法不同,LBM 从微观层次出发,采用统计物理方法得出流体的宏观特性,而且在可操作性方面,它计算方便,编程易于实现,边界易于处理等优点已经得到广泛地证实.由于心血管疾病多集中于具有复杂几何形状和具有复杂流动特性的区域,流动区域和剪切应力的分布对理解、诊断和治疗这种疾病有很重要的作用.近年来,LBM 在血液动力学方面的应用越来越受到重视[326].本文的主要工作是用格子Boltzmann 方法模拟二维轴对称狭窄血管内脉动流的流动特性.首先对狭窄血管内定常流特性进行了研究,模拟比较不同狭窄模型和不同Reynolds 数对管壁切应力、压强和压力梯度分布的影响.然后对二维轴对称狭窄血管内脉动流的流动特性进行了研究,模拟比较在改变Reynolds 数、Womersley 数时动脉血流的流动特性,找到动脉血流的非定常性对狭窄血管中流场速度、压强和剪切应力分布的影响,从而对常见的心血管疾病发展机制给出物理解释,为进一步分析动脉粥样硬化的形成、发展及其影响提供新的研究方法和理论参考.1 二维轴对称狭窄血管内定常流特性的研究111 管壁几何模型 假定血管的狭窄处为轴对称,如图1所示,狭窄形状采用常用的余弦形状,即y =h2[1+co sπL(x -x 0)],(1)图1 二维轴对称余弦狭窄模型 其中h 是狭窄的最大高度,对应于x =x 0处,L 是狭窄总长度的一半,L x 是血管段的长度,L y 是狭窄发生前的血管宽度.112 数值计算 模拟中,计算网格选为N x ×N y =300×40,狭窄中心处为x 0=121,通过调整h 和L 来控制血管狭窄程度.血管出入口采用压强边界条件[7],管壁边界采用Mei 改进的曲线边界条件[8].为了研究不同狭窄情况下管壁的切应力、压强和压强梯度的变化规律,我们选择3个不同的狭窄模型,如表1.表1 不同的狭窄模型狭窄模型M1M2M3狭窄高度h L y /8L y /4L y /4狭窄长度2L16h 8h 16h 在保证Reynolds 数(Re =ρUL y μ=UL yν,ν=μ/ρ为流体运动学黏滞系数,U 为入口附近的平均速度)一定时,计算得3种模型管壁切应力、压强和压强梯度见 140 北京师范大学学报(自然科学版)第46卷 图2~4.Re =114,狭窄中心x 0=121.图2 3种狭窄模型下管壁切应力分布 从图2中可以看出,管壁切应力振荡的负峰值在靠近狭窄中心(x 0=121)的上游,这个峰值达到一定值后,该部位血管内皮组织易发生机械应力损伤.当狭窄长度一定时,狭窄高度越大,切应力的负峰值越大,如图2中的M1和M2;当狭窄高度一定时,狭窄长度越短,切应力的负峰值越大,如图2中的M2和M3.同时也可以看出在狭窄处的下游切应力变小,特别是M2,血液容易在此处发生流体分离.模拟得到狭窄区域的压强和压强梯度分布如图3和4所示.在相同狭窄长度下,狭窄高度越大,血管狭窄上游压强下降越大,下游压强上升越大,同时狭窄区域前后的压强落差越大,如图3中的M1和M2.另一方面,在相同狭窄高度下,狭窄长度越长,血管狭窄上游压强下降越大,同时狭窄区域前后的压强落差越大,如图3中的M2和M3.压强梯度在狭窄区域波动加图3 管壁上压强分布(Re =114),p 0是狭窄发生前的压强,u 0是x =20处的中心流速 图4 管壁上的压强梯度分布(Re =114) 剧,压强梯度波动最大的是狭窄模型M2(图4),其对应的切应力负峰值也为最大值,狭窄部位管壁切应力与压强梯度的变化规律具有相似性.选择模型M2,比较管壁切应力和狭窄附近的流场分布随Re 的变化规律,如图5和6.从图5中可以看出,狭窄模型一定时,随着Re 的增加,管壁切应力增大,在狭窄区域的下游,切应力的增加相对减小,这是由于出现了流体分离,如图6的流场分布.图6显示了模型M2在不同Re 下狭窄附近的流场分布,可以看出,随着Re 的增大,在狭窄下游管壁处出现流动分离区,且Re 越大,流动分离区越大.113 分析与结论 通过改变参数,我们获得了大量有关狭窄血管中的流场的信息.模拟结果表明,血管局部图5 管壁切应力随Reynolds 数的变化曲线(狭窄模型M2) 第2期张立换等:格子Boltzmann方法模拟二维轴对称狭窄血管内的脉动流141图6 不同Re下的流场分布(M2,Re=114、215、318)狭窄会对血液的流动状态产生明显的影响,从而带来一系列的生理和病理方面的复杂变化.例如,动脉硬化斑块主要发生在几何形状急剧变化和高Re流动状态的血管内.在动脉硬化斑块发展的初期,血管狭窄度比较小,对于黏度是常数的血液流体,其Re比较小,无流动分离,管壁切应力可能达到临界应力值,对狭窄上游血管壁内皮细胞造成损伤,使壁面进一步异常增生,导致血管狭窄度增加,进而导致此处流动Re的增加.当血管狭窄增大到一定值时,在狭窄下游管壁附近就会有流动分离区形成,在该区域内血液会发生滞留,血液中的血小板和纤维蛋白就会沉积,并在血管壁处形成网络结构致使血液中的脂质颗粒沉积,而最终导致动脉粥样硬化现象的出现.同时,狭窄度较大时,对应的压力梯度的值也会较大,也可以反映病变血管的异常血液流动情况.2 二维轴对称狭窄血管内脉动流的流动特性选择模型M2为研究对象,模拟中选取周期T=10000,流动的Womersley数(α=L y2ων,ω=2πf=2π/T是脉动的角频率)为α=31357,入口压强随时间周期性变化,即p(0,t)=Δp cosωt+p out,Δp为一常量,出口压强pout设为定值,图7显示一个周期8个不同时刻的脉动流管道中心中轴线上的压强分布.从图7中可以看出,中轴线上的压强不是线性变化,在靠近狭窄部位压强下降幅度明显增加,在最大狭窄处附近压强出现极小值,狭窄下游压强又逐渐回升,远离狭窄后,压强变化逐渐恢复类直管变化趋势,并且压强随时间的波动存在一定的滞后,如图中1/8T和7/8T,2/8T和6/8T以及3/8T和5/8T不完全重合.狭窄中心x0=121,狭窄长度为78.图7 iT/8时刻中轴线上的压强分布 142 北京师范大学学报(自然科学版)第46卷 脉动流前半周期的流场分布如图8所示.从图中可以看出,在T/4时刻,在狭窄下游管壁附近开始出现流动分离区,且分离区逐渐扩大,如3T/8时刻,接着又缓慢消失,如T/2时刻,流体平滑地流过凸包.图8 脉动流在前半周期内不同时刻的流场分布 需要注意的是心脏的周期性泵血作用使动脉中的血液以脉动的形式流动,动脉中血液流动的参量———压强、流量等流动参数也会随时间变化,虽然动脉中血液的流动是脉动流而不是定常流,但动脉中血流的方向平均来说却是始终不变的,即总是从动脉流向毛细血管,再流向静脉.因此,可以把由心脏收缩和舒张所引起的动脉中的脉动流看作是一定常流分量与一振荡分量的叠加,即在图8所示的流场分布中叠加上一个定常流,最终倒流的出现时间将非常短暂,且流速很小.对应于一个周期中的不同时刻,我们发现,管壁切应力的随时间的波动也存在一定的滞后.如图9给出前半周期的切应力分布.3 结束语我们讨论了二维余弦狭窄血管中血液流动的切应力、流场速度、压强和压强梯度在不同狭窄模型和不同图9 前半周期内管壁切应力的变化曲线Re下的分布规律,所得结论与用其他实验,理论和数值模拟得到的结论相同[9211],但用LBM方法编程简单,参数易于选择,从分布函数就可以得到所有主要宏 第2期张立换等:格子Boltzmann方法模拟二维轴对称狭窄血管内的脉动流143观量,证实了LBM在此模型下的适用性.考虑到血液流动的脉动性,研究了一个脉动周期中流场的变化特点,并与定常流动比较,分析其差异.由于Womersley数的选择在血流参数范围内,故认为上述结论具有参考性.值得注意的是,流动分离区并不同于定常流动所述那样在管壁处停留,而是随着时间的演化,流动分离区间歇性的出现,如对α=710797的流场分布模拟显示,与α=31357的不同点是流动分离区在管壁附近产生后,随着时间的推移,又会向管轴附近发展.与定常流情况下在Re达到300后才出现明显的分离区不同,对于脉动流,在Re较小时,就已经可以观察到明显的流动分离区了.4 参考文献[1] Qian Y H,d’Humieres D,lallemand ttice B GKmodels for Navier2Stokes equation[J].Europhys Lett, 1992,17:479[2] Chen H,Chen S,Matthaeus W H.Recovery of theNavier2Stokes using a lattice2gas Boltzmann method[J].Phys Rev A,1992,45:R5339[3] Artoli A M,Kandhai D,Hoef sloot H ttice B GKsimulations of flow in a symmetric bif urcation[J].FutureG eneration Computer Systems,2004,20:909[4] Boyd J,Buick J,Cosgrove J A,et al.Application of thelattice Boltzmann model to simulated stenosis growth in a two2dimensional carotid artery[J].Phys Med Biol,2005, 50:4783[5] Li H B,Fang H P,Lin Z ttice Boltzmannsimulation on particle suspensions in a two2dimensional symmetric stenotic artery[J].Phys Rev E,2004,69: 031919[6] 康秀英,刘大禾,周静,等.用格子Boltzmann方法模拟动脉分叉流场[J].北京师范大学学报:自然科学版,2005, 41(4):364[7] Z ou Q,He X.On pressure amd velocity boundaryconditions for the lattice Boltzmann B GK model[J].Phys Fluids,1997,9(6):1591[8] Mei R,L uo L S,L uo Shyy W.An accurate curvedboundary treatment in the lattice Boltzmann method[J].J Comput Phys,1999,155:307[9] 姚力,李大治.刚性轴对称狭窄血管内压强及其梯度的研究[J].应用数学和力学,2006,27(3):311[10] 刘国涛,王先菊,艾保全,等.狭窄动脉血管中Poiseuille流动对管壁切应力的影响[J].中山大学学报:自然科学版,2004,4(6):29[11] 秦杰,刘辉,孙利众,等.刚性狭窄管内血流压力分布的研究[J].生物力学,1989,4(6);57SIMU LATING B LOOD FLOW IN A TWO2DIMENSIONALSYMMETRIC STENOTIC ARTER Y BYTHE LATTICE BOL TZMANN METH ODZHAN G Lihuan KAN G Xiuying J I Yupin(Depart ment of Physics,Beijing Normal University,100875,Beijing,China)Abstract In t his st udy t he lattice Boltzmann met hod has been applied to a two2dimensional symmet ric stenotic artery.The velocity,p ressure and shear st ress distribution of blood flow were simulated and compared when p ulsatio n over t he blood was added.We have observed t he impact of blood flow when changing t he steno sis struct ure,Reynolds number and Womersley number.These data provide a p hysical explanation for blood vessel lesions and arterio sclero sis.K ey w ords lattice Boltzmann met hod;Reynolds number;Womersley number;p ulsating blood;steno sed artery。

不可压缩navier-stokes方程

不可压缩navier-stokes方程

以下是不可压缩navier-stokes方程的具体描述和数学表达方式:
不可压缩Navier-Stokes方程是描述不可压缩流体运动的方程。

它是由法国数学家Navier和Stokes在19世纪初期研究流体运动时提出的。

不可压缩Navier-Stokes方程包含了流体运动的连续性方程和动量方程。

连续性方程描述了流体的质量守恒,即流体在任意时刻体积不变。

动量方程则描述了流体的动量守恒,即流体的加速度与施加于它的力成正比。

不可压缩Navier-Stokes方程的数学表达式如下:
连续性方程:
$$\nabla \cdot \boldsymbol{v} = 0$$
动量方程:
$$\rho \frac{\partial \boldsymbol{v}}{\partial t} + \rho (\boldsymbol{v} \cdot \nabla)\boldsymbol{v} = -\nabla p + \mu \nabla^2 \boldsymbol{v} + \boldsymbol{f}$$
其中,$\boldsymbol{v}$是流体速度矢量,$\rho$是流体密度,$p$是压力,$\mu$是粘度系数,$\boldsymbol{f}$是外力源矢量。

不可压缩Navier-Stokes方程的求解非常困难,因为它包含了非线性项和高阶微分方程。

目前,只有一些特殊情况下的解析解可用,而大多数情况下需要使用数值方法进行求解。

二维定常不可压缩N-S方程无量纲分析

二维定常不可压缩N-S方程无量纲分析

.二维定常不可压缩N-S方程无量纲分析一、引言计算流体力学的控制方程通常认为是N-S(Navier-Strokes)方程组,包含了能量方程、动量方程、连续性方程等方程组的总称。

当考虑流体的黏性时,作用在流体质点上的力除了质量力、法向应力(垂直于作用面的压力)外,还有与作用面相切的切向力,N-S方程建立了流体微团的动量变化率与作用在微团上的惯性力,压力以及粘性剪切力之间的关系,反映了黏性流体运动的基本规律,对计算流体力学有着十分重要的意义。

本文旨在对二维定常不可压缩N-S方程进行无量纲化,方便简化计算和分析相似实验。

量纲分析就是对有量纲的物理方程进行参数的组合,实现参数和方程的无量纲化,将方程无量纲化有以下几点好处:(1)方程形式可以得到简化并且可能减少方程个数,进而提高实际计算速度;(2)通过无量纲化尽可能的减少方程中的常数运算,将这些常数转化为某个特征参数,这样可以降低计算难度;(3)防止方程中的物理参数在数量级上造成差异,从而降低精度损失;(4)将方程中的物理量无量纲化后容易实现计算中的相似模拟。

流体力学中的相似通常可以分为几何相似、运动相似和动力相似。

流动相似的概念来源于几何相似的概念,两个流动如果相似,例如模型流动与实际流动相似,则其流场中相应点上各同类物理量将具有各自固定的比例关系,也即可将模型实验的成果应用于实际流动中。

相似原理指出,两个流动若相似必满足一定条件,即满足几何相似、运动相似、动力相似,这些条件还应包括边界条件和初始条件相似。

根据相似原理,两个流动现象只要同时满足上面的相似条件,它们之间就存在相似关系,其对应物理量都成一定的比例关系。

在应用中,首先需要分析所要研究的流体,找出影响流动问题的作用力,我们只需要满足一个主要作用力相似,而不必计较其它作用力是否达到相似。

例如对于一些流动现象,只要流动的雷诺数不是很大,一般其相似条件都依赖于雷诺数。

雷诺数是用来判断流体流动特性的无量1 / 6.纲量,对于封闭环境内的流动,当雷诺数小于2300时的流动为层流,能用N-S 方程表示;当雷诺数大于4000时的流动为湍流,不能用N-S方程表示。

7不可压缩流体动力学基础解析

7不可压缩流体动力学基础解析
v x x v x y v x z x 2 y 2 z 2
图7-2 流体微团运动速度分量
如图 7-2 所示,在流场中任取一平行六面体的流体微 团,以该流体微团的运动速度为讨论对象。已知 t 瞬时中 心点 O(x, y, z)的速度 v v x i v y j v z k 。在该流体
具体分析如下:
v y x v y y vy x 2 y 2
v y x v y y vy x 2 y 2
v x v x y vx x x 2 y 2
v y x v y y vy x 2 y 2
D y δy C x
A
vx
vr 1 v 1 v 2vr v cot 0 r r r sin r r
式中 r 为径矩;为纬度; 为径度。
【例】
已知不可压缩流体运动速度
v 在 x ,y
两个轴方向的分量
为 vx 2x 2 y , v y 2 y 2 z 。且在 vz 。 方向的速度分量 【解】对不可压缩流体连续性方程为:
首先推导在笛卡儿坐标系中微分形式的连续性方程。
如图7-1
微元六面体
设该微元六面体中心点O(x, y, z)上流体质点的速度
密度为 ,于是和 在
x 轴垂直的两个平面上的质量流量如图所示。
dx vx vx dydzdt x 2
v x v y vz
方向上, dt 时间通过EFGH面流入的流体质量为: x ( a)
它适用于理想流体和粘性流体、定常流动和非定常流动。 在定常流动中,由于
0 t
v x v y v z 0 x y z
对于不可压缩流体(

8第八讲 不可压缩的Navier-Stokes方程的解法

8第八讲 不可压缩的Navier-Stokes方程的解法

不可压缩Navier-Stokes方程的数值方法不可压缩Navier-Stokes方程的特点人工压缩性方法(求解定常方程)投影法涡量-流函数方法(二维问题)SIMPLE方法的连续性方程,得到2uu u ±=±迎风差分,建议采用高阶的)''(,,1j i j i p p −+α带入离散的连续性方程:/)(/)(12/1,12/1,1,2/11,2/1=Δ−+Δ−+−+++−++y v v x u u n j i n j i n j i n j i )''(,,1*,2/11,2/1j i j i j i n j i p p u u −+=++++α)''(,1,*2/1,12/1,j i j i j i n j i p p v v −+=−+++β得到离散的压力Poisson 方程:)5(ˆ'''''1,41,3,12,11,cp c p c p c p c p j i j i j i j i j i ++++=−+−+求解后,得到压力修正值:ji p ,'(4)带入(4)时得到n+1时刻的速度具体步骤:1)已知n 时刻的速度、压力2)预估压力(可取为n 时刻的压力)3)带入(1)(2)式,解出(隐格式,需迭代求解)4)求解压力的修正方程(5)得到修正压力5)带入(4)式,得到n+1时刻的速度及压力6)推进求解直到给定时刻(或收敛)*p **,v u 如该步改用显格式,则为(离散型)投影法'*1p p p n +=+楼群的三维视图来流的速度分布假定速度在10m 高处的大小为,其余高度的分布则采用风廓线分布:其中Z 是距离地面的高度,是地面粗糙系数,我们在模拟的过程中取它为0.28.10U 10(/10)U U Z α=α算例1. 北风s mU/310=5m高的速度分布15m高的速度分布35m高的速度分布70m高的速度分布80m高的速度分布5m高的流线示意图南北向截面流线示意图近壁面压力分布汽车的表面三角网格模拟来流:正前方u=20m/s(72km/h)汽车基本参数:总宽:900mm总长:2600mm总高:600mm车轮半径:180mm车轮厚度:120mm顶层长度, 宽度:800mm, 620mm流速分布图:低于21m/s的截面图流速分布图:y=0m二维流线:y=0m三维速度面,u=14m/s压力分布图。

第五章-不可压缩流动的数值方法_73005466

第五章-不可压缩流动的数值方法_73005466

的动量方程,但由于数值积分的误差,其结果将由于选择不同的积分路径,而得 到不同的积分值。即所谓压力漂移现象,所以求解压力场 P,通常并不推荐积分
方法,而是导出压力场满足的微分方程,进行压力场的数值求解。
⎧ ∂u 1 ∂p 1 2 ∂u ∂u ⎪ ∂t + u ∂x + υ ∂y = − ρ ∂x + Re ∇ u 方程: ⎪ 1 ∂p 1 2 ∂υ ∂υ ⎪ ∂υ +u +υ =− + ∇υ ⎨ Re t x y y ρ ∂ ∂ ∂ ∂ ⎪ ⎪ ∂u ∂υ =0 ⎪ + ⎩ ∂ x ∂y
ς=
∂u ∂υ − ∂y ∂x
这是一个对流-扩散(输运)方程加上 Poisson 方程的偏微分方程组问题,其差分 数值方法包括以下几部分内容: (一) 涡输运方程的求解(略) (二) 流函数的 Poisson 方程的求解(略) (三) 在求解涡运输方程时的涡量边界条件的确定问题 (四) 在求出了 ς ,ψ ,后,如何确定压力场 P。 此处讨论后两个问题,即: * 涡量的边界条件提法问题:固壁边界的涡量条件提法-壁涡公式 由粘性条件:壁面上:u=v=0
( ρu ) e
Φ + ΦW Φ − ΦW ΦE + ΦP Φ − ΦP − ( ρu ) w P = [Γe E − Γw P ] 2 2 (δx ) e (δx ) w
整理:
ΦP[
Γ Γ Γ Γ ( ρ u )e ( ρu) w ( ρu )e ( ρu ) w + e − + w ] = ΦE[ e − ] + ΦW [ w + ] 2 (δx )e 2 (δx ) w (δs )e 2 (δx ) w 2
{

格子玻尔兹曼方程

格子玻尔兹曼方程

格子玻尔兹曼方程格子玻尔兹曼方程(GridvonNeumannequation)也称为离散热力学方程,是二维或三维的热动力学模型,由John von Neumann在1950年提出,并且被用来模拟计算机和集群系统的并行性能。

它的基本原理是描述数学模型中流体的凉却过程,也就是总体热量经过散布和消耗而转变为低温和低能量。

根据格子玻尔兹曼方程,热量在不同空间不同方向传播,当物体热量达到某种最终温度时,可以得出一些实际结论和应用。

格子玻尔兹曼方程有两种形式:N维和N+1维。

N维表示热量穿越网格的每一个方向,而N+1维表示热量穿越整个网格体系。

格子玻尔兹曼方程的另一个有趣特性是,经济上改变温度,其最终状态仍可能让热量更均匀地分布,以及有一定的实用价值。

因此,基于格子玻尔兹曼方程的研究可改变计算机系统中热量分布,以提高系统性能。

格子玻尔兹曼方程是物理机器学习中一个重要的方法。

它可以用来学习物理过程,在仿真模拟的时候,它可以用来代替算法中的模型,以表示物理过程。

此外,它还可以预测某种物理现象,如内涡流、热传输等,并且可以用来评估经济和政治行为的结果。

格子玻尔兹曼方程可以用于量化复杂非线性系统的热力学行为,比如气体的复杂流动。

基于它,我们可以估计某一种物质在不同环境条件下多维空间中的热量变化和流动,从而提供优化设计方案以提高效率或者降低能源消耗。

格子玻尔兹曼方程也可以应用在决策环境中,比如模拟投资和金融决策,为投资者提供最佳结果。

这种模拟可以基于金融市场的历史数据来进行,以帮助投资者获得最佳投资组合,并能够有效预测未来投资变化的结果。

格子玻尔兹曼方程的另一个应用领域是仿生工程。

仿生工程是一种机器人技术,将自然界的动物模拟出来,从而实现更有效的机器人。

根据格子玻尔兹曼方程,机器人可以根据温度分布,来实现依据周围环境温度的控制,从而实现仿生学研究。

格子玻尔兹曼方程也可以应用于计算机硬件的设计,如数据中心的温度控制,这可以更好地控制数据中心的内部温度,以避免内部温度过高而影响系统性能的情况。

不可压缩流体二维边界层概述课件讲解

不可压缩流体二维边界层概述课件讲解

流体力学性质与边界层传输特性的关系
1 黏性
2 流量
影响了流体的运动状态和 流动结构。边界层的形成 和发展也与黏性密切相关。
决定着流体的运动量,它 也是边界层流动的核心参 数之一。
3 温度
热传递过程对流体的运动 有很大的影响,决定着边 界层内流体的温度分布。
非定常边界层流动
非定常边界层流动指的是流体速度和边界的运动状态随时间而变化的流动, 研究它所涉及到的问题比较复杂,但对于某些实际应用非常重要。
旋转
旋转会引起流体的湍流,对于 边界层流动也有着重要的影响。
边界层流动在工业生产中的应用
• 与工业生产过程中的传热、传质等过程密切相关,因此被广泛应用。 • 可用于优化能量利用和资源消耗,提高生产效率,降低生产成本。 • 在航空、航天等领域中的应用也非常重要。
相关流动情况,观 察流体的运动规律。
流动可视化技术
通过颜料、激光等工具追踪流体 的运动,观察其流速、涡旋等特 征。
传热实验
通过分析流体传热规律和温度分 布情况,研究边界层流动中的传 热特性。
二维边界层的无量纲分析
无量纲分析是在流体力学中广泛应用的一种方法,可以将复杂的流动过程简化为一组基本的无量纲关系。
边界层压力梯度方程
边界层压力梯度方程是研究边界层内流体压力变化的方程。它主要涉及压力 的传递和液体的运动状态。
二维边界层的概念
二维边界层是指在二维平面内的边界层流动。它有着更加简洁优美的数学模型,更容易理解和求解。
二维边界层的流场结构
1
加速段
2
流速迅速增加,流线弯曲。
3
减速段
4
流率减慢,流线逐渐变直。
不可压缩流体二维边界层 概述课件讲解

不可压缩流体的数学和物理解析

不可压缩流体的数学和物理解析

不可压缩流体的数学和物理解析
不可压缩流体是指在流动过程中密度基本保持不变的流体。

从数学和物理角度来解析不可压缩流体,我们可以从以下几个方面来进行讨论:
1. 数学解析:
在数学上,不可压缩流体通常可以通过连续方程来描述。

连续方程是描述流体运动的基本方程之一,它表示了流体质点的密度随时间和位置的变化关系。

对于不可压缩流体,连续方程可以简化为一个常数,即密度不随时间和位置的变化而变化。

这样的数学特性使得不可压缩流体的数学分析相对简化,可以更容易地进行数值模拟和解析求解。

2. 物理解析:
从物理角度来看,不可压缩流体的性质使得其在流动过程中密度变化很小,因此可以近似为密度恒定。

这样的性质在实际工程和科学问题中有着重要的应用,例如在空气动力学、水力学和石油工程等领域。

物理上,不可压缩流体的性质使得其流动更为稳定,
且在一些情况下可以简化流体力学方程的求解和分析过程。

3. 数值模拟:
对于不可压缩流体的数学和物理解析,数值模拟也是一种重
要的方法。

通过数值模拟,可以对不可压缩流体的流动进行模拟和
分析,得到流场的各种参数和特性。

数值模拟在工程实践中有着广
泛的应用,可以帮助工程师和科学家更好地理解不可压缩流体的行
为和特性。

综上所述,不可压缩流体的数学和物理解析涉及到连续方程、
物理特性以及数值模拟等多个方面。

通过综合运用数学和物理知识,可以更好地理解和分析不可压缩流体的性质和行为。

求解二维对流扩散方程格子Boltzmann方法

求解二维对流扩散方程格子Boltzmann方法

1引言格子Boltzmann方法(简称LBM)是20世纪80年代末提出来的一种流体力学计算方法[1-3],与传统的宏观连续方法和微观的分子动力学方法不同,格子Boltzmann 方法是一种介于宏观和微观之间的介观方法。

与传统的数值模拟方法相比,格子Boltzmann方法物理背景清晰,处理简单,具有天生的并行特性,而且边界容易处理,程序易于实施,从而得到国内外众多学者的关注。

近二十余年来,格子Boltzmann方法已经成功运用到流体运动的数值模拟[4-6],在许多传统模拟方法难以胜任的领域,如微尺度流动、多孔介质、生物流体、磁流体、晶体生长等,格子Boltzmann方法都可以进行有效的模求解二维对流扩散方程的格子Boltzmann方法彭碧涛1,郑洲顺1,2,刘红娟1,汤慧萍2,王建忠2PENG Bitao1,ZHENG Zhoushun1,2,LIU Hongjuan1,TANG Huiping2,WANG Jianzhong21.中南大学数学与统计学院,长沙4100832.金属多孔材料国家重点实验室,西安7100161.School of Mathematics and Statistics,Central South University,Changsha410083,China2.State Key Laboratory of Porous Metal Materials,Xi’an710016,ChinaPENG Bitao,ZHENG Zhoushun,LIU Hongjuan,et ttice Boltzmann method for solving two dimensional convection-diffusion puter Engineering and Applications,2015,51(23):68-73.Abstract:In view of two-dimensional convection diffusion equation,it deduces the conditions of equilibrium distribution function must to be satisfied in every velocity directions based on D2Q4Lattice velocity and gives the specific expres-sions of the equilibrium distribution function.Through Chapman-Enskog multi-scale analysis technology,taking the time scale directly to the second order and the spatial scale to the first order,the equilibrium distribution function can recover the original convection-diffusion equation,thus the new D2Q4Lattice Boltzmann(LB)model for solving two-dimensional convection diffusion equation is ing the model,it implements a diffusion equation and two convection dif-fusion equation with different initial and boundary conditions,the numerical results are in good agreement with analyticsolutions,furthermore the boundary error is very low compared with related document,therefore the effectiveness of the new model is verified.Key words:convection-diffusion equation;lattice Boltzmann method;D2Q4model;multi-scale analysis technology摘要:针对二维对流扩散方程,基于D2Q4格子速度,用Chapman-Enskog多尺度分析技术,将时间尺度取为二阶,空间尺度取为一阶,推导了各个速度方向上的平衡态分布函数所满足的条件,给出了简单且对称的平衡态分布函数表达式,所得到的平衡态分布函数能正确地恢复出二维对流扩散方程,从而构建了一种新的求解二维对流扩散方程的D2Q4格子Boltzmann(LB)模型。

基于FTM方法的二维Kelvin-Helmholtz不稳定性数值模拟

基于FTM方法的二维Kelvin-Helmholtz不稳定性数值模拟

基于FTM方法的二维Kelvin-Helmholtz不稳定性数值模拟尚文强;张莹;陈昭奇;袁志平;董博恒【摘要】使用界面跟踪法FTM(Front Tracking Method)对二维不混溶、不可压缩流体的K-H (Kelvin-Helmholtz)不稳定性进行数值模拟.研究表明,速度梯度层越厚,界面在水平分量中移动越快,卷起越少;初始水平速度差越大,界面卷起越多,内扰动增长速度越快,K-H不稳定性的特征形式更加明显;此外,在Neumann边界条件(即无滑移边界条件)下界面的扰动发展得比Dirichlet边界条件(即对称边界条件)下的扰动快.由于Dirichlet边界中的边界层,在开始时刻涡量扩展到两侧,影响了K-H 不稳定性的生长速率;而在Neumann边界条件下涡量由于初始水平速度差,在界面中心聚集.最后,研究了不同边界条件下各种理查德森数对K-H不稳定性的影响.【期刊名称】《计算力学学报》【年(卷),期】2018(035)004【总页数】7页(P424-430)【关键词】Kelvin-Helmholtz不稳定性;速度梯度层厚度;Neumann边界条件;Dirichlet边界条件;界面跟踪法【作者】尚文强;张莹;陈昭奇;袁志平;董博恒【作者单位】南昌大学机电工程学院,南昌 330031;南昌大学机电工程学院,南昌330031;南昌大学机电工程学院,南昌 330031;南昌大学机电工程学院,南昌330031;南昌大学机电工程学院,南昌 330031【正文语种】中文【中图分类】TP273;O242.11 引言K-H(Kelvin-Helmholtz instability)[1]不稳定性是自由剪切流的无粘不稳定性。

一般认为,不互溶的混合流体中,上层流体与下层流体间的速度差异足够大时,界面上存在的微小扰动会使界面变形,出现Kelvin-Helmholtz不稳定性。

K-H不稳定性的数值研究一直是计算流体力学的经典问题。

div2

div2

计算物理,3卷2期,1986, 217-226解不可压缩流体力学问题的降阶法 Ⅱ. 空间vh的基函数和误差估计于 欣(中国科学院力学研究所) 1985年9月6日收到摘 要本文对于一大类数值求解二维Navier-Stokes 方程边值问题的有限元格式给出了零散度空间V h 的一组简单基函数,讨论了速度数值误差对压力的数值解的影响,并提出一个改进算法。

一、前 言设2R ⊂Ω是k 连通的有界区域。

考虑Navier-Stokes 方程边值问题:+∇⋅+∆-u u u )(μgrad f p =, 中在Ω, )1.1(1………………………..div u =0, 中在Ω, )1.1(2………………………… 0=u , 上在Ω∂,)1.1(3其中u=(21,u u )是速度向量,p 是压力,常数μ>0。

设T h 是Ω的一个可允许的有限剖分。

hT k h ∈=Ω K 。

速度的近似空间为h U 。

压力的近似空间为hP (参见第二章中的定义)。

解边值问题(1.1)的混合有限元问题为[3][4]求(hu ,hp )hhP U ⨯∈,使得)2.1(0h h h h h h h h h h U v v f p v b v u A ∈∀〉〈=+,,),(),(, )2.1(1 0),(=h h h q u b ,,h h p q ∈∀ )2.1(2其中),,(),(),(10h h h h hh hh h h v u u a v u a v u A +=,,,dx v f v f h h h h ⋅⎰=〉〈Ω〈=μ),(0h h h v u a grad h u ,grad h h v 〉 .h h i h i h i h i i h hhhh gradv u gradu v w v u w a 〉-〈=∑=,21),,(211h h hh hhdiVv v b q q 〉-〈=,),(在实际数值求解中,有限元问题(1.2)的计算量是非常大的。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

定义 2.2:设 (t, x,v) 是任意一点,且 ( x,v) ∉γ 0 ,令 (t0, x0,v0 ) = (t, x,v ) ,我们有如下定义:
( ) ( ) tk +1, xk +1,vk +1= tk − tb ( xk ,vk ), xb ( xk ,vk ),vk +1 ,
{ } ∑ X= cl (s;t, x,v ) ( ) ( ) 1[tk+1,tk ) s xk + s − tk vk , k
Pure Mathematics 理论数学, 2020, 10(2), 128-138 Published Online February 2020 in Hans. /journal/pm https:///10.12677/pm.2020.102019
Incompressible Limit of the Two Dimensional Boltzmann Equation
Tingting Gao
School of Mathematics, South China University of Technology, Guangzhou Guangdong
本文主要结果如下:
定理 1.1:设 Ω 是 2 中的有界开集,边界 ∂Ω 属于 C3 ,若 Φ ∈C1 (Ω) ,且 Φ L2 (Ω) 1 ,则
对于 0 < ε 1 ,(1)存在唯一解 F= µ + ε µ f ,f 满足
v ⋅∇x f + ε2
1 µ
Φ

∇v
µ f + ε −1Lf = Γ( f , f ) + εΦ ⋅ v
定常数。
估计(7)的右边需要下列 Sobolev-Gagliardo-Nirenberg 不等式:
对 于 1 ≤ p ≤ N 和 有 界 1 区 域 Ω ⊂ N , 若 u ∈W 1, p (Ω) , 则 对 任 意 的 p ≤ p* =Np , 有
N−p
1
∫ ( )
u
p*
p*
≤C
N, P,Ω
= u ( x) 0= ,θ ( x) 0 ,在 ∂Ω 上,其中 u = (u1,u2 ) 。
2. 预备知识和主要引理
首先引进一些基本记号:
1
(∫ ) ( ) 我们定义下列范数: ⋅ ≡ ν 1 2 ⋅ =
ν
2
( ) ν 1 2 ⋅ L2 Ω×3
;= f p,±
= f 1γ ± p
f
γ
x,v 1± p dγ
Open Access
1. 引言
本文考虑如下玻尔兹曼方程的不可压缩极限:
ε −1v ⋅ ∇xF + εΦ ⋅ ∇vF =ε −2Q ( F, F ) , ( x,v) ∈ Ω × 3 ,
{ } ∫ ( ) F x,v n(x)⋅v <0 ( = 2πµ n(x)⋅u >0 F x,u) n ( x) ⋅ u du , x ∈ ∂Ω ,
γ+
fψ −
γ−
fψ .
(7)
∫∫ ( = −ε −1τ Ω ψ ×3 L I − P) f ∫∫ + ψ Ω×3 g
令a−
f
= a ,则 P f =
a
+
v⋅b
+
c
v
2 −3 2
µ ,分以下三步证明。
第一步 估计 c,对于足够小的 ε > 0 ,我们将证明如下结论:
( ) c 4 < o(1) Pf
2
定义 2.1:记 Ω × 3 的边界为 γ = ∂Ω × 3 ,我们将 γ 分为以下三种情况:
{ } = γ + ( x,v) ∈∂Ω × 3 : n ( x) ⋅ v > 0 ,
{ } = γ − ( x,v) ∈∂Ω × 3 : n ( x) ⋅ v < 0 ,
{ } = γ 0 ( x,v) ∈∂Ω × 3 : n (= x) ⋅ v 0 .
其中,对于 ( x,v) ∈ Ω × 3 , tb ( x= ,v) inf {t > 0 : x − tv ∉ Ω} , xb ( x,v)= x − tb ( x,v)v ∉ Ω 。
( ) ( ) 引理 2.3:假设 Φ ∈ L∞ Ω × 3 ,存在 g ∈ L2 Ω × 3 使得
∫∫Ω×3 g ( x,v) µdxdv = 0 ,

u W1,p (Ω) ,
并且 W 1, p (Ω N = 2 时,我们想要 p* = 2 ,所以由 1 − 1 ≤ 1 ,得 p ≥ 1 ,这里我们取 p = 4 ,则对于任意的
(2)
且 f 是下述方程在分布意义下的解:
λ
+ (1 −τ )ε −1ν

1 2
ε


v
f
+ v ⋅∇x f
+ ε 2Φ ⋅ ∇v f
+ ε −1τ Lf
= g ,在 Ω × 3 上,
(3)
f− = Pγ f ,在 γ − 上,
则对足够小的 λ ≥ 0 以及趋于 1 的τ ∈[0,1] ,
( )
2
{ } { } ∫∫ ∫ ∫ Ω×3 v ⋅ ∇x f + ε 2Φ ⋅ ∇v f ψ + v ⋅ ∇xψ + ε 2Φ ⋅ ∇v= ψ f γ + fψ − γ − fψ .
结合上式以及 LPf = 0 ,由(3)得到
∫∫ ∫ ∫ ω Ω×3 τ fψ − v ⋅ f ∇xψ − ε 2 f Φ ⋅ ∇vψ +
Keywords
Boltzmann Equation, Navier-Stokes-Fourier Equation, Hydrodynamic Limit
二维玻尔兹曼方程的不可压缩极限
高婷婷
华南理工大学,数学学院,广东 广州
收稿日期:2020年1月29日;录用日期:2020年2月17日;发布日期:2020年2月24日
DOI: 10.12677/pm.2020.102019
129
理论数学
高婷婷
( ) ( ) 最后,当 ε → 0 ,f在 L2 Ω × 3 上弱收敛于 f1 = u ⋅ v + θ v 2 − 5 2 µ ,而 ( p,u,θ ) 是具有狄利克
雷边界条件以及外力场 Φ 的稳态INFS方程的唯一弱解: u ⋅ ∇xu + ∇x p = σ∆u + Φ , ∇x ⋅ u =0 ,在 Ω 上, u ⋅ ∇xθ = κ∆θ ,在 Ω 上,
v ∈ 2,vˆ ∈ 。定义 µ =
1
− |v|2
3e 2
,B (V ,ω=)
V ⋅ω 。Φ ( x) = (Φ1,Φ2,0) 表示给定的外力,F 是稀薄
(2π)2
气体分子的分布函数, ε 表示气体分子的平均自由程。
文献[2] [3]在一些先验假设下研究了玻尔兹曼方程的不可压缩的 Navier-Stokes-Fourier (简称 INFS)极 限。对于重整化解的收敛极限,完整证明由[4]给出。对于稳态 Boltzmann 的研究则比较少,正如文献[5] 所指 出 , 尽管 稳态 Navier-Stokes-Fourier 方 程在 应 用中 很 重 要, 但 从稳 态 Boltzmann 推 导出 稳态 Navier-Stokes-Fourier 方程一直是一个待解决的重要问题。最近文献[6]通过 L2-L∞方法结合 L6 估计证明了 三维稳态玻尔兹曼方程的不可压缩极限。在此基础上我们研究二维稳态玻尔兹曼方程的不可压缩极限。 需要指出,相比于三维情形,对于二维稳态玻尔兹曼方程,我们需要 L2-L∞方法结合新的 L4 估计来研究, 由此导致了不同估计和困难。
µ ,在 Ω 上,
f = Pγ f ,在 γ − 上,
并且
( ) f
+
2
Pf
+ ε −1
4
(I − P) f
ν
+ ε −1 2
1 − Pγ
f
+ε1 2
2,+
wf
∞ 1,
其中 Pγ f = 2π µ (v) ∫n(x)⋅u>0 f (u) µ (u){n ( x) ⋅ u}du , w(v) = eβ v 2 , 0 < β 1 。
p,
( ) ( ) = dγ
n(x)⋅v
dS ( x)dv ;⋅
表示

L∞
Ω × 3
范数或者 L∞ (Ω) 范数;⋅ ∞ 表示 L∞
∂Ω × 3
范数或者 L∞ (∂Ω)
范数。
X < Y 等价于 X < CY ,C 是与 X ,Y 无关的常数;定义 P=f a µ + v ⋅ b µ + c v 2 − 3 µ 。
4 + ε −1 ( I − P) f
+
ν
(1 − P) f
+
4
1 − Pγ
f+ 2,+
g. ν2
(8)
( ) 这里择测试函数= ψ ψ c,4 ≡ v − βc µv ⋅ ∇ ϕx c,4 ( x) ,其中 −∆xϕc,4 ( x ) = c3 ( x ) ,ϕc,4 ∂Ω = 0 ,βc 是一待
文章引用: 高婷婷. 二维玻尔兹曼方程的不可压缩极限[J]. 理论数学, 2020, 10(2): 128-138. DOI: 10.12677/pm.2020.102019
关键词
玻尔兹曼方程,Navier-Stokes-Fourier方程,流体动力学极限
高婷婷
Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). /licenses/by/4.0/
相关文档
最新文档