大气科学专业流体力学第二章基本方程

大气科学专业流体力学第二章基本方程1第一页，共一百三十七页，2022年，8月28日第二章 基本方程

流体运动同其他物体的运动一样，同样遵循质量守恒、动量守恒和能量守恒等基本物理定律。本章将介绍描述流体运动的连续方程、运动方程和能量方程。2第二页，共一百三十七页，2022年，8月28日主要内容：第一节 连续方程第二节 作用于流体的力、应力张量第三节 运动方程第四节 能量方程第五节简单情况下的N-S方程的准确解第二章 基本方程

3第三页，共一百三十七页，2022年，8月28日第一节 连续方程

连续方程是流体力学的基本方程之一，它是在质量守恒定律在流体力学中的应用。流体运动的连续方程，反映流体运动和质量分布的关系，重点讨论几种不同表现形式的流体连续方程。4第四页，共一百三十七页，2022年，8月28日1、拉格郎日(Lagrange)观点下的流体连续方程Lagrange观点下质量守恒定律：某一流体块（流点）在运动过程中，尽管其体积和形状可以发生变化，但其质量是守恒不变的。拉格郎日型连续方程5第五页，共一百三十七页，2022年，8月28日Lagrange观点下连续方程的物理意义？6第六页，共一百三十七页，2022年，8月28日

