平板边界层内流速分布实验

平板边界层实验报告

流体力学实验 平板边界层实验报告 班级 姓名 实验日期 指导教师 北京航空航天大学流体力学研究所

流体力学实验 平板边界层实验报告 一、实验目的 测定平板边界层内的流速分布，并比较层流边界层及紊流边界层的速度分布的差别。 二、实验设备 本实验使用的是一个二维开路闭口低速风洞，在该风洞实验段中装有两块平板，以分别测量层流及紊流边界层的速度分布。为测量速度分布，在平板板面上安装有总压排管及静压管。这些测压管分别用橡皮管连接到多管压力计上，通过测量多管压力计液柱高度推算出速度来，具体原理见后。为测出实验段风速，在实验段侧壁上装有风速管，风速管的总压孔及静压孔也分别用橡皮管连接于多管压力计上，装备情况见图1。 图1 三、实验原理 当气流流过平板时由于粘性作用使紧贴平板表面处的流速为零，离开板面速度就逐渐增大，最后达到相当于无粘时的气流速度。对平板来说，就等于来流速度了。由于空气粘性很小，只要来流速度不是很小时，流速变化大的区域只局限在靠近板面很薄的一层气流

内，这一薄层气流通常叫作边界层。人为地规定，自板面起，沿着它的法线方向，至达到99%无粘时的速度处的距离，称为边界层厚度δ。 不可压流场中，每一点处的总压P 0，等于该点处的静压和动压 1 2 2ρv 之和。 p p v 021 2=+ ρ 则 v p p = -20() ρ （1） 因此只需测出边界层内各点处的静压p ，总压p 0，就可计算出各点的速度来。但考虑到垂直平板方向的静压梯度等于零（即??p y /=0），我们只需在平板表面开一静压孔，所测的静压就等于该点所在的平板法线方向上各点的静压。要测边界层内的速度分布就只要测出沿平板法线上各点的总压即可。 p i 0──为各测点的总压。 p i ──为各测点的静压。 v i ──为各测点的速度。 γ ──为多管压力计所使用的液体重度（公斤/米3）。 ?h i ──为各测点总压管与静压管的液柱高度差。 ρ ──为空气的密度，实验时可依据当时室温及大气压强由表查出。 φ ──为多管压力计的倾斜角。 根据（1）式，边界层内各测点处的速度为 v h i i = 2 ρ γφ?sin （2） 通常边界层内的速度分布用无量纲的形式表示为 v v f y i i 1=()δ y i 为各测点至板面的高度，δ 为边界层厚度，v 1为边界层外边界上的速度，对平板 来说即为来流速度。 v 1可通过风速管的静压管和总压管在多管压力计上的液柱高度差?h 1，由下式算出： v h 112 =ρ γφ?sin （3） 由（2）式和（3）式，可得 v v h h i i 1 1 =?? （4）

什么是边界层

