浅话边界条件与初始条件
偏微分方程的定解条件与解的存在唯一性
偏微分方程的定解条件与解的存在唯一性偏微分方程（Partial Differential Equation, 简称PDE）是数学领域中的重要研究对象，广泛应用于物理学、工程学、金融学等领域。
在求解偏微分方程时，我们需要考虑定解条件，以确保解的存在和唯一性。
本文将探讨偏微分方程的定解条件，并讨论解的存在唯一性。
一、偏微分方程的定解条件在求解偏微分方程之前，我们需要明确的是问题的定解条件。
定解条件是指在区域Ω上关于未知函数u及其偏导数的附加条件。
常见的定解条件包括初始条件和边界条件。
1. 初始条件（Initial Condition）初始条件是在区域Ω的某个子集Ω₀上给定的函数值及其偏导数，常用符号表示为u(x, t₀) = g(x, t₀)，其中g(x, t₀)为已知函数，t₀为给定的初始时间。
2. 边界条件（Boundary Condition）边界条件是在区域Ω的边界上给定的函数值及其偏导数，常用符号表示为u(x, t) = f(x, t)，其中f(x, t)为已知函数。
在一些情况下，还需要考虑特殊的边界条件，如周期性边界（Periodic Boundary Conditions）和运动边界（Moving Boundary Conditions）等。
二、解的存在唯一性偏微分方程的解的存在唯一性是指在给定的定解条件下，方程是否有解以及解是否唯一。
1. 解的存在性对于某些偏微分方程，我们可以通过适当的数学工具（如变分法、分离变量法、线性化等）证明其存在解。
然而，并非所有的偏微分方程都具备解的存在性，存在着某些无解的情况。
因此，对于求解偏微分方程问题，我们需要首先考虑其解的存在性。
2. 解的唯一性在一些情况下，即使偏微分方程存在解，其解也不一定是唯一的。
对于线性偏微分方程，我们可以通过使用变分法或利用极大模原理来证明解的唯一性。
而非线性偏微分方程的唯一性则比较复杂，通常需要借助于更加深入的分析和数学工具。
初始条件与边界条件
u 0. n S
过程可知，有边界条件
当物体与外界接触的表面 S 上各单位面积在单位 时间内流过的热量已知时，由傅立叶定律，在 S dQ u k 上有 dSdt n，这表明温度沿外法线方向的方 向导数是已知的，故边界条件可以表示为
u M,t n S
x 0
二阶偏微分方程
2u 2u 2u u u a11 2 2a12 a 22 2 b1 b2 cu f xy x y x y
可简写为 L[u] f . 定解条件
u x
g
x 0
可简写为 B[u ] g.
叠加原理 1 若 ui 满足线性方程
热传导方程的Cauchy问题
utt a 2 uxx 0 u |t 0 ( x ) u | ( x ) t t 0
( x , t 0) ( x )
波方程的Cauchy问题
由偏微分方程和相应边界条件构成的定解问题称 为边值问题。
§1.3
定解问题的提法
初始条件和边界条件都称为定解条件。 定解问题是指偏微分方程和相应定解条件的结合体。 偏微分方程和相应初始条件构成的定解问题称为初 值问题或者柯西(Cauchy)问题。
ut a 2uxx 0 u |t 0 ( x ) ( x , t 0) ( x )
定解条件。也就是说，当定解条件有微小变动时，
引起解的变动是否足够小。若是，则称解是稳定的，
否则称解是不稳定的。
例 设弦的两端固定于x=0 和x=l，弦的初始位移 如下图，初速度为零，求弦满足的定解问题。 解：
2u 2u a2 0 x l , t 0; 2 2 0 t x u u x l 0; x 0 l x, 0 x u 2 , 0 ut 0 t t 0 l x, l x l 2
Airpak模拟的边界条件和初始条件
Airpak模拟的边界条件和初始条件各位同学，大家好，我是七师兄，今天我们来学习Airpak高级班的第三节课。
在第二节课中，我们着重介绍了在CFD模拟过程中要遵循的一些控制方程。
那么这几课我们就来看下，如何求解这些方程。
我们在模拟中用到的这些方程比如，质量守恒，能量守恒，动量守恒等，这些方程可以组成一个方程组，但是这些方程组并不能很好的求出解来，要想求出解，也就数数学上所说的让方程组闭合,必须有所谓的定解条件才能封闭上述方程组,才能得出问题的解。
对于一个一般性的非稳态问题,定解条件包括边界条件和初始条件。
边界条件和初始条件，他是让我们方程闭合的前提。
首先我们来看下边界条件，我们在工程热力学中，学过三大边界条件，那么其实下面，我们说的也就是三大边界条件。
1.边界条件1)给出变量τ中的值,如壁面的温度,非滑动壁面的速度分量为零等。
2)给出τ沿某方向的导数值,如已知壁面的热流量。
3）给出时间和传热量的关系式,如通过表面传热系数以及周围流体温度而限定壁面的换热量等。
那么这里讲的是理论的边界条件，那么在我们CFD模拟的时候，具体的边界条件有哪些呢，我们来看下。
在CFD模拟计算时,基本的边界类型包括以下几种:(1)入口边界条件入口边界条件：就是指定入口处流动变量的值。
常见的入口边界条件有速度、压力、质量流量入口边界条件。
速度入口:用于定义流动速度和流动入口的流动属性相关的量。
这一边界条件适用于不可压缩流,如果用于可压缩流会导致非物理结果,这是因为它允许驻点条件浮动。
应注意不要让速度入口靠近固体妨碍物,因为这会导致流动入口驻点属性具有太高的非一致性。
压力入口:用于定义流动入口的压力及其他标量属性。
它既适用于可压流也可用于不可压流。
压可用于压力已知但是流动速度或速率未知的情况。
这一情况可用于很多实际问题,如浮力驱动的流动。
压力入口边界条件也可用来定义外部或无约束流的自由边界。
质量流量入口：用于已知入口质量流量的可压缩流动(2)出口边界条件压力出口边界条件:压力出口边界条件需要在出口边界处指定表压。
1-2_初始条件与边界条件chen
练习题： 考虑长为 l 的均匀杆的导热问题． 若 １．杆的两端温度保持零度；
２．杆的一端为恒温零度，另一端绝热． 写出杆在上面两种情况下的边界条件． 答案： u
= 0, u = 0;
u = 0, u x = 0;
x=0
x=l
x=0
x=l
数学物理方程与特殊函数 主页 上一页 下一页
退出
退出
偏微分方程要给出 n 个初始条件才能确定一个特 解．
初始条件的个数的确定：关于时间 t 的 n 阶
边界条件的个数的确定：关于空间变量 x 的 n 阶 偏微分方程不一定要给出 n 个边界条件．
数学物理方程与特殊函数 主页 上一页 下一页
退出
练习题： 若 考虑长为 l 的均匀杆的导热问题． １．杆的两端温度保持零度； ２．杆的一端为恒温零度，另一端绝热． 写出杆在上面两种情况下的边界条件． 答案：
∂u = f . 外法线方向的方向导数,即 ∂n s 2
第二类边界条件.
第一类边界条件. 三是在边界 S 上给出了未知函数 u 及其沿 S
的外法线方向的方向导数某种线性组合的
⎛ ∂u ⎞ + σ u ⎟ = f3 . 值,即 ⎜ ∂n 第三类边界条件. ⎝ ⎠s
数学物理方程与特殊函数 主页 上一页 下一页
M
初始条件:
u(M ,t)
t=0
= ϕ (M )
退出
数学物理方程与特殊函数 主页 上一页 下一页
泊松方程与拉 普拉斯方程
泊松方程与拉 普拉斯方程都是描 述稳恒状态的,与初 始状态无关,不提初 始条件.
数学物理方程与特殊函数 主页 上一页 下一页
退出
初始条件和边界条件
11. 初始条件和边界条件
它们也是数值方法成败的关键
纷繁复杂的天体现象满足几乎同一方程，这很大程度上归因于定解条件：初边值
举例
Chen & Shibata (2000)
解释日冕物质抛射
Chen et al. (2002)
解释EIT波
（一）定态解析解
2
2
202
220)1(2,)1()1(2ky x k kx B B ky x k ky B B y x ++−=+++=势场
2
02
0cos ,sin r B B r B B r ϕϕϕ−==3
3
0330sin ,cos 2r
a B B r a B B r θ
θθ==直角坐标系下的拱形场
柱坐标系下
线偶极子场
球坐标系子午面偶极子场
（三）观测值作初条
通常是部分观测量，如太阳大气温度的分布。
v
(5.4) 特征线方法在数值计算中的应用
借助特征线方法一阶偏微分方程组化为特征形式后，方程维数减1，因此该方法对一维非定常问题的处理特别有效，因为原方程退化为常微分方程，可沿特征线直接积分。
但对多维问题，此优点并不突出，且计算量繁重。
偏微分方程中的边界条件与初始条件
偏微分方程中的边界条件与初始条件在偏微分方程的求解过程中，边界条件和初始条件是非常重要的。
边界条件定义了方程在空间边界上的行为，而初始条件则规定了方程在时间初始时刻的状态。
这两个条件的正确选择和准确给定对于得到准确的解是至关重要的。
一、边界条件的选择和应用在偏微分方程中，边界条件通常指定在空间边界上。
根据不同的问题和方程形式，可以有不同类型的边界条件。
1. Dirichlet 边界条件：Dirichlet 边界条件指定了方程解在边界上的具体数值。
例如，在一个热传导问题中，可以通过指定边界上的温度值来应用Dirichlet 边界条件。
2. Neumann 边界条件：Neumann 边界条件指定了方程解在边界上的法向导数。
例如，在一个扩散问题中，可以通过指定物质的流量或粒子的扩散速率来应用Neumann边界条件。
3. Robin 边界条件：Robin 边界条件是一种将 Dirichlet 边界条件和Neumann 边界条件结合在一起的条件形式。
它通常用于描述具有热传导和对流作用的问题。
在实际问题中，选取合适的边界条件需要根据问题的物理特性进行合理的选择。
对于特定的问题，可能需要根据实际情况进行数值模拟或实验来确定最合适的边界条件。
二、初始条件的选择和应用初始条件是指在时间初始时刻系统状态的描述。
在时间发展过程中，系统的状态随着时间的推移而变化，初始条件的准确给定对于求解偏微分方程的初值问题非常重要。
对于偏微分方程中的初始条件，一般需要给定系统在时间初始时刻的值或概率分布。
例如，在一个传热问题中，可以通过给定系统在时间初始时刻的温度分布来应用初始条件。
同样，初始条件的选择和给定需要根据具体问题的特性进行合理的确定。
在实际应用中，初始条件的准确性和合理性会直接影响到模拟结果的可靠性和有效性。
三、边界条件与初始条件的影响边界条件和初始条件的选择和应用直接影响着偏微分方程求解结果的准确性和可靠性。
不恰当或不准确的边界条件和初始条件可能导致解的不合理或不符合实际的情况。
数理方程初始条件与边界条件
2
解记为 u1 ( x, t )
（可由达朗贝尔公式给出）
utt a 2u xx f ( x, t ), t 0, x , （C） 解记为 u2 ( x, t ) u ( x,0) 0, ut ( x,0) 0.
由叠加原理可知
u( x, t ) u1 ( x, t ) u2 ( x, t ).
• 定解问题：泛定方程加上适当的定解条件就构成一个定 解问题，即定解问题=泛定方程+定解条件。
1.3
定解条件
1、初始条件——描述系统的初始状态
A、 波动方程的初始条件
u |t 0 ( x) 系统各点的初位移 u ( x) t 系统各点的初速度 t 0
B、热传导方程的初始条件 初始时刻的温度分布： u(M , t ) |t 0 (M ) C、泊松方程和拉普拉斯方程的初始条件 不含初始条件，只含边界条件条件
哈密尔顿算子或梯度算子，读作nabla ˆ ˆ ˆ i j k x y z
与梯度算子有关的场论运算
gradu u
divA A
rotA A
2 2 2 ห้องสมุดไป่ตู้ 拉普拉斯算子 3 2 2 2 2 x y z
作变量代换
x x at
u a 0
解为：u f ( x at)
f
为任意函数
7
举例（未知函数为二元函数）
3.
2u 0 xt
解为： u g ( x) h(t )
2 2u u 2 a 0 4. 2 2 t x
变换
x at x at
2、边界条件——描述系统在边界上的状况
数值天气预报第四章_初始条件与边界条件
兰州大学大气科学学院
初始条件与边界条件
2、初始化发展历程(续)
Miyakoda and Moyer（1968），Nitta and Hovermale提出新的初始化方法，称为动力初 始化。基本原理是
平衡方程(4.3)可以写为
( ) ∇2ψ =
1 f
2.平衡初值
平衡初值是采用平衡方程作为风场和气压场之间的协调
关系。平衡方程为：
fζ −βu +2J (u,v) =∇2Φ
(4.2)
兰州大学大气科学学院
初始条件与边界条件
一、静力初始化(续)
假定水平无辐散，引入流函数 ψ，则平衡方程为：
( ) f
∇2ψ
+ ∇f
⋅ ∇ψ
+
2
ψ xxψ
yy
−ψ
2 xy
初始条件与边界条件
2、初始化发展历程
最早的初始化过程基于准地转理论。Charner(1955) 建议用非线性平衡由分析的位势高度场计算流函数， 他认为这种平衡的初始状态可以有效地抑制惯性重力 波。
Hinkelmunn（1959）和Phillips（1960）论证了仅仅 利用非线性平衡方程还不足以达到上述目的。他们建 议运用 ω方程给出初始的速度位势χ。这些早期的 准地转初始化过程可以看成是准地转约束的扩展。在 求解ω方程，特别是非线性平衡方程时有一些技术困 难，最大的困难是“椭圆型”问题，一般采取选代算法 来解决。准地转初始化技术在中高纬度地区相当成 功，但在低纬低区就不适宜。
数理方程：第2讲典型方程的定解条件
弹性力 k u xl
张力
T u x xl
u
u
x
xl
0
( T )
k
(2) 热传导问题（端点自由冷却）
散失的热量
dQ1 h(u u1)dSdt
内部流到边界的热量
dQ2
k
u dSdt n
dQ1 dQ2 k nu h(u u1 )
即
(u
u )
x xl
u1
( k )
h
§3 定解问题
x
, t
0)
的解等于问题(I)和问题(II)的解之和
(I)
utt u t0
a2uxx 0
( x), ut
t0
(
(x)
x
, t
0)Байду номын сангаас
(II) utt a2uxx f ( x, t ) u t0 0, ut t0 0
( x ,t 0)
如
Lu
a11
2u x 2
2a12
2u xy
a22
2u y 2
b1
u x
b2
u y
c
u
u x
x0
二阶偏微分方程
a11
2u x 2
2a12
2u xy
a22
2u y 2
b1
u x
b2
u cu y
f
可简写为
L[u] f
叠加原理 1 若ui 满足线性方程 L[ui ] fi ，i 1,2,, n
➢包含初值条件和边界条件的定解问题称为混合问题 (初边值问题)
uutt
0
a2(uxx
(x, y
u yy ,z)
微分方程 重点条件
微分方程重点条件1 前言微分方程是数学中一门重要的学科，被广泛应用于物理、工程、生物等领域。
在求解微分方程时，我们需要考虑各种条件，这些条件被称为微分方程的重点条件。
本文将介绍微分方程的重点条件，并给出一些例子。
2 初始条件初始条件是微分方程求解中最常见的一种条件。
对于一个一阶常微分方程y'=f(x,y)，它的初始条件为y(x0)=y0，其中x0和y0是已知的常数。
这个条件告诉我们在解的曲线上点（x0,y0）是曲线的起点。
对于高阶微分方程，初始条件也是类似的。
举例：考虑一阶常微分方程y'=2x，初始条件y(0)=1。
我们可以求得其解为y=x^2+1。
而这个初始条件告诉我们这个解曲线经过点（0,1）。
3 边界条件边界条件是微分方程在一定边界范围内的条件，通常为两种形式：Dirichlet边界条件和Neumann边界条件。
Dirichlet边界条件指定微分方程在边界上的函数值，形式为y(a)=C1，y(b)=C2。
Neumann边界条件指定微分方程在边界上的导数值，形式为y'(a)=C1，y'(b)=C2。
举例：考虑二阶常微分方程y''+y=x，Dirichlet边界条件为y(0)=0，y(\pi/2)=1。
我们可以求得y=sin(x)。
4 周期条件周期条件是微分方程在某个周期内成立的条件，通常用于描述周期性的现象。
举例：考虑一阶常微分方程y'+y=sin(x)，周期条件为y(0)=y(2π)。
我们可以求得一个解为y=-cos(x)+A*sin(x)，其中A为常数。
5 约束条件约束条件是微分方程所满足的某些特殊条件，通常会限制解的取值范围或形式。
常见的约束条件有非负性、有界性、单调性等。
举例：考虑一阶常微分方程y'=y^2，约束条件为y(0)=1。
我们可以求得其解为y=1/(1-x)，但这个解不满足约束条件y≥0，因此实际上只有在x<1时才是合法解。
第9章边界初始条件
ui |t =t0 =ψi (x, y, z) t
(i=1,2,…,n) …
式中ψ 的已知函数。 式中ψi是x，y，z的已知函数。 对于稳定的物理过程 稳定的物理过程， 对于稳定的物理过程，此时刻的状态和彼时刻 的状态都是相同的。因此，便不须有(也提不出) 的状态都是相同的。因此，便不须有(也提不出)初 始条件了。 始条件了。
P P |x=0 = |x=xmax = 0 x x
对于初始的静水 环境，则有： 环境，则有：
21
2
11
P |t =t0 = ρ f gZ
1
应用流体、固体、渗流、热力学方面的知识， 应用流体、固体、渗流、热力学方面的知识，推 导出遵守质量守恒定律、动量定律、 导出遵守质量守恒定律、动量定律、能量守恒定律的 数学表达式，称为基本方程式。 数学表达式，称为基本方程式。 基本方程式 但基本方程式的个数往往比未知量的个数少， 但基本方程式的个数往往比未知量的个数少，不 封闭，需要补充一些方程，使其构成为封闭方程组。 封闭，需要补充一些方程，使其构成为封闭方程组。 后者须由实验事实或物理概念的基本假设出发来构造 或选择具体过程的物理模型，如流动是等熵的、 或选择具体过程的物理模型，如流动是等熵的、材料 物理模型 是弹塑性的、渗流是等温的等等。 是弹塑性的、渗流是等温的等等。这些热力学状态和 特殊物理假定的数学表达式，就是所谓 “补充方程”。 特殊物理假定的数学表达式， 补充方程”
如下图所示，区域 为正常压力区， 为超压区。 如下图所示，区域2为正常压力区，1为超压区。 的上边界处， gZ（静水压力）； ）；假 在1的上边界处，有P|11 =ρfgZ（静水压力）；假 11 若沉积底界为不渗水基底，则有： (P若沉积底界为不渗水基底，则有：k(P- ρfgZ) =0 。若 盆地边界为无水平流动侧面， 盆地边界为无水平流动侧面，则：
PFC2D学习笔记之边界条件与初始条件
PFC2D学习笔记之边界条件与初始条件边界条件墙边界一般，颗粒组是在一组墙内创建并压缩。
这些墙也可以当作边界约束，按一定速度运动监测其反力；或保持反力一定，伺服控制其速度。
但是不能在墙上直接施加力。
当墙为多段线时，要考虑到墙与球在线段转角处的接触。
当转角处为凹时，两段墙可能会同时与一个球接触，但PFC2D只允许每个墙与球的接触最多一个，因此要在此处将墙断开，设置成两个墙；当转角为凸时，两段墙不肯能同时与一个球接触，因此不需要将墙断开。
墙也可以是圆形或弧或点。
颗粒边界可以创建一个颗粒串，并使用这些颗粒串作为边界条件。
固定速度的颗粒边界用 FISH函数得到边界颗粒，并固定其平动自由度，然后删除墙，并施加加速度。
从而可以得到内部颗粒的速度。
边界颗粒的速度保持初始值不变；该模拟对应与应变控制式试验。
有两种方法可以得到边界颗粒。
一种是，遍历每个球的接触列表并探测与墙接触的球，作为边界；另一种是，使用RANGE，将落在一个空间范围内的所有颗粒作为边界颗粒。
第一种方法的核心程序为：bp = ball_headloop while bp # null ; scan all ballssectioncp = b_clist(bp)loop while cp # null ; scan ball’s contactsif c_nforce(cp) # 0.0 thenb2 = c_ball2(cp)if pointer_type( b2 ) = 101 then ; b2 is a wallb_xfix(bp) = 1 ; fix original ball in x,yb_yfix(bp) = 1b_color(bp) = 1 ; identify by colorexit section ; all done for this ballend_ifend_ifif c_ball1(cp) = bpcp = c_b1clist(cp)elsecp = c_b2clist(cp)end_ifend_loopend_sectionbp = b_next(bp)end_loop第二种方法参考RANGE的用法，较为简单，此处略去。
浅话边界条件与初始条件
浅话边界条件与初始条件边界条件在说边界条件之前，先谈谈初值问题和边值问题。
初值和边值问题：对一般的微分方程，求其定解，必须引入条件，这个条件大概分两类---初始条件和边界条件，如果方程要求未知量y(x)及其导数y′(x)在自变量的同一点x=x0取给定的值，即y(x0 )=y0，y′(x0)= y0′，则这种条件就称为初始条件，由方程和初始条件构成的问题就称为初值问题;而在许多实际问题中，往往要求微分方程的解在在某个给定的区间a ≤ x ≤b的端点满足一定的条件，如y(a) = A , y(b) = B 则给出的在端点(边界点)的值的条件，称为边界条件，微分方程和边界条件构成数学模型就称为边值问题。
三类边界条件:边值问题中的边界条件的形式多种多样，在端点处大体上可以写成这样的形式，Ay+By'=C，若B=0，A≠0,则称为第一类边界条件或狄里克莱(Dirichlet)条件;B≠0,A=0，称为第二类边界条件或诺依曼(Neumann)条件;A≠0,B≠0，则称为第三类边界条件或洛平 (Robin)条件。
总体来说，第一类边界条件：给出未知函数在边界上的数值;第二类边界条件：给出未知函数在边界外法线的方向导数;第三类边界条件：给出未知函数在边界上的函数值和外法向导数的线性组合。
对应于comsol，只有两种边界条件：Dirichlet boundary(第一类边界条件)—在端点，待求变量的值被指定。
Neumann boundary(第二类边界条件)—待求变量边界外法线的方向导数被指定。
再补充点初始条件：初始条件，是指过程发生的初始状态，也就是未知函数及其对时间的各阶偏导数在初始时刻t=0的值.在有限元中，好多初始条件要预先给定的。
不同的场方程对应不同的初始条件。
总之，为了确定泛定方程的解，就必须提供足够的初始条件和边界条件.边界条件与初始条件是控制方程有确定解的前提。
边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。
初始条件边界条件和约束.完美版PPT
边界条件
局部坐标系的定义： *DEFINE_COORDINATE 定义一局部坐标系 指定坐标系ID号(CID) *DEFINE_COORDINATE_NODES »3节点: 局部坐标系原点，沿x轴, 局部 x-y 平面内 *DEFINE_COORDINATE_SYSTEM »三点 的x, y, z 坐标(与 NODES方式一样) *DEFINE_COORDINATE_VECTOR »2 个矢量: 局部 x轴, 局部x-y平面内矢量
初始条件
• 起爆点和动量 • 初始应力/应变 • 初始温度 • 初始速度 缺省状态下初始应力、温度和速度为零 边界条件高于初始条件
初始条件
*INITIAL _DETONATION 和*INITIAL_MOMENTUM 用于模拟施加在体单元上的一种脉冲载荷 起爆点:引爆炸药材料(parts) »
动量: 单元上施加一个初始动量
第章初始条件边界条件和约束
概述
• 若FEA模型已具有节点，单元，材料特性和 parts等。则可以： »施加：初始条件、边界条件、载荷和约束 »必须具有：boxes，曲线，sets，矢量等
• 当parts间发生碰撞，或与其它问题撞击时： »刚墙 »接触 (第三章)
• 怎样测定载荷/动量 »接触力 »横截面分析
载荷
• 目的：定义施加的“力 ” »梁 »体 »爆炸效果 »热，温度 »节点和刚体 »壳（压力）
• 避免单点集中载荷 »物理上无意义 »防止沙漏模式 • 避免阶跃载荷 • 要求一载荷曲线 • 载荷能缩放
载荷
*LOAD_BEAM • 沿梁单元局部轴(r,s,t)定义牵 引载荷 • 梁单元或系列梁单元 • 每单位长度上的力
爆炸球的例子
爆炸球的例子
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
浅话边界条件与初始条件
边界条件
在说边界条件之前，先谈谈初值问题和边值问题。
初值和边值问题：
对一般的微分方程，求其定解，必须引入条件，这个条件大概分两类---初始条件和边界条件，如果方程要求未知量y(x)及其导数y′(x)在自变量的同一点x=x0取给定的值，即y(x0 )=y0，y′(x0)= y0′，则这种条件就称为初始条件，由方程和初始条件构成的问题就称为初值问题;
而在许多实际问题中，往往要求微分方程的解在在某个给定的区间a ≤ x≤b 的端点满足一定的条件，如y(a) = A , y(b) = B则给出的在端点(边界点)的值的条件，称为边界条件，微分方程和边界条件构成数学模型就称为边值问题。
三类边界条件:
边值问题中的边界条件的形式多种多样，在端点处大体上可以写成这样的形式，Ay+By'=C，若B=0，A≠0,则称为第一类边界条件或狄里克莱(Dirichlet)条件;B≠0,A=0，称为第二类边界条件或诺依曼(Neumann)条件;A≠0,B≠0，则称为第三类边界条件或洛平(Robin)条件。
总体来说，
第一类边界条件：给出未知函数在边界上的数值;
第二类边界条件：给出未知函数在边界外法线的方向导数;
第三类边界条件：给出未知函数在边界上的函数值和外法向导数的线性组合。
对应于comsol，只有两种边界条件：
Dirichlet boundary(第一类边界条件)—在端点，待求变量的值被指定。
Neumann boundary(第二类边界条件)—待求变量边界外法线的方向导数被指定。
再补充点初始条件：
初始条件，是指过程发生的初始状态，也就是未知函数及其对时间的各阶偏导数在初始时刻t=0的值.在有限元中，好多初始条件要预先给定的。
不同的场方程对应不同的初始条件。
总之，为了确定泛定方程的解，就必须提供足够的初始条件和边界条件.边界条件与初始条件是控制方程有确定解的前提。
边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。
对于任何问题，都需要给定边界条件。
初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况，对于瞬态问题，必须给定初始条件，稳态问题，则不用给定。
对于边界条件与初始条件的处理，直接影响计算结果的精度。
在CFD模拟时，基本边界条件有：
1流动进口边界
包括速度进口边界，压力进口边界，质量进口边界(可压流动)。
在使用流动进口边界时，需要涉及到某些流动参数，如绝对压力，湍动能及耗散率，这些参数要做特殊考虑。
关于参考压力，在流场数值计算中，压力总是按相对值表示的，实际求解的压力并不是绝对值，而是相对于进口压力而言的。
在有些情况下，可以通过设定进口压力为0，求解其他点的压力。
还有时，为了减小数字截断误差，往往故意抬高或降低参考压力场的值，可使其余各处的计算压力场与整体数值计算的量级相吻合。
2流动出口边界
一般选在离几何扰动足够远的地方来施加。
在这样的位置，流动是充分发展的，沿流动方向没有变化。
该边界只有在进入计算域的流动是以进口边界条件给定时才使用，而且在只有一个出口的计算域中使用。
3给定压力边界
在流动分布的详细信息未知，但边界的压力值已知的情况下，使用恒压边界条件。
应用该边界条件的典型问题包括：物体外部绕流，自由表面流，自然通风及燃烧等浮力驱动流和有多个出口的内部流动。
4壁面边界
对于湍流计算，在近壁面区，要用壁面函数来做处理。
壁面函数法的前提是：流动平行于壁面而且只在垂直壁面的方向变化，流动方向不存在压力梯度，壁面处不存在化学作用，雷诺数足够大。
若这些条件之一不满足，则使用壁面函数法的预测精度会大大降低，甚至完全不可用。
5对称边界
在对称边界上，垂直边界的速度取为0，而其他物理量的值在该边界内外是相等的，即计算域外紧邻边界的节点的值等于对应的计算域内紧邻边界的节点的值。
6周期性(循环)边界
流出循环边界出口的所有流动变量的通量等于进入循环边界的对应变量的通量。
使用边界条件要注意的地方：
①选择合理的边界条件，边界条件的4种组合方式：只有壁面，壁面、进口和至少一个出口，壁面、进口和至少一个恒压边界，壁面和恒压边界。
②为得到准确的结果，出口边界位置的选取要位于最后一个屏障物后10
倍于屏障高度或更远的位置。
③近壁面网格。
要保证壁面函数法有效，就需要使离壁面最近的一内节点位于湍流的对数律层中，即y+必须大于11.63，这就相当于给最靠近壁面的网格到壁面的距离△yp设定了一个下限。
初始条件
在瞬态问题中，除了要在计算开始之前初始化相关的数据外，不需要其他特殊处理。
给定初始条件时要注意的是：要针对所有计算变量，给定整个计算域内各单元的初始条件;初始条件一定是物理上合理的，要靠经验或实测结果。
来源：元计算官网。