对于不可压缩流体，它在流动过程中每个流点的密度始终保持不变，应有，此时流体的连续性方程为：7第七页，共一百三十七页，2022年，8月28日例2-1-1判断下列流体运动是否为不可压缩？8第八页，共一百三十七页，2022年，8月28日利用欧拉控制体积法导出流体的连续方程的微分形式。在空间上选取一无限小的控制体，如图所示。2、欧拉(Euler)观点下的流体连续方程（一）yxz单位时间内通过左侧面流入控制体的流体质量为：单位时间内通过右侧面流出控制体的流体质量为：单位时间内x方向上流体通过控制体的质量净流出量为：9第九页，共一百三十七页，2022年，8月28日类似可得到y、z方向上的表达式，单位时间内通过整个控制体的流体净流出量为：单位时间内，该控制体内的质量减少为：根据质量守恒定律，对于固定的控制体，单位时间内流出控制体的流体质量应等于单位时间内该控制体内质量的减少，由此得到：10第十页，共一百三十七页，2022年，8月28日2、欧拉(Euler)观点下的流体连续方程（二）拉格郎日型连续方程欧拉型连续方程11第十一页，共一百三十七页，2022年，8月28日欧拉型连续方程的物理意义单位体积的流体质量通量12第十二页，共一百三十七页，2022年，8月28日对于流体的定常运动，有流体的连续性方程可写为：可知，在定常运动中，通过任意控制体表面流体质量的净流入量等于零，即单位时间内流出控制体表面的质量等于流进控制体表面的质量。13第十三页，共一百三十七页，2022年，8月28日对于沿流管的定常流动，设流速与截面垂直，且密度和流速在任意截面内为定值，则沿流管的连续方程：14第十四页，共一百三十七页，2022年，8月28日3、具有自由表面的流体连续方程通常把自然界中水与空气的交界面称为水面或水表面。这种因流动而伴随出现的可以升降的水面，在流体力学中称之为自由表面。实际物理现象：当水面向某处汇集时，该处水面将被拥挤而升高；反之，当该处有水向四周散开时，将使得那里的水面降低。水空气交界面15第十五页，共一百三十七页，2022年，8月28日假设流团密度为，考虑流体运动为二维的，即满足：，取流向方向为x轴。设流体自由表面高度为，即h在各处高低不同且可以随时间变化。具有自由表面的流体连续方程的导出：16第十六页，共一百三十七页，2022年，8月28日在流体中，选取一个以为底的长方形柱体，该柱体是一底面固定不动的空间区域，称为控制区。流体可以通过控制区的侧面，沿x轴方向流出、流入该柱体。17第十七页，共一百三十七页，2022年，8月28日经流体柱后侧流入的流体质量应为：同时，经流体柱前侧流出的质量为：考虑柱体内流体的质量为：流入质量=流出质量=18第十八页，共一百三十七页，2022年，8月28日流出质量减去流入质量柱体内的净流出量＝柱体内质量的减少。流出质量=流入质量=柱体内流体的质量减少为：19第十九页，共一百三十七页，2022年，8月28日 ***积分上限h为x，y，t的函数，可变上限的积分规则：对上式两项展开，左端项为：20第二十页，共一百三十七页，2022年，8月28日 ***积分上限h为x，y，t的函数，可变上限的积分规则：右端项为：21第二十一页，共一百三十七页，2022年，8月28日考虑到与z无关，并消掉等式两端公共项可得：22第二十二页，共一百三十七页，2022年，8月28日可以得到：考虑水为不可压缩的，根据连续方程有： 23第二十三页，共一百三十七页，2022年，8月28日讨论时流向仅取x轴。如流向取平面上的任意方向，上式可写为：这就是用自由表面高度所表示的连续方程。进一步有：均匀流体自由表面附近的流体（浅流体）24第二十四页，共一百三十七页，2022年，8月28日具有自由表面的流体连续方程欧拉型连续方程水空气25第二十五页，共一百三十七页，2022年，8月28日具有自由表面的流体连续方程的物理意义？通常流向取平面上的任意方向它是讨论水面波动及简单的大气动力学问题所经常用到的。26第二十六页，共一百三十七页，2022年，8月28日1、作用于流体的力质量力流体的作用力表面力分析对象：流体中以界面包围的体积为的流体块第二节作用于流体的力、应力张量27第二十七页，共一百三十七页，2022年，8月28日质量力1定义：质量力是指作用于所有流体质点的力。如重力、万有引力、电磁力等。2特征：（1）质量力是一种长程力：质量力随相互作用的元素之间的距离的增加而减小，但对于一般流体的特征运动距离而言，质量力均能显示出来。（2）质量力是一种分布力，分布于流体块的整个体积内，流体块所受的质量力与其周围有无其他流体无关。通常情况下，作用于流体的质量力通常就是指重力。28第二十八页，共一百三十七页，2022年，8月28日如果表示单位质量的流体的质量力：其中是作用在质量为的流体块上的质量力。不难看出，可以看做质量力的分布密度。例如：对处于重力作用的物体而言，质量力的分布密度就是重力加速度。29第二十九页，共一百三十七页，2022年，8月28日表面力1定义：表面力是指流体内部之间或者流体与其他物体的接触面上所受到的相互作用力。如流体内部的粘性力和压力、流体与固体接触面上的摩擦力等。30第三十页，共一百三十七页，2022年，8月28日表面力的特征：（1）表面力是一种短程力：源于分子间的相互作用。表面力随相互作用元素之间的距离增加而迅速减弱，只有在相互作用元素间的距离与分子距离同量级时，表面力才显现出来。（2）流体块内各部分之间的表面力是相互作用而相互抵消的。（3）表面力也是一种分布力，分布在相互接触的界面上。31第三十一页，共一百三十七页，2022年，8月28日定义单位面积上的表面力（即：表面应力）为： 其中是作用于某个流体面积上的表面力例如：流体受到的表面力为压力，就是压强。32第三十二页，共一百三十七页，2022年，8月28日

矢量是质量力的分布密度，它是时间和空间点的函数，因而构成了一个矢量场。而矢量为流体的应力矢量，它不但是时间和空间点的函数，并且在空间每一点还随着受力面元的取向不同而变化。所以要确定应力矢量，必须考虑点的矢径、该点受力面元的方向（或者说面元的法向单位矢量）以及时间t。确切地说应力矢是两个矢量（、）和一个标量的函数t。质量力和表面力的比较质量力和表面力有着本质的差别。33第三十三页，共一百三十七页，2022年，8月28日2、应力张量取如图所示的流体四面体元，分析其受力情况。MxyzABC质量为质量力为表面力？？？34第三十四页，共一百三十七页，2022年，8月28日MxyzABC为了区分不同面元所受到的表面力，将应力矢量的下标取其受力面元的外法向方向，并且规定为外法向流体对另一部分流体施加的应力。35第三十五页，共一百三十七页，2022年，8月28日根据牛顿第二定律，MxyzABC根据作用力与反作用力原理36第三十六页，共一百三十七页，2022年，8月28日 根据作用力与反作用力原理，方程可以写成如下形式：37第三十七页，共一百三十七页，2022年，8月28日四面体体积取极限时：上式为作用于流体微元的应力矢量之间的相互关系。38第三十八页，共一百三十七页，2022年，8月28日MxyzABC考虑面元与的关系：PPAMKx39第三十九页，共一百三十七页，2022年，8月28日xyzABC考虑各面元间的关系：40第四十页，共一百三十七页，2022年，8月28日将其在直角坐标系中展开，则有：41第四十一页，共一百三十七页，2022年，8月28日引进应力张量：