什么是边界层?广义讲:在流体介质中,受边界相对运动以及热量和物质交换影响最明显的那一层流体。具体到大气边界层,是指受地球表面摩擦以及热过程和蒸发显著影响的大气层。大气边界层厚度,一般白天约为1.0km,夜间大约在0.2km左右,地表提供的物质和能量主要消耗和扩散在大气边界层内。大气边界层是地球-大气之间物质和能量交换的桥梁。全球变化的区域响应以及地表变化和人类活动对气候的影响均是通过大气边界层过程来实现的。 什么是湍流?英文湍流为“turbulence”,日文为“乱流”,湍流简单定义:流体微团进行的有别于一般宏观运动的不规则的随机运动,从宏观上看,它没有稳定的运动方向,但它能够象分子运动一样通过其随机运动过程有规律地传递物质和能量。从1915年由Taylor[1]提出大气中的湍流现象到1959年Priestley[2]提出自由对流大气湍流理论,可以说,到20世纪50年代以前经典的湍流理论基本上已经形成。以后,湍流理论基本上再没有出现大的突破。1905年Ekman[3]从地球流体力学角度提出了著称于世的Ekman螺线,在此基础上形成了行星边界层的概念,他的基本观点仍沿用至今。1961年,Blackadar[4]引入混合长假定,用数值模式成功地得到了中性时大气边界层具体的风矢端的螺旋图象。行星边界层的提出使人们认识到了大气边界层在大气中的特殊性和一些奇妙的规律。从20世纪50年代开始,由于农业、航空、大气污染和军事科学的需要,掀起了大气边界层研究的高潮。1954年, Monin和Obukhov[5]提出了具有划时代意义的Monin—Obukhov相似性理论,建立了近地层湍流统计量和平均量之间的联系。1982年,Dyer[6]等利用1976年澳大利亚国际湍流对比实验ITCE对其进行完善使得该理论有了极大的应用价值。1971年Wyngaard[7]提出了局地自由对流近似,补充了近地面层相似理论在局地自由对流时的空白。从20世纪70年代开始,随着大气探测技术和研究方法的发展,特别是雷达技术,飞机机载观测, 系留气球和小球探空观测以及卫星遥感和数值模拟等手段的出现,大气边界层的研究开始从近地层向整个边界层发展。简洁地概括,对大气边界层物理结构研究贡献最突出的是两大野外实验和一个数值实验,即澳大利亚实验的Wanggara和美国的Min-nesota实验以及Deardorff的大涡模拟实验。相似性理论是大气边界层气象学中最主要的分析和研究手段之一,在建立了比较成熟的用于描述大气近地面层的Monin—Obukhov相似性理论以后,人们开始寻求类似的全边界层的相似性理论。国际上,除Neuwstadt[8]、Shao[9]等做了大量工作外,我国胡隐樵等以野外实验验证了局地相似性 理论,并建立了各种局地相似性理论之间的关系。张强等还对局地相似性理论在非均匀下垫面近地面层的适应性做了一些研究。自1895年雷诺平均方程建立以来,该方程组的湍流闭合问题是至今未解决的一个跨两个世纪的科学难题。人们发展湍流闭合理论,以达到能够数值求解大气运动方程,实现对大气的数值模拟。闭合理论有一阶局地闭合理论即K闭合。1990年HoIt-sIag[12]在1972年理论框架的基础上,用大涡模拟资料对K理论做了负梯度输送的重大修正。为更精确地求解大气运动方程,也为了满足中小尺度模式,特别是大气边界层模式刻画边界层湍流通量和其它高阶矩量的目的,高阶湍流闭合技术也开始被模式要求。由于大气边界层研究是以野外探测实验为基础的实验性很强的科学,我国以往由于经济落后,无法得到第一手的实验资料,研究相对落后,与国外相比,总体上差距在20a左右,但我国学者在大气边界层的研究中也有其特殊贡献:1940年周培源先生[13]提出的湍流应力方程模式理论,被认为是湍流模式理论开始的标志,这一工作奠定了他在国际湍流研究领域的崇高地位。苏从先等在上世纪50年代给出的近地面层通量廓线与当时国外同类研究同步,被国外学者称为“苏氏定律”,在上世纪80年代苏从先等首次发现了干旱区边界层的绿洲“冷岛效应”结构。上世纪70年代周秀骥[16]提出的湍流分子动力学理论也很有独特的见解。1981年周 明煜[17]提出的大气边界层湍流场团块结构是对湍流结构的新认识。上世纪80~90年代赵鸣[18]对边界层顶抽吸作用的研究是对Charney—Eiassen公式的很好发展。在20世纪90年代的“黑河实验”中,胡隐樵等和张强[19]首次发现了邻近绿洲的荒漠大气逆湿,并总结提出了绿洲与荒漠相互作用下热力内边界层的特征等等。国内外有关大气边界层和大气湍流的专著

matlab求解平板边界层问题

《粘性流体力学》程序 平板边界层问题求解

1.1编程思路 平面边界层问题可以归结为在已知边界层条件下解一个高阶微分方程,即解0 f。Matlab提供了解微分方程的方法，运用换+ff '''= 5.0 '' 元法将高阶微分方程降阶，然后运用“ode45”函数进行求解。函数其难点在于如何将边界条件中1 η运用好，由四阶龙格-库塔方 →f ,→ ' ∞ 法知其核心是换元试算匹配，故在运用函数时通过二分法实现 η是可行的。 ∞ →f ,→ ' 1 1.2m函数 function dy = rigid(x,y) dy = zeros(3,1); dy(1) = y(2); dy(2) = y(3); dy(3) = -0.5*y(1)*y(3); %main程序 [X, Y] = ode45('rigid',[0 5],[0 0 0]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') %二分法试算f’’的初始值以满足f’趋向无穷时的边界条件，图像上可以清晰看出f’无穷时的结果 >> [X, Y] = ode45('rigid',[0 5],[0 0 1]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') >> [X, Y] = ode45('rigid',[0 5],[0 0 0.5]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') >> [X, Y] = ode45('rigid',[0 5],[0 0 0.25]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') >> [X, Y] = ode45('rigid',[0 5],[0 0 0.375]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') >> grid on >> [X, Y] = ode45('rigid',[0 5],[0 0 0.3125]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') >> grid on >> [X, Y] = ode45('rigid',[0 5],[0 0 0.34375]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') grid on >> [X, Y] = ode45('rigid',[0 5],[0 0 0.328125]); plot(X, Y(:,1),'-',X, Y(:,2),'*',X, Y(:,3),'+') grid on

大气探测学复习思考题版

大气探测学复习思考题（2011版）一、写出下列云状的国际简写或由国际简写写出云状学名 浓积云Cu cong 碎积云Fc 淡积云Cu hum 秃积雨云Cb calv 鬃积雨云Cb cap 荚状层积云Sc lent 堡状层积云Sc cast 透光层积云Sc tra 积云性层积云Sc cug 蔽光层积云Sc op 层云St 碎层云Fs 雨层云Ns 碎雨云Fn

透光高层云As tra 蔽光高层云As op 透光高积云Ac tra 蔽光高积云Ac op 堡状高积云Ac cast 荚状高积云Ac lent 积云性高积云Ac cug 絮状高积云Ac flo 毛卷云Ci fil 密卷云Ci dens 伪卷云Ci not 钩卷云Ci unc 匀卷层云Cs nebu 毛卷层云Cs fil 卷积云Cc 二、解释名词 大气科学、大气探测、气象资料的代表性、气象资料的准确性、气象资料的比

较性、云、、云量、天气现象、气象能见度、气象光学距离、气温、摄氏温标、华氏温标、热电现象、热滞系数、百叶箱、湿度、露点温度、盖﹒吕萨克尺度、气压、本站气压订正、海平面气压订正、风、阵风、降水量、蒸发量、积雪、太阳常数、直接辐射、雾、环日辐射、散射辐射、全辐射、净辐射、日照时数、高空测风、单经纬仪定点测风、双经纬仪基线测风、一次雷达、二次雷达、测风雷达的测角原理、等信号强度法、自动气象站、遥感、主动式大气遥感探测、被动式大气遥感探测、激光雷达、声雷达、可见光探测、红外辐射探测、微波探测、大气边界层探测、气象塔、对比视感阈 三、简述或论述下列各题 1．为什么要提出气象观测资料的“三性”？ 2．什么是观测资料的测站代表性和区域代表性？ 3．怎样来衡量观测资料的代表性和准确性？它们之间有何关系？怎样保证比较性？ 4.淡积云、浓积云、秃积雨云、鬃积雨云，它们之间的区别界限是什么？ 5.碎积云、碎层云、碎雨云，它们之间在外形及成因上有何不同？ 6.卷层云和高层云、高层云和雨层云、雨层云和层云，各有何异同之处？ 7.卷积云和高积云、高积云和层积云，各有何异同之处？

流体力学平板边界层内的流速分布实验报告电子版

平板边界层内的流速分布实验 实验日期 2011-5-21 小组成员：李超，郭静文（93班）等 报告人 周楠 能动95 09031125 实验目的 1） 测量离平板前缘任意截面边界层内的速度分布； 2） 根据速度分布确定边界层厚度； 3） 了解风洞结构及测量仪器。 仪器设备 吸入式风洞、大气压强计、温度计、微压计、U 型测压管、平板模型、总压探针及三维坐标架。 其中仪器的重要参数包括： (1)吸入式低速风洞P max =P a , 工作截面尺寸300mm ×300mm; (2)风洞的气体流速u max <25m/s, M<0.3,所以风洞内气体流动可以看成二维不可压缩流动即ρ=ρa (3)平板尺寸325mm ×200mm (4)总压探针头部直径：d=0.9mm 实验原理 1 流体在大雷诺数下绕物体流动时，由于流体粘性的作用，与物体表面接触的流体速度为零，然后沿法向很快增至主流速度，这层贴近物体表面，沿着法向有很大速度梯度的流动薄层，称为边界层； 2 在边界层内，速度梯度很大，不能忽略流体的粘性，因此流动作实际流动u x 和p o 都在变化且u x

大气探测Word版

第一章 1.什么是大气探测？可以划分为哪几部分? 2.简述大气探测的目的是什么? 3.地面气象观测场的要求有哪些？ 4.简述大气探测有哪几种方法? 5.大气探测仪器的性能包括哪几个? 6.如何保证大气探测资料的代表性和可比性? 7.气象探测环境? *大气探测 ．．．．，是指利用科技手段对大气层和近地层的各种物理过程、化学过程、生物过程等进行系统的观察和测量。是对表征大气状况的气象要素、天气现象及其变化过程进行系统的、连续的观察和测定。大气探测资料精确性直接影响到气象预报、气候预测和气候变化研究等的准确性。大气探测可划分为地面气象观测，高空探测，特种探测，遥感探测 *．大气探测的科技手段 ．．．．．．．．．：机械方法：风向标等。光学方法：经纬仪光谱仪等。热学方法：蒸发器等。电磁方法：雷达卫星信息传输等 * 大气探测的空间范围， ．．．．．．．．．．从海平面往上到几百公里甚至一千多公里 范围不断扩大！ * 大气边界层， ．．．．．．大气与地面之间充分湍流化的气层，厚度约为1km，大气要素变化剧烈并随时间和地表条件变化，是地面气象观测和大气边界层探测的主要大气层，该层对局地天气、气候产生重要影响。

近地面层（常通量层，表面层），直接与地表接触、受地面强烈影响， 地面气象观测在此高度内。 *影响大气边界层探测的物理过程，．．．．．．．．．．．．．．．辐射传输过程：短波辐射，长波 辐射。热力传输过程：显热，潜热。动力作用：平流，垂直运动，局 地环流。湍流运动：无规则，分类（机械、热力），特征，尺度谱， 串级传递能量过程，湍流输送过程。 *．大气探测学研究目的 ．．．．．．．．．，．1.大气探测是从事大气科学教学、科研的基 础。为天气、气候诊断分析、预报及环境保护部门、国家及全球气象 资料网络系统等提供大气观测资料。2.随着科学技术的发展，大气探 测的要素量和空间范围越来越大。大气探测分为近地面层大气探测、 高空大气层探测和专业性大气探测。3.. 是推动大气科学进步的重要 力量大气科学重大理论的提出和完善得益于大气探测（虽然数值模拟 在其中也发挥着重要作用）。新的探测原理和观测平台促进大气科学 的前进（如GPS水汽探测、卫星平台） *大气探测 ．．．，业务需求，社会需求，科研需求。．．．．同时满足 ．．．．3．种需求 *地面气象观测场， ．．．．．．．．观测场四周空旷平坦，所取得的资料应具有较好 的代表性；经纬度（精确到分）和海拔高度（精确到0.1m）刻在石 碑上，埋设在场内；观测场为25mX25m的平整场地，保持均匀草坪， 草高不超过20cm，不准种植作物；设 1.2m高稀疏围栏，内设 0.3m-0.5m宽小路，且只准在小路上行走，小路下建线缆沟或埋设 线缆管。 *气象探测环境 ．．．．．．，是指为避开各种干扰保护气象探测设施准确获得气

平板边界层速度剖面的测定讲义2

2009年04月20~22日平板附面层速度剖面与厚度的测定 一、实验目的： 1．熟悉附面层速度分布和厚度的测量方法。 2．具体测定平板附面层层流与湍流附面层的速度分布及其厚度。 3．把实验结果与理论计算结果进行比较，分析其差异产生的原因。 二、实验原理： 粘性匀质不可压缩流体，测量边界层内的速度，仍利用风速管（皮托管）测风速的原理，即测出某点的总压P0和静压P后再换算成该点的速度，因为边界层很薄，其厚度往往只有几mm到十几mm，因而只能用极细的探针去探测边界层内的压力。 由于在边界层内部满足?(P)/?(Y)=0，即静压P沿着平板的法线方向不变，因此，可以用壁面上的静压P来表示边界层内法线上所有不同高度的静压。于是，本实验将一根微总压管装在一标架上，使微总压管以很小的间距上下移动，测出不同高度处的总压P0(y)后，即可算出法线上离壁面y处的速度。 实验时，把总压管由壁面逐步往上移动，则测出的总压越来越大。当移动到某一高度以后，再继续往上移动几个间距，这时所测到的总压已不再随高度的变化而变化。记录下数据，经软件分析后可得速度边界层厚度和速度剖面，并与理论曲线对照。 理论分析中总是假定从平板（或物体）的前缘（或驻点）就开始形成层流或湍流边界层。实际上绕流体的运动常常是组合边界层问题，即在物体的前部分首先形成层流边界层，在它的后部分形成湍流边界层，在它们之间还有一个过渡段。 过渡段从层流的失稳点（层流不稳定点）开始直到流动成为完全湍流之点（湍流过渡点）结束。性质介于两者之间。 为了读出压力的微小变化,本实验采用压力传感器,采用总压和静压之差，将其采集的压力信号转换成电信号，再通过放大器进行信号放大后，输入A/D转换器，由计算机直接计算出速度值。 由于速度剖面是以无量纲形式画成的，因此，不需要计算一点的速度，只要计算出速度的相对值就可以了。计算各高度上的u y/v和y/δ的值，以y/δ为纵

实验 两平行平板间的缝隙流动试验

实验 两平行平板间的缝隙流动试验 一、实验目的 1. 两平板之间（平板之间没有相对运动）充满了不可压缩流体，测量其流量，温度，黏度，密度，及流速及其雷诺系数，分析流体流动的动态特性以及不可压缩流体以均匀的速度U 沿二元平板作恒定流动时边界层的厚度和壁面切应力的分步规律。 2. 通过紊流对平板的作用力和平板对紊流的反作用力验证不可压缩流体定常流动的动量方程。 3.不同的边界情况下平行平板缝隙流，写出截面上每个点的切应力分布，平板上的摩察应力，摩察系数。 二、实验装置 实验装置如图2.1所示。 雷诺试验装置主要由稳压溢流水槽、试验导管和转子流量计等部分组成，如图1所示。自来水不断注人并充满稳压溢流水槽。稳压溢流水槽的水流经试验导管和流量计，最后排入下水道。稳压溢流水槽的溢流水，也直接排入下水道。 图2.1 动量方程实验装置简图 三、实验原理 经许多研究者实验证明：流体流动存在两种截然不同的型态，主要决定因素为流体的 密度和粘度、流体流动的速度，以及设备的几何尺寸（在圆形导管中为导管直径）。 将这些因素整理归纳为一个无因次数群，称该无因次数群为雷诺准数（或雷诺数），即 ()1 u d R e μ ρ= 式中d 一导管直径，m ρ一流体密度，kg ·m -3 ； μ一流体粘度，Pa · s ； u 一流体流速，m · s -1 ； 大量实验测得：当雷诺准数小于某一下临界值时，流体流动型态恒为层流；当雷诺数 大于某一上临界值时，流体流型恒为湍流。在上临界值与下临界值之间，则为不稳定的过 渡区域。对于圆形导管，下临界雷诺数为2000，上临界雷诺数为10000。一般情况下，上临 界雷诺数为400O 时，即可形成湍流。 应当指出，层流与湍流之间并非是突然的转变，而是两者之间相隔一个不稳定过渡区 域，因此，临界雷诺数测定值和流型的转变，在一定程度上受一些不稳定的其他因素的影 响。 四、实验步骤及注意事项

流体力学——平板边界层编程

对于本次编程编程作业，小组运用matlab 和c++两种程序对平板边界层问题和绕过楔形体边界层流动问题进行分析研究。以下是运用matlab 解决问题的过程。 一、 平板边界层问题 该问题可以归结为在已知边界层条件下解一个高阶微分方程,即解0''5.0'''=+ff f 。Matlab 提供了解微分方程的方法，运用换元法将高阶微分方程降阶，然后运用“ode45”函数进行求解。函数其难点在于如何将边界条件中1',→∞→f η运用好，由四阶龙格-库塔方法知其核心是换元试算匹配，故在运用函数时通过二分法实现1',→∞→f η是可行的。程序如下： 第一问m 函数 function dy = rigid(t,y) dy = zeros(3,1); dy(1) = y(2); dy(2) = y(3); dy(3) = -0.5*y(1)*y(3); %第一问main 程序 [T,Y] = ode45('rigid',[0 5],[0 0 0]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') %二分法试算f ’’的初始值以满足f ’趋向无穷时的边界条件，图像上可以清晰看出f ’无穷时的结果 >> [T,Y] = ode45('rigid',[0 5],[0 0 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> [T,Y] = ode45('rigid',[0 5],[0 0 0.5]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> [T,Y] = ode45('rigid',[0 5],[0 0 0.25]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> [T,Y] = ode45('rigid',[0 5],[0 0 0.375]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> grid on >> [T,Y] = ode45('rigid',[0 5],[0 0 0.3125]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> grid on >> [T,Y] = ode45('rigid',[0 5],[0 0 0.34375]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on >> [T,Y] = ode45('rigid',[0 5],[0 0 0.328125]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on >> [T,Y] = ode45('rigid',[0 10],[0 0 0.328125]);%当f ’’为0.328125时，逼近结果已经很好，在0到5的变化范围内已经非常接近精确解 plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on >> [T,Y] = ode45('rigid',[0 5],[0 0 0.335975]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on

平板边界层实验

平板边界层实验（一） （一）实验目的 1．测定平板边界层内的流速分布，从而确定流速分布指数规律、边界层名义厚度δ、位移厚度1δ、动量厚度2δ、能量厚度3δ。 2．掌握毕托管和测压计的测速原理和量测技能。 （二）DQS 系列空气动力学多功能实验装置： 该装置相当于小型风洞，为组装式结构。由主机和多种易更换实验段组成，流量可以控制。风机提供气流，在压出段设有流量调节阀门，气流通过风道进稳压箱流速减慢进入阻尼网，阻尼网由二层细密钢丝网构成，可将流体较大尺度的旋涡破碎，使气流均匀地进入收缩段，经过收缩段可将收缩段进口的速度不均匀度缩小n 2倍，n 为收缩比，本收缩段的收缩比较大。收缩曲线应用波兰人维托辛斯基曲线。收缩段出口接各种实验段，实验排放的气流由实验台面的孔口进吸音箱回到风机入口，如图1所示。 多管测压计，设有可改变角度的测压排管及调平设置，当测某点压强时取与大气连通的测压管与该点测压管的读数差，即为测点的压强水头，如图2所示。 1.稳压箱 1.测压管 2.收缩段 2.角度盘 3.风道 3.支架 4.调节阀门 4.联通管 5.通风机 5.输液管 6.吸音箱 6.酒精库 7.阻尼网 7.通气管 图 1 图 2 （三）实验段简图 稳压箱内的气流经过阻尼网及收缩段均匀进入实验段，在实验段轴心位置安装一块一面光滑一面粗糙的平板，平板可沿轴线滑动，在实验段的出口装有精致的鸭咀形毕托管，其

头部厚度仅有0.3㎜,并配有千分卡尺，灯光显示设置和多管测压计，见图1-1。 （三）实验原理及计算式 1．平板紊流边界层的流速分布 实际流体因存在粘性，紧贴壁面的流体将粘附于固体表面，其相对速度为零，沿壁面法向随着与壁面距离的增加，流体的速度逐渐增大，当距离为δ时，其速度达到未受扰动的主流流速∞u ，这个厚度为δ的薄层称为边界层，通常规定从壁面到∞=u u x 99.0处的距离作为边界层的厚度。 边界层的厚度沿平板长度方向是顺流渐增的，在平板迎流的前段是层流边界层，如果平板足够长，则边界层可以过渡到紊流，判别过渡位置的特征值是雷诺数x Re ，如图1-2所示。 若量测断面坐标为x ，则该断面x Re 为 ν x u x ∞= Re (本装置用0u 代表∞u ) ( 1－1 ) 其中ν为空气运动粘滞系数，α ρμν= μ为动力粘滞系数，αρ为空气密度。 n T T ??? ? ??=00μμ 15.2880=T °Ｋ时，250/·10789.1m S N -?=μ K ℃t T 15.273/+= 90°K ＜T ＜300°K 时，n 取8/9，代入 29 85/· 15.28810789.1m S N K T ?? ? ????=-μ 紊流边界层内的流速分布用指数律表示为 图 1 - 1 n x y u u 10?? ? ??=δ ( 1－2 ) 式( 1－1 ) 、( 1－2 )中

边缘检测实验报告

图像边缘提取实验报告 一、实验目的 通过课堂的学习，已经对图像分割的相关理论知识已经有了全面的了解，知道了许多图像分割的算法及算子，了解到不同的算子算法有着不同的优缺点，为了更好更直观地对图像分割进行深入理解，达到理论联系实际的目的，特制定如下的实验。 二、实验原理 检测图像边缘信息，可以把图像看做曲面，边缘就是图像的变化最剧烈的位置。这里所讲的边缘信息包含两个方面：一是边缘的具体位置，即像素的坐标；而是边缘的方向。微分算子有两个重要性质：定域性(或局部性)、敏感性(或无界性)。敏感性就是说，它对局部的函数值变化很敏感，但是因其对变化过于敏感又有了天然的缺陷——不能抵抗噪声。局部性意思是指，每一点的导数只与函数在该点邻近的信息有关。 主要有两大类基于微分算子的边缘检测技术：一阶微分算子边缘检测与二阶微分算子边缘检测。这些检测技术采用以下的基本步骤： (1) 将相应的微分算子简化为离散的差分格式，进而简化为模板(记为 T)。 (2) 利用模板对图像f(m,n)进行运算，获得模板作用后的结果Tf(m,n)。 (3) 提出阈值h,在采用一阶微分算子情形记录下高于某个阈值h 的位置 坐标 }),(|),{(h n m Tf n m S h ≥= (而采用二阶微分算子情形，一般是对某个阈值0>ε确立

}),(|),{(ε≥=n m Tf n m S h ) (4) 对集合h S 进行整理，同时调整阈值h 。 Roberts 算子 Roberts 算子是一种利用局部差分算子寻找边缘的算子，两个模板分别为 ??????-=1001x R ?? ? ???-=0110y R 则，),(j i f R x =)1,1(),(++-j i f j i f ),(j i f R y =)1,(),1(+-+j i f j i f 算法的步骤为： (1) 首先用两个模板分别对图像作用得到f R x 和f R y ； (2) 对2 2 ),(y x R R j i Tf +=，进行阈值判决，若),(j i Tf 大于阈值则相应的点 位于便于边缘处。 对于阈值选取的说明：由于微分算子的检测性能受阈值的影响较大，为此，针对具体图像我们采用以下阈值的选取方法，对处理后的图像统计大于某一阈值的点，对这些数据求平均值，以下每个程序均采用此方法，不再做说明。 Sobel 算子 Sobel 算子采用中心差分，但对中间水平线和垂直线上的四个邻近点赋予略高的权重。两个模板分别如下： ????? ??---=101202101x S ???? ? ??---=121000121 y S

平板边界层速度分布测量

平板边界层速度分布测量实验指导书 实验目的： 通过零迎角平板流动的流速测量，获取流速沿物面法向分布。 学习总压管测速。 实验装置和仪器： （1）风洞：回流开口小型风洞，试验段见 右图，矩形有机玻璃管道中夹放一 金属板，来流沿管道被该板分开， 从出口流出。出口截面的静压为大气 压。 （2）偏平总压探针头：偏平总压探针头顶可 在出口截面内水平移动，移动量由微分尺控制。 （3）酒精斜管压力计：斜角θ=30o，系数K=1.0, 一头通大气，另一头接总压探头。 实验原理： 测量原理，就是伯努利定理：不计重力，气流的动压和静压之和为总压。 设总压为P 0，则 )(])()([2 1)(220y P y v y u y P ++=ρ （1） y 为探头中心距平板的距离，u 、v 分别为平行于平板的流速和平板法向的流速， p 为当地静压，ρ为气流的密度。 因为 a P y P ≡)( ， u v << 由（1）可得 ρ])([2)(0a P y P y u -= （2） 实验步骤 : 图 风洞试验段示意图

（1）实验室大气参数读取和记录； （2）探头零位确定； （3）压力计底座水平调解，测压管液面零刻度调节； （4）风洞开车； （5）调节好探头距平板的距离y ，从压力计读取并记录相应的压力值Po-Pa 实验要求: 测压时，每移动探头至新位置，应等待几秒钟，在压力平衡后再读取数据。测量中，观察随探头离开平板距离的增大，压力的变化趋势。 实验报告要求： （1）实验参数：大气压P a （毫米汞柱） ，大气温度t (?C ) ， 大气密度 ) (15.273)(464.0C t mmHg P o a +?=ρ （公斤/米3） 。 测量为之举平板前缘的距离X ； （2）测压原始数据，及由（2）是换算成流速，给出曲线y y u -)(； （3）找出不随距离y 而变的速度值，记为U 1，并找出满足u(y)= U 1的最小的y 值作为δ，给出曲线δ//)(1y U y u -。并给出雷诺数Re=ρU 1X/μ。

基于小波变换的风廓线雷达大气边界层高度探测

基于小波变换的风廓线雷达大气边界层高度探测 戈书睿1)艾未华1)魏浩1)徐芬2) （1 解放军理工大学气象海洋学院，南京，211101） （2江苏省气象科学研究所，南京，211101） 大气边界层是地表与大气相互接触的一层大气，响应地面作用的时间尺度为一小时或者更短，大气边界层的厚度虽然只有几公里，但它是与人类关系最为密切的大气层。大气边界层的高度是污染物扩散模式、气候模式、大气模式的一个重要输入参数，边界层高度的变化对数值预报中的物理过程，天气预报的诊断分析，城市污染物的监控也有相当重要的作用。然而，边界层高度的连续监测缺乏有力的手段。风廓线雷达凭借其高时间分辨率和空间分辨率，加上其能够连续探测等优点，可以成为连续监测边界层高度的有效工具。利用协方差小波变换对风廓线雷达距离订正后信噪比数据进行分析，来确定边界层的高度，并与无线探空仪确定的边界层高度进行比较。采用不同的小波基函数对典型的风廓线雷达距离订正后的信噪比数据廓线进行变换，比较不同变化尺度下几种常用小波基函数变换的稳定性，得到最优的反演边界层高度的小波基函数，为实际应用提供参考。结果表明，协方差小波变换法与探空结果具有较好的一致性。 关键词：协方差小波变换；风廓线雷达；边界层高度；探测 参考文献： 1 Stull R B. An introduction to boundary layer meteorology[M]. Springer, 1988. 2 Ottersten H. Atmospheric structure and radar backscattering in clear air[J]. Radio Science, 1969, 4(12): 1179-1193. 3 Cohn S A, Angevine W M. Boundary layer height and entrainment zone thickness measured by lidars and wind-profiling radars[J]. Journal of Applied Meteorology, 2000, 39(8): 1233-1247. 4 Lippmann M. Health effects of tropospheric ozone[J]. Environmental science & technology, 1991, 25(12): 1954-1962. 5 Bell M L, McDermott A, Zeger S L, et al. Ozone and short-term mortality in 95 US urban communities, 1987-2000[J]. Jama, 2004, 292(19): 2372-2378. 6 Kane E S, Valentine D W, Schuur E A G, et al. Soil carbon stabilization along climate and stand productivity gradients in black spruce forests of interior Alaska[J]. Canadian Journal of Forest Research, 2005, 35(9): 2118-2129. 7 Cha J S, Choi J C, Ko J H, et al. The low-temperature SCR of NO over rice straw and sewage sludge derived char[J]. Chemical Engineering Journal, 2010, 156(2): 321-327. 8 Angevine W M, White A B, Avery S K. Boundary-layer depth and entrainment zone characterization with a boundary-layer profiler[J]. Boundary-Layer Meteorology, 1994, 68(4): 375-385. 9 A. B. White, C. J. Senff, and R. M. Banta, 1999: A Comparison of Mixing Depths Observed by Ground-Based Wind Profilers and an Airborne Lidar.J. Atmos. Oceanic Technol.,16, 584–590. 10 Heo B, Jacoby-Koaly S, Kim K, et al. Use of the Doppler Spectral Width to Improve the Estimation of the Convective Boundary Layer Height from UHF Wind Profiler Observations[J]. Journal of Atmospheric & Oceanic Technology, 2003, 20(3):408-424. 11 Bianco L, Wilczak J M. Convective Boundary Layer Depth: Improved Measurement by Doppler Radar Wind Profiler Using Fuzzy Logic Methods[J]. Journal of Atmospheric & Oceanic Technology, 2002, 19(11):1745-1758.