42第四十二页，共一百三十七页，2022年，8月28日¤¤对应力分量的下标作如下规定：第一个下标表示受力面元的外法向方向；第二个下标表示受到的应力矢量所投影的方向。应力分量的物理含义：

例2-2-1说明应力、表示的物理含义。43第四十三页，共一百三十七页，2022年，8月28日法应力和切应力

通常应力矢量也可以表示为：切应力法应力44第四十四页，共一百三十七页，2022年，8月28日例2-2-1已知流体中某点的应力张量为试求作用于通过该点，方程为的平面上的法应力和切应力。45第四十五页，共一百三十七页，2022年，8月28日例2-2-2流体中的应力张量为试求位于点（1，2，3）的法应力。46第四十六页，共一百三十七页，2022年，8月28日其中为反映流体粘性的粘性系数或内摩擦系数；而流体与其他物体的粘性系数则称为外摩擦系数。牛顿粘性假设牛顿粘性定律建立了粘性应力与流速分布之间的关系。3、应力张量与流体运动状态间的关系47第四十七页，共一百三十七页，2022年，8月28日广义牛顿粘性假设牛顿粘性定律建立了粘性应力与流速分布之间的关系，但它的不足在于仅仅适用与流体直线运动。牛顿将以上的粘性应力与形变率的关系推广到任意粘性流体运动，即广义牛顿粘性假设：48第四十八页，共一百三十七页，2022年，8月28日说明：根据广义牛顿粘性假设的应力张量计算得到的应力包含了流体压力和流体粘性力两部分即：不可压流体49第四十九页，共一百三十七页，2022年，8月28日牛顿粘性流体的概念：满足牛顿广义粘性假设的流体。给定流体的粘性系数和流体运动流速场，根据牛顿粘性假设，就可以计算得到流体的粘性应力。50第五十页，共一百三十七页，2022年，8月28日51第五十一页，共一百三十七页，2022年，8月28日例2-2-4设速度场为：，试求位于的单位质量长方体（高为）作用在顶面和底面上的粘性应力。52第五十二页，共一百三十七页，2022年，8月28日第三节运动方程流体的运动方程（普遍形式）纳维-斯托克斯（N-S）方程（具体形式）欧拉方程（理想流体的运动方程）静力方程（最简单情形的运动方程）

53第五十三页，共一百三十七页，2022年，8月28日在运动流体中选取一小六面体体元，其边长分别为：为了导出流体的运动方程，首先来分析小体元的受力情况。一、流体的运动方程根据牛顿第二定律： xyz54第五十四页，共一百三十七页，2022年，8月28日x方向质量力分析x方向的质量力55第五十五页，共一百三十七页，2022年，8月28日小体元所受到前后侧面的沿x方向上表面力合力：x方向受到的表面力合力分析周围流体对小体元的六个表面都有表面力的作用后侧面：x?前侧面：xyz56第五十六页，共一百三十七页，2022年，8月28日因此，周围流体通过六个侧面作用于小体元沿x方向的表面力合力为：右左侧面：上下侧面：xyz57第五十七页，共一百三十七页，2022年，8月28日根据牛顿运动定律：小体元受到的合外力等于其质量与加速度的乘积。x方向合力分析单位质量流体在x方向的运动方程方程可以简化为：58第五十八页，共一百三十七页，2022年，8月28日单位质量流体在y方向的运动方程单位质量流体在z方向的运动方程同理可得：59第五十九页，共一百三十七页，2022年，8月28日矢量形式或者：流体运动方程的普遍形式60第六十页，共一百三十七页，2022年，8月28日分析对象：流体中以界面包围的体积为的流体块根据牛顿第二定律流体运动方程的普遍形式61第六十一页，共一百三十七页，2022年，8月28日应用奥—高公式，将以上曲面积分转化为体积分，则有：当曲面面元向内无限收缩时，即体积元趋向于零：62第六十二页，共一百三十七页，2022年，8月28日二、纳维-斯托克斯（Navier-Stokes）方程流体运动方程的普遍形式纳维-斯托克斯方程广义牛顿粘性假设63第六十三页，共一百三十七页，2022年，8月28日流体运动方程的普遍形式广义牛顿粘性假设这就是适合牛顿粘性假设的流体运动N-S方程。法国工程师Navier英国数学家Stokes64第六十四页，共一百三十七页，2022年，8月28日65第六十五页，共一百三十七页，2022年，8月28日定义流体运动学粘性系数，记作。直角坐标系中形式为：对于不可压流体N-S方程简化为：66第六十六页，共一百三十七页，2022年，8月28日其中是单位质量流体的加速度，为单位质量流体所受的质量力。①压力梯度力②粘性（粘滞）力方程物理意义的讨论：①②67第六十七页，共一百三十七页，2022年，8月28日方程右端的第二项，对于某一流体块，有

从而得到：

即为周围流体通过单位质量流点的表面，对其所产生的压力的合力矢量，将其称为压力梯度力。68第六十八页，共一百三十七页，2022年，8月28日仅考虑压力梯度力的作用高压中心低压中心大气的运动形式？69第六十九页，共一百三十七页，2022年，8月28日仅考虑流体作直线运动对于某一流体块其受到的粘滞力U小U大当四周流体速度大于所考虑的流体块时，粘滞力为曳力；当四周流体速度小于所考虑的流体块时，粘滞力为阻力；70第七十页，共一百三十七页，2022年，8月28日东亚副热带急流中心受到的粘滞力大于0，小于0，等于0？仅考虑流体作直线运动71第七十一页，共一百三十七页，2022年，8月28日仅考虑流体作直线运动上层流体运动，图中处于静止状态的流体块受到的粘滞力大于0，小于0，等于0？72第七十二页，共一百三十七页，2022年，8月28日仅考虑流体作直线运动流体运动达到稳定状态（定常），图中流体块受到的粘滞力大于0，小于0，等于0？73第七十三页，共一百三十七页，2022年，8月28日3、欧拉方程理想流体（不考虑流体粘性），则纳维-斯托克斯方程：可以简化，相当于去掉方程中含有粘性的项。于是，方程简化为： 欧拉方程：理想流体的运动方程74第七十四页，共一百三十七页，2022年，8月28日例2-3-1已知流场u=ay，v=bx，w=0，其中a、b为常数，试根据不计质量力和流体粘性的运动方程，导出等压线方程。75第七十五页，共一百三十七页，2022年，8月28日

如流体静止时，即流体的速度和加速度的个别变化均为零，作用于流体的力应该达到平衡。此时，可得如下形式方程：

即所谓的静力方程。它表明了流体的粘性只与流体的运动状态有关，或者说流体的粘性只有在相对运动时才体现出来。4、静力方程76第七十六页，共一百三十七页，2022年，8月28日假设流体所受的质量力就是重力，静力方程可以变化为：上式表明：当流体静止时，作用于单位截面积流体柱的顶面、底面上的压力差，正好等于流体柱的重力；静力方程应用：77第七十七页，共一百三十七页，2022年，8月28日静力平衡条件下A点受到的压力？xz0hA在大尺度大气运动中，垂直运动速度很小（1/100m/s），大气科学中常用到静力方程(静力平衡)：静力平衡的应用：78第七十八页，共一百三十七页，2022年，8月28日已知为定义在某物质体上的标量，试证明：79第七十九页，共一百三十七页，2022年，8月28日