平板边界层测量

二、平板边界层速度剖面测量 1实验目的： 了解平板边界层特性，学习测量平板边界层速度剖面的方法。 在离平板前缘不同位置处，测量平板边界层内速度分布，确定边界层厚度，并和理论值进行比较。 2实验装置： 图 1 实验装置示意图 图 2 平板边界层测量原理 （1）平板：在三维小风洞中安装一块宽240毫米、长750毫米的尖前缘平板。平板表面光滑，零攻角安装。沿平板中线有若干静压孔（见上图）。 （2）总压管：头部直径1毫米的总压管，用于测量边界层内总压分布。总压管安装在坐标架上，总压管前端与静压孔齐平，小孔对准气流轴线且与平板平行。 （3）坐标架：安装在风洞上方，用于调节总压管位置。 （4）压力扫描测试仪：用于测量压差。使用时需注意仪表初始读数，以便对测量值进行修正。仪表拨盘位置与平板上测点相对应。 3实验步骤：

（1）安装好平板，并使其表面与风洞轴线平行。安装好总压管，使其对准气流方向并与平板平行。 （2）将总压管、静压孔分别与压力扫描测试仪相连。 （3）记录当天大气压和温度和仪表初读数。 （4）将总压管降到刚好与平板表面接触（必须反复调整总压管数次，以求找到最佳位置）。 这时总压管中心离平板表面的高度为y1=h/2 (h为总压管，外径=1mm)，此时坐标架的位置高度应为0.5毫米。 （5）启动风洞，调整到设定风速（变频器频率植）。记录仪表读数。 （6）上下移动坐标架，改变总压管位置，重复测量边界层内压力分布和总压管高度。由于总压管较细、管道较长，压力平衡需要一定时间。实验中要等到压力平衡后再读数。 总压管上下移动步长为1mm。 （7）重复步骤（6），直到压力计读数不再随总压管位置不同而改变为止。这时表明总压管已经到达边界层外面。由于接近边界层外边界时速度变化很小，所以必须再要往上移动总压管若干次，确认总压管已经到达边界层外部。 （8）改变总压管水平位置，同时转动压力扫描测试仪拨盘，使指针指向对应静压孔位置。 重复上述步骤，测量3-5个边界层速度剖面。 （9）风洞停车。 （10）整理实验数据，按照要求完成实验报告。 注：将总压管处于边界层外均匀区时测得的P0－P对应的气流速度作为来流速度。 4数据处理： （1）计算边界层内速度分布，速度边界层厚度。 由于边界层速度剖面是以无量纲形式画出的，不需要计算出每一点的速度，只要计算出相对速度就可以了。设y处的速度u y为 边界层外缘的速度U为 其中p0为总压管测得的压力，p为静压孔测得的压力，（p0一p）y是边界层内测得的读数。（p0一p）表示气流均匀区测得的读数，ρ为气体密度。相对速度为 当u y / U =0.99 时, 总压管的高度就是该处平板边界层的厚度。 (2) 计算边界层位移厚度及动量厚度.