单位时间总外力的作功率内能动能单位质量的物质：第四节能量方程吸收或者释放的热量q研究对象流体中以界面包围的体积为的流体块单位时间总外力的作功率内能动能吸收或者释放的热量80第八十页，共一百三十七页，2022年，8月28日外界对系统所作的功率＋吸收或释放的热量（内能+动能）的变化率流体中以界面包围的体积为的流体块研究对象81第八十一页，共一百三十七页，2022年，8月28日方程变换总能量的变化项：热流量的变化率82第八十二页，共一百三十七页，2022年，8月28日表面力作功率项：83第八十三页，共一百三十七页，2022年，8月28日可以改写为：单位质量流体的能量方程，它是能量守恒定律在流体运动中的具体表现形式。流体块的能量守恒方程84第八十四页，共一百三十七页，2022年，8月28日动能方程根据流体的运动方程上式两端同乘速度矢量右端第二项展开后，则有：85第八十五页，共一百三十七页，2022年，8月28日单位质量流体微团的动能方程利用广义牛顿粘性假设86第八十六页，共一百三十七页，2022年，8月28日单位质量流体微团的动能方程物理意义：①②①质量力作功率②表面力作功率外力作功率引起的动能变化87第八十七页，共一百三十七页，2022年，8月28日E恒为正值粘性耗散项动能内能？88第八十八页，共一百三十七页，2022年，8月28日③③膨胀、收缩在压力作用下引起的能量转换项：膨胀收缩动能内能？动能内能？流体压缩性89第八十九页，共一百三十七页，2022年，8月28日热流量方程用能量方程减去动能方程反映内能变化率的热流量方程90第九十页，共一百三十七页，2022年，8月28日对于理想流体，即考虑无粘性，热流量方程简化为：“热力学第一定律”——能量转换和守恒定律在大气科学中所用的的形式。91第九十一页，共一百三十七页，2022年，8月28日例2-4-1设不可压缩流体平面无旋运动，试证明：在运动平面上任取周长线为S所围的单位厚度的流体块的动能可写为：92第九十二页，共一百三十七页，2022年，8月28日伯努利方程的适用条件：(1)无粘性流体(2)不可压缩流体(3)定常流动(4)质量力为有势力（保守力）伯努利方程理想不可压缩流体在重力作用下作定常运动时，流体的总机械能（动能、重力势能、压力能之和）沿着流线或迹线守恒。93第九十三页，共一百三十七页，2022年，8月28日对于理想流体，动能方程简化为：理想流体动能的变化，仅仅是由质量力和压力梯度力对流体微团作功造成的，而与热能不发生任何转换。故最终理想流体的动能方程可以写成：又因为94第九十四页，共一百三十七页，2022年，8月28日假设质量力是有势力，且质量力位势为，即满足：如考虑为一定常场，则有：95第九十五页，共一百三十七页，2022年，8月28日理想流体的动能方程假设质量力是有势力且为定常场96第九十六页，共一百三十七页，2022年，8月28日理想流体微团的动能方程：不可压缩定常97第九十七页，共一百三十七页，2022年，8月28日等式左端括号内部分的个别变化为零，即：理想不可压缩流体在重力作用下作定常运动时，流体的总机械能（动能、重力势能、压力能之和）沿着迹线守恒。98第九十八页，共一百三十七页，2022年，8月28日定常运动:流体运动的迹线和流线是重合于是沿流体运动的流线也有：伯努利方程99第九十九页，共一百三十七页，2022年，8月28日伯努利方程100第一百页，共一百三十七页，2022年，8月28日定常不可压缩各项点乘速度矢量101第一百零一页，共一百三十七页，2022年，8月28日例2-4-2理想不可压流体，所受质量力仅为重力的情况下作定常运动时，其中一流管如图所示，已知O点压力和速度均为零，讨论此时图中处于同一流线上A、B两点的流速VA、VB及压力PA、PB间的相对大小。O102第一百零二页，共一百三十七页，2022年，8月28日皮托管，又名“空速管”，“风速管”，英文是Pitottube。皮托管是确定气流速度的一种管状装置，由法国H.皮托发明而得名。下图是皮托管的结构示意图。它是由两个同轴细管组成,内管的开口在正前方,如图中A所示。外管的开口在管壁上,如图中B所示。两管分别与U型管的两臂相连,在U型管中盛有液体(如水银),构成了一个压强计,由U型管两臂的液面高度差h确定气体的流速。

VA＝0VB？皮托管示意图103第一百零三页，共一百三十七页，2022年，8月28日例：求定常条件下水从容器壁小孔中流出时的速率。解：水从小孔中流出时的流速可以根据伯努利方程求解。设ABC为一条流线。A和B分别是这条流线在水面和小孔处的两点,其中水面上点A和孔口处点B都与大气接触,所以那里的压强都等于大气压p0。容器的横截面比小孔的截面大得多,根据连续性方程,VA<<VB,故可以认为VA