北京城区大气边界层的激光雷达观测_杨辉

第31卷第2期2005年3月光学技术 OPT ICAL T ECHN IQU E Vol.31No.2 M ar.2005 文章编号:1002-1582(2005)02-0221-03 北京城区大气边界层的激光雷达观测X 杨辉1,刘文清1,陆亦怀1,谢品华1,徐亮1,赵雪松1,虞统2,于建华2 (1.中国科学院安徽光学精密机械研究所,安徽合肥230031; 2.北京市环境保护监测中心,北京100044) 摘要:激光雷达是边界层大气气溶胶和云的一个高效探测工具。利用M P L激光雷达测量出的数据反演出了测站地域上空大气边界层气溶胶消光系数的垂直分布,利用气溶胶消光系数的垂直分布可决定测站上空的大气边界层高度。测量数据表明:大气边界层内气溶胶的含量较为稳定,有明显的气溶胶多层结构。 关键词:激光雷达;气溶胶;大气边界层;边界层高度 中图分类号:T N958.98;X513文献标识码:A PBL observations by lidar at Peking YAN G Hui1,LIU Wen-qin g1,LU Y-i huai1,XIE Ping-hua1,XU Lian g1 ZH AO Xue-shong1,YU Tong2,YU Jian-hua2 (1.Anhui Institute of Optics and Fine M echanics,Chinese A cademy of Sciences,Hefei230031,China) (2.Beijing M unicipal Env ironmental M onitor ing Center,Beijing100044,China) Abstract:T he ver tical di stribution of PBL aerosol extinction coefficients over the obser vation site is show n,and t he heights of PBL are discussed.T he data indicates that the aeroso l o ver t he measurement site is trapped in PBL and is r at her stable,and the mult-i layer str ucture of the aerosol distribution is obvious. Key words:lidar;aerosol;PBL(planetary boundary layer);PBL height 1引言 大气边界层(PBL,planetary boundary layer)通常是指大气的最低部分受地面的直接影响,并与地面有直接作用的气层,是地球表面与自由大气间进行物质、能量、热量和水汽交换必经的气层。它主要是通过摩擦阻力、蒸发和蒸腾、热量输送、污染物排放以及影响气流变化的地形等与地面的作用。大气边界层是与人类关系最为密切的一层,是由人类活动和各项生态环境构成的主要气层。由于热力作用导致的强烈的日变化是大气边界层的一个重要特征。边界层高度随地表特征、季节和天气背景的不同而不同,每天可在几十米至几千米内变化[5,6]。 气溶胶是大气边界层的垂直结构和时空变化观测的天然载体。气溶胶是指液态或固态微粒均匀地散布在气体中形成的相对稳定的悬浮体系,它在地气系统的收支平衡中扮演着十分重要的角色。气溶胶对全球环境的影响也日趋严重,如气相物质经过光化学反应后急剧向固态微粒转化产生光化学烟雾,并进而生成城市气溶胶。这些光化学过程大大加剧了大气污染的程度,会直接损害人们的身体健康,影响正常的生产和工作秩序,会破坏大自然的固有生态平衡[1,7]。大气边界层内的气溶胶对地面气候、大气温度和气压场、风场和垂直运动场等都有明显的影响[8]。 激光雷达已逐渐成为大气边界层探测的主要工具。如拉曼(Ram an)激光雷达可探测边界层中水汽含量的时空分布;差分吸收(DIAL)激光雷达可探测大气边界层中污染气体,如NO2,SO2,O3等含量的时空分布;多普勒激光雷达可探测大气边界层内风场的时空分布;Mie散射激光雷达可连续地探测大气边界层中气溶胶粒子的光学特性以及气溶胶粒子和大气边界层高度的时空分布。 北京奥运会即将在2008年举行,利用激光雷达对其城区大气边界层中的气溶胶粒子进行监测是一项具有极高的科学意义和实用价值的工作。本文给出了影响北京市大气质量的主要因素以及在大气环境容量研究中的激光雷达对大气边界层的观测结果。 2实验中的激光雷达系统 221 X收稿日期:2004-06-07E-mail:hyang@https://www.wendangku.net/doc/6c8537973.html, 作者简介:杨辉(1974-),男,四川省德阳市人,中国科学院安徽光学精密机械研究所博士研究生,主要从事激光雷达大气探测、光电子技术等方面的研究。