=0。将以上条件代入上式,即可求得小孔处的流速,为h104第一百零四页，共一百三十七页，2022年，8月28日第五节简单情况下的N－S方程的准确解流体力学的基本方程组：运动方程连续方程考虑流体为均匀不可压缩（=常数），且粘性系数为常数（=常数）的情况下，方程组是闭合的。105第一百零五页，共一百三十七页，2022年，8月28日求解流体力学问题的一般方法，就是求解这样的闭合的方程组并使之适合应当的初始条件和边界条件。由于流体运动方程含有如平流加速度的非线性项，它是一个非线性方程组，在数学上求解这样一个非线方程组是难以做到的。求解方程组前，对初始条件和边界条件进行介绍。本节通过简单问题的求解－－了解基本方法

106第一百零六页，共一百三十七页，2022年，8月28日在初始时刻，基本方程组之解所应满足的既定条件，即在时，在定常流场的情况下，所有的流场参数均与时间无关，因而不存在初始条件的问题。（1）初始条件107第一百零七页，共一百三十七页，2022年，8月28日当流体流经固体壁时，必须满足不可穿透条件和无滑脱条件。 （2）边界条件而当固体壁运动时，则满足：当固体壁静止时，满足：固体壁边界流体与固体分界面上的条件108第一百零八页，共一百三十七页，2022年，8月28日（2）边界条件在自由表面上，两种流体在边界面上的法向速度应该相等，即：

另外，如果不考虑表面张力，两种流体质点在边界面上的法向应力应该相等，即:流体空气流体与流体分界面上的条件109第一百零九页，共一百三十七页，2022年，8月28日一、平面库埃托流动（PlaneCouetteFlow）h

h

Uu？zx考虑如下简单流动，设流体在两相距为2h的无界平行平板间，沿x轴作定常直线平面运动，此时满足：试确定流体的速度分布。上平板匀速运动下平板静止110第一百一十页，共一百三十七页，2022年，8月28日考虑了xoz平面的运动，则 假设流体是不可压缩的：可见，u仅仅是z的函数作定常直线平面运动111第一百一十一页，共一百三十七页，2022年，8月28日N-S方程简化为：积分沿x轴作定常直线平面运动质量力为重力：流体是不可压缩的。112第一百一十二页，共一百三十七页，2022年，8月28日方程第一式可以得到：积分上式可以得到：113第一百一十三页，共一百三十七页，2022年，8月28日设在x方向的压力分布均匀，即：已知边界条件：

最终可以得到：上式即给出了平面库埃托流动的流速分布，流速沿z轴呈线性分布。114第一百一十四页，共一百三十七页，2022年，8月28日二、平面普瓦瑟耶流动（PlanePoiseuilleFlow）h

h

zx在平面库埃托流动的基础上，假定沿x方向的压力梯度不为零，而上、下板处于静止状态。下平板静止上平板静止115第一百一十五页，共一百三十七页，2022年，8月28日此时，边界条件为： 即为平面Poiseuille流动的流速分布，流速沿z轴方向呈抛物线分布。将边界条件代入方程通解中，可以得到：116第一百一十六页，共一百三十七页，2022年，8月28日埃克曼流动

（EkmanFlow）1893～1896年，挪威海洋调查船“前进”号横越北冰洋时，F.南森观察到冰山不是顺风漂移，而是沿着风向右方20º～40º的方向移动。1905年，V.W.埃克曼研究了这种现象，得出了著名的埃克曼漂流理论。117第一百一十七页，共一百三十七页，2022年，8月28日三、埃克曼流动埃克曼螺线：风向随高度增大而向右旋转（北半球），风速随高度增加而增大，不同高度的风速矢量末端的连线为一螺线。118第一百一十八页，共一百三十七页，2022年，8月28日考虑粘性系数和密度均为常数的流体，在旋转角速度为的旋转坐标系中的运动，此时出现了地转偏向力的作用。而地转偏向力为：假设流体作平面运动，该平面绕轴转动，则流速表示为：119第一百一十九页，共一百三十七页，2022年，8月28日假设流体相对于旋转参考系无加速度，且无质量力作用，其运动方程（N-S方程）为：地转偏向力为：120第一百二十页，共一百三十七页，2022年，8月28日假设p与x，y，z无关地转偏向力与粘滞力相平衡121第一百二十一页，共一百三十七页，2022年，8月28日考虑u、v仅是z的函数，即满足：；