计算流体力学---报告
工程流体力学实验
1.40
0
3.500
7.500
3.500
11.000
6.100
19
1.40
0
3.500
6.208
3.500
9.708
0.000
1.292
毕托管测速计算表
编号
6
8
12
14
16
18
备注
测速管读数
44.75
23.9
22.5
14.7
12
11
测压管读数
14.7
15.05
11.85
6.9
9.6
3.5
点流速u（cm/s）
三、使用仪器、材料
自循环供水器、恒压水箱、溢流板、稳水孔板、可控硅无级调速器、实验管道、流量调节阀、接水阀、接水盒、回水管测压计。
四、实验步骤
1、熟悉实验仪器，分清普通测压管和测速管及两者功能上的区别。
2、打开电源，启动供水系统，水箱供水至溢流，排净实验管道内的空气后关闭流量调节阀。检查所有的测压管液面是否齐平，若不平需查明原因并排除气体。
8、在均匀流断面上，推求测速管处的流速，将测试与计算成果列于表中。
水箱面高程 =47.60cm直径
实验装置图：
五、实验过程原始记录(数据、图表、计算等)
测点液面读数于断面能量转换的测算表单位：cm
测点
管径d
位置水头Z
压强水头p/γ
流速水头
测压管水头z+ p/γ
总水头H
测压管水头差△(z+ p/γ)
水头损失h=-=
如果自由表面压强p0与当地大气压pa压强相等时，液体内任一点相对压强可表示为：
式中：h为液体自由表面下任一点液体深度。
cfd仿真分析报告
CFD仿真分析报告1. 研究背景和目的本报告旨在使用计算流体力学（Computational Fluid Dynamics, CFD）方法对某个特定问题进行仿真分析，并提供详细的模型设置、计算方法和结果分析。
2. 模型设置2.1 几何模型在本次仿真分析中，我们选取了某个具体的几何模型进行研究。
该几何模型包括了各种流体流动的关键部分，并且具有一定的复杂性，以保证仿真结果的可靠性和准确性。
2.2 边界条件为了模拟真实的流体流动情况，我们设置了合适的边界条件。
这些边界条件包括了流体的入口速度、出口压力、壁面摩擦等。
通过调整这些边界条件，我们可以模拟不同工况下的流体流动情况。
2.3 流体属性在仿真过程中，我们需要指定流体的物理属性，如密度、粘度等。
这些属性值对于模拟流体流动的精确性和可靠性非常重要。
3. 计算方法3.1 Navier-Stokes方程在本次仿真分析中，我们采用了Navier-Stokes方程来描述流体流动的行为。
Navier-Stokes方程是基于质量守恒和动量守恒原理推导出来的，可以有效地描述流体流动的各种现象。
3.2 数值方法为了求解Navier-Stokes方程，我们采用了某种数值方法。
在本次仿真分析中，我们选择了合适的数值方法，并结合计算资源的可用性进行计算。
3.3 网格划分为了进行计算，我们需要将几何模型划分为离散网格。
这些网格用于存储流体的各种物理量，并且可以在计算过程中进行相应的更新。
4. 结果分析4.1 流体流动分布通过仿真计算，我们得到了流体流动的分布情况。
根据模拟结果，我们可以观察到不同位置的流速、压力等参数的变化规律，进而分析流体在几何模型中的流动情况。
4.2 流体特性分析基于仿真结果，我们可以对流体的特性进行深入分析。
这些特性包括了流体的速度分布、压力分布、湍流情况等，可以为相关领域的研究提供重要的参考依据。
4.3 敏感性分析在某些情况下，我们可能需要对模型参数进行敏感性分析。
流体力学实验报告
实验一 柏努利实验一、实验目的1、通过实测静止和流动的流体中各项压头及其相互转换，验证流体静力学原理和柏努利方程。
2、通过实测流速的变化和与之相应的压头损失的变化，确定两者之间的关系。
二、基本原理流动的流体具有三种机械能：位能、动能和静压能，这三种能量可以互相转换。
在没有摩擦损失且不输入外功的情况下，流体在稳定流动中流过各截面上的机械能总和是相等的。
在有摩擦而没有外功输入时，任意两截面间机械能的差即为摩擦损失。
流体静压能可用测压管中液柱的高度来表示，取流动系统中的任意两测试点，列柏努利方程式：∑+++=++f h p u g Z P u g Z ρρ2222121122对于水平管，Z 1=Z 2，则 ∑++=+f h p u p u ρρ22212122若u 1=u 2, 则P 2<P 1；在不考虑阻力损失的情况下，即Σh f =0时，若u 1=u 2, 则P 2=P 1。
若u 1>u 2 , p 1<p 2；在静止状态下，即u 1= u 2= 0时，p 1=p 2。
三、实验装置及仪器图2－2 伯努利实验装置图装置由一个液面高度保持不变的水箱，与管径不均匀的玻璃实验管连接，实验管路上取有不同的测压点由玻璃管连接。
水的流量由出口阀门调节，出口阀关闭时流体静止。
四、实验步骤及思考题3、关闭出口阀7，打开阀门3、5，排出系统中空气；然后关闭阀7、3、5，观察并记录各测压管中的液压高度。
思考：所有测压管中的液柱高度是否在同一标高上？应否在同一标高上？为什么？4、将阀7、3半开，观察并记录各个测压管的高度，并思考：（1）A、E两管中液位高度是否相等？若不等，其差值代表什么？（2）B、D两管中，C、D两管中液位高度是否相等？若不等，其差值代表什么？5、将阀全开，观察并记录各测压管的高度，并思考：各测压管内液位高度是否变化？为什么变化？这一现象说明了什么？五、实验数据记录.液柱高度 A B C D E阀门关闭半开全开实验二 雷诺实验一、实验目的1、 观察流体在管内流动的两种不同型态，加强层流和湍流两种流动类型的感性认识；2、掌握雷诺准数Re 的测定与计算；3、测定临界雷诺数。
流体力学报告
流体力学结课报告摘要：流体力学主要研究在各种力的作用下，流体本身的状态，以及流体和固体壁面、流体和流体间、流体与其他运动形态之间的相互作用的力学分支。
流体力学发展简史人类同洪水斗争的历史，可以追溯到遥远的上古时期。
在中国古代的典籍中，就有相传4000多年以前大禹治水，“疏壅导滞”使滔滔洪水各归于河的记载。
先秦时期（公元256～公元251）在四川岷江中游建都江堰，从此成都平原“水旱从人，不知饥馑，时无荒年”。
隋朝自文帝始，历二世（公元584～610），修浚并贯通南北大运河，“自是天下利于转输”，“运漕商旅，往来不绝”。
又如隋大业年间（公元605～公元617），工匠李春在交河上建赵州桥，这座石拱桥的跨径37.4米，拱背上还有4个小拱，既减轻了主拱的负载，又可泄洪，迄今为止1380年依然完好。
历史上，这些伟大的工程，皆因“顺应水性”，，才能跨江河逾千年而不毁。
对流体力学学科的形成作出第一个贡献的是古希腊的阿基米德，他建立了包括物理浮力定律和浮体稳定性在内的液体平衡理论，奠定了流体静力学的基础。
此后千余年间，流体力学没有重大发展。
直到15世纪，意大利达·芬奇的著作才谈到水波、管流、水力机械、鸟的飞翔原理等问题；17世纪，帕斯卡阐明了静止流体中压力的概念。
但流体力学尤其是流体动力学作为一门严密的科学，却是随着经典力学建立了速度、加速度，力、流场等概念，以及质量、动量、能量三个守恒定律的奠定之后才逐步形成的。
流体力学成为一门独立的科学，是在16世纪欧洲文艺复兴以后至18世纪中叶。
这个时期，作为近代自然科学基础的经典力学已相当成熟，为流体力学的建立，奠定了理论基础。
资本主义工商业的崛起，航海、造船、水利以及城市建设等新兴产业的要求，是流体力学建立和发展的推动力。
20世纪以来，随着航空技术的发展，以及大型水利工程、环境工程的需要，流体力学得到了空前的发展。
近年来，由于科学技术的飞速进步，流体力学与其他学科相互渗透，形成了一系列边缘学科，如电磁流体力学、化学流体力学、生物流体力学、高温气体力学、爆炸力学等等。
流体力学实验报告(建能)
流体力学实验指导书与报告姓名：学号：班级：组员：2019年9月目录1 静（水）压强实验 (1)2 能量方程实验..... . (8)3局部阻力系数测定实验 (18)4沿程阻力系数测定实验 (24)5 毕托管（测速）实验 (30)1 静(水)压强实验一、实验目的和要求1.掌握用测压管测量流体静压强的技能；2.验证不可压缩流体静力学基本方程；3.测定油的密度；4.通过对诸多流体静力学现象的实验观察分析，加深流体静力学基本概念理解，提高解决静力学实际问题的能力。
二、实验装置1．实验装置简图实验装置及各部分名称如图1所示。
图.1 流体静力学综合型实验装置图1. 测压管2. 带标尺测压管3. 连通管4. 通气阀5. 加压打气球6. 真空测压管7. 截止阀8. U型测压管9. 油柱10. 水柱11. 减压放水阀说明：下述中的仪器部件编号均指实验装置图中的编号，如测管2即为图1中“2. 带标尺测压管”。
后述各实验中述及的仪器部件编号也均指相应实验装置图中的编号。
2. 装置说明(1) 流体测点静压强的测量方法之一——测压管流体的流动要素有压强、水位、流速、流量等。
压强的测量方法有机械式测量方法与电测法，测量的仪器有静态与动态之分。
测量流体点压强的测压管属机械式静态测量仪器。
测压管是一端连通于流体被测点，另一端开口于大气的透明管，适用于测量流体测点的静态低压范围的相对压强，测量精度为1mm 。
测压管分直管型和“U ”型。
直管型如图1中管2所示，其测点压强p gh ρ=，h 为测压管液面至测点的竖直高度。
“U ”型如图中管1与管8所示。
直管型测压管要求液体测点的绝对压强大于当地大气压，否则因气体流入测点而无法测压；“U ”型测压管可测量液体测点的负压，例如管1中当测压管液面低于测点时的情况；“U ”型测压管还可测量气体的点压强，如管8所示，一般“U ”型管中为单一液体（本装置因其它实验需要在管8中装有油和水两种液体），测点气压为p g h ρ=∆，∆h 为“U ”型测压管两液面的高度差，当管中接触大气的自由液面高于另一液面时∆h 为 “+”，反之∆h 为“-”。
流体静力学实验报告-流体静力学实验步骤
(3)在U形测压管5中,自水面向下深度为12.8cm的一段水柱也是真空区.这段高度与测压管液面低于水箱液面的高度相等,亦与真空测压管液面高于水杯液面高度相等,均为12.8 cm.
3．若再备一根直尺，试采用另外最简便的方法测定 。
2．当 时，试根据记录数据确定水箱内的真空区域。
以p0<0时,第二次B点测量数据(表1-1)为例,此时
pB/γ=-0.6cm<0,相应容器的真空区域包括以下三部分:
(1)过测压管液面(Δh=1.2cm)做一水平面,由等压原理知,相对测压管及水箱内的水体而言,该水平面为等压面,均为大气压强,故该平面以上由密封的水/气所占的空间区域,均为真空区域.
1.掌握用液式测压及测量流体静压强的技能。
2.验证不可压缩流体静力学基本方程，加深对位置水头，压力水头和测压管水头的理解。
3.观察真空度（负压）的生产过程，进一步加深对真空度的理解。
4.测量油的相对密度。
5.通过对诸多流体静力学现象的实验分析，进一步提高解决静力学实际问题的能力。
二、实验原理
1.在重力作用下不可压缩流体静力学基本方程。
最简单的方法,是用直尺分别测量水箱内通大气的情况下,U形测压管油水界面至水面和油水界面至油面的垂直高度hw和ho,由式γwhw=γoho,从而求得do.
4．如测压管太细，对测压管液面的读数将有何影响?
一般说来,当玻璃测压管内径大于10mm时,毛细管影响可以忽略不计,另外,当水质不好时,σ减小,毛细高度较纯净水也减小,当采用有机玻璃管作为测压管时,浸润角θ较大,其h较普通玻璃小.如果用同一根测压管测量液体相对压差值,则毛细现象无任何影响.因为测量高低压强时均有毛细现象,在计算差时,相互抵消了.
板式塔流体力学性能测定-实验报告
化学实验教学中心
实验报告
化学测量与计算实验Ⅱ实验名称：板式塔流体力学性能测定实验报告
学生姓名：学号：
院（系）：年级：级班
指导教师：研究生助教：
实验日期： 2017.05.25 交报告日期： 2017.06.01
图1 筛孔塔板干板压头降Δh d 与筛孔速度u a 之间的关系图2 板式塔的Δh 与空塔速度的关系曲线
塔板上形成稳定液层后，塔板上气液两相的接触和混和状态，也将随着气速的改变而发生变化。
当气速较较小时，气体以鼓泡方式通过液层。
随着气速增大，鼓泡层逐渐转化为泡沫层，并在液面上形成的雾沫层也将随之增大。
对传质效率有着重要作用的因素是充气液层的高度及其结构。
充气液层的结构通常用其平均密度大小来表示。
如果充气液层的气体质量相对于液体质量可略而不计，则
h fρf= h1ρl(4)
调节阀和孔板流量计进入塔底。
通过塔板的尾气由塔顶排出。
气体通过塔板的压力降由压差计显示。
图3 筛板塔
1.塔体；
2.筛孔塔板；
3.漏液排放口；
4.温度计；
5.溢流装置
图4 板式塔流动特性实验装置流程
空气源；2.放空阀；3.消声器；4.孔板流量计；5.U型水柱压差计；6. U型汞柱压差计；
7.板式塔；转子流量计；9. U型水柱塔压差计；10.高位槽；11.排水管。
流体力学实验报告材料
流体力学实验组班级化33吴凡灿学号2013011925成绩 实验时间第6周周日同组成员芦琛琳、董晓锐一、实验目的1、观察塔板上气液两相流动状况，测量气体通过塔板的压力降与空塔气速的关系；测定雾沫夹带量、漏液量与气速的关系；2、研究板式塔负荷性能图的影响因素，作出筛孔塔板或斜孔塔板的负荷性能图；比较筛孔塔板与斜孔塔板的性能；3、观察填料塔气液两相流动状况，测定干填料及不同液体喷淋密度下填料层的阻力降与空塔气速的关系；4、测定填料的液泛气速，并与文献介绍的液泛关联式比较；5、测定一定压力下恒压过滤参数K 、q e 和t e ；6、测定压缩性指数S 和物料特性常数K 。
二、实验原理1．板式塔流体力学特性测定 塔靠自下而上的气体和自上而下的液体逆流流动时相互接触达到传质目的，因此，塔板传质性能的好坏很大程度上取决于塔板上的流体力学状态。
当液体流量一定，气体空塔速度从小到大变动时，可以观察到几种正常的操作状态：鼓泡态、泡沫态和喷射态。
当塔板在很低的气速下操作时，会出现漏液现象；在很高的气速下操作，又会产生过量液沫夹带；在气速和液相负荷均过大时还会产生液泛等几种不正常的操作状态。
塔板的气液正常操作区通常以塔板的负荷性能图表示。
负荷性能图以气体体积流量（m 3/s ）为纵坐标，液体体积流量（m 3/s ）为横坐标标绘而成，它由漏液线、液沫夹带线、液相负荷下限线、液相负荷上限线和液泛线五条线组成。
当塔板的类型、结构尺寸以及待分离的物系确定后，负荷性能图可通过实验确定。
传质效率高、处理量大、压力降低、操作弹性大以及结构简单、加工维修方便是评价塔板性能的主要指标。
为了适应不同的要求，开发了多种新型塔板。
本实验装置安装的塔板可以更换，有筛板、浮阀、斜孔塔板可供实验时选用，也可将自行构思设计的塔板安装在塔上进行研究。
筛板的流体力学模型如下： 1) 压降l c p p p ∆+∆=∆式中，Δp —塔板总压降，Δp c —干板压降，Δp l —板上液层高度压降， 其中20)(051.0c u g p v c ρ=∆式中ρv —气相密度，kg/m 3；g —重力加速度，m/s 2，u 0—筛孔气速，m/s ，c 0—筛孔流量系数，筛板上因液层高度产生的压降Δp l 即液层有效阻力h l ：l l l gh p ρ=∆式中ρl —液相密度，kg/m 3，g —重力加速度,m/s 2，h l —液层有效阻力，m 液柱。
流体力学动量方程实验报告
流体力学动量方程实验报告流体力学动量方程实验报告引言：流体力学是研究流体运动规律的学科，其中动量方程是描述流体运动的基本方程之一。
本实验旨在通过实验验证流体力学动量方程，并探究不同因素对流体运动的影响。
实验设备与方法：1. 实验设备：本实验使用的设备包括流体力学实验装置、流速计、压力计等。
2. 实验方法：首先，将流体力学实验装置设置在水平台面上，并校准流速计和压力计。
然后，通过调节装置中的阀门控制流体的流速和压力。
在实验过程中，记录不同条件下的流速和压力数据，并进行数据处理。
实验结果与分析：1. 流体速度与动量的关系：在实验中，我们通过改变流体的流速，记录了不同流速下的动量数据。
结果显示，流体的动量与流速成正比关系。
这符合流体力学动量方程中的基本原理，即动量等于质量乘以速度。
2. 流体压力与动量的关系：在实验中，我们通过改变流体的压力，记录了不同压力下的动量数据。
结果显示，流体的动量与压力成正比关系。
这也符合流体力学动量方程中的基本原理，即动量等于质量乘以速度。
3. 流体密度与动量的关系：在实验中，我们通过改变流体的密度，记录了不同密度下的动量数据。
结果显示，流体的动量与密度成正比关系。
这同样符合流体力学动量方程中的基本原理。
4. 流体粘度对动量的影响：在实验中，我们通过改变流体的粘度，记录了不同粘度下的动量数据。
结果显示，流体的动量与粘度成反比关系。
这是因为高粘度的流体阻力大，导致动量的损失较大。
结论：通过本实验，我们验证了流体力学动量方程，并研究了不同因素对流体运动的影响。
实验结果表明，流体的动量与流速、压力、密度和粘度等因素密切相关。
这对于理解和预测流体运动具有重要意义，也为相关工程应用提供了理论依据。
未来展望：在今后的研究中，可以进一步探究其他因素对流体运动的影响，如温度、流道形状等。
同时，可以结合数值模拟方法，对流体力学动量方程进行更深入的研究，以提高流体力学理论的准确性和应用价值。
结语：通过本实验，我们对流体力学动量方程有了更深入的理解。
流体力学实验报告(全)参考模板
工程流体力学实验报告实验一流体静力学实验实验原理在重力作用下不可压缩流体静力学基本方程或(1.1)式中：z被测点在基准面的相对位置高度；p被测点的静水压强，用相对压强表示，以下同；p0水箱中液面的表面压强；γ液体容重；h被测点的液体深度。
另对装有水油（图1.2及图1.3）U型测管，应用等压面可得油的比重S0有下列关系：(1.2)据此可用仪器（不用另外尺）直接测得S0。
实验分析与讨论1.同一静止液体内的测管水头线是根什么线？测压管水头指，即静水力学实验仪显示的测管液面至基准面的垂直高度。
测压管水头线指测压管液面的连线。
实验直接观察可知，同一静止液面的测压管水头线是一根水平线。
2.当P<0时，试根据记录数据，确定水箱内的真空区域。
B，相应容器的真空区域包括以下三部分：(1)过测压管2液面作一水平面，由等压面原理知，相对测压管2及水箱内的水体而言，该水平面为等压面，均为大气压强，故该平面以上由密封的水、气所占的空间区域，均为真空区域。
(2)同理，过箱顶小水杯的液面作一水平面，测压管4中，该平面以上的水体亦为真空区域。
(3)在测压管5中，自水面向下深度某一段水柱亦为真空区。
这段高度与测压管2液面低于水箱液面的高度相等，亦与测压管4液面高于小水杯液面高度相等。
3.若再备一根直尺，试采用另外最简便的方法测定γ。
最简单的方法，是用直尺分别测量水箱内通大气情况下，管5油水界面至水面和油水界面至油面的垂直高度h和h0，由式，从而求得γ0。
4.如测压管太细，对测压管液面的读数将有何影响？设被测液体为水，测压管太细，测压管液面因毛细现象而升高，造成测量误差，毛细高度由下式计算式中，为表面张力系数；为液体的容量；d为测压管的内径；h为毛细升高。
常温(t=20℃)的水，=7.28dyn/mm，=0.98dyn/mm。
水与玻璃的浸润角很小，可认为cosθ=1.0。
于是有(h、d单位为mm)一般来说，当玻璃测压管的内径大于10mm时，毛细影响可略而不计。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
.word资料. 计算流体力学
—— 计算机编程报告
班级: 09083201 学号: 0908320125 : 马海涛
一.题目：非定常热传导方程 .word资料.
∂u∂t=γ∂2u
∂x2 (γ>0)
初值问题的有限分差求解。求解域(x,t)∈[0,1]×[0,∞],初值条件和边界条件为
{u(x,0)=f(x)u(0,t)=a(t)u(1,t)=b(t) 其中，γ=1。处置边界条件具体取法为
F(x)= {0 0<𝑥<0.3 1 0.3≤x≤0.7−103x+103 0.7≤x≤1.0 取网格点数Mx=101,要求计算t=0.01,0.1,1,10时的数据。计算中时间步长的取法为𝑥=γ∆t∆x2=0.1,0.5,1。 差分格式：FTCS格式，BTCS格式 二.分析：编制通用过程，给定参数由手动输入获得，将给定参数带入相应方
程得到不同格式不同参数下的数值解，再将数据输出进行处理。程序界面如下：
三.编程过程:
（1）FTCS格式:根据FTCS格式解法第一步利用边界条件置n=0，为差分格式赋初值，编程步骤为： .word资料.
For i = 1 To 100 '赋初值 x(i) = x(i - 1) + 0.01 u(1, i) = 0 u(0, i) = 0 Next For i = 0 To 100 '求函数 If x(i) < 0.3 Then f(i) = 0 ElseIf x(i) >= 0.3 And x(i) <= 0.7 Then f(i) = 1 ElseIf x(i) > 0.7 And x(i) < 1 Then f(i) = (-10 / 3) * x(i) + 10 / 3 End If Next For i = 0 To 100 '为数组u赋初值 u(i, 0) = f(i) Next
第二步：利用循环求解点数值解，并将相应边界条件带入，在计算t
值大于等于输入t值时停止循环，因为迭代循环的一个数据只与它的前一个数据有关，为了节省存在编程时循环当中不断地用本次得到的数据赋给上一个变量，以达到释放前n-2个变量所占用的存。具体方式为： Do '迭代循环 u(0, 1) = 0 u(100, 1) = 0 for k = 1 To 99 '迭代方程 u(k, 1) = c * u(k + 1, 0) + (1 - 2 * c) * u(k, 0) + c * u(k - 1, 0) Next For k = 0 To 100 u(k, 0) = u(k, 1) ‘变量替换 Next k t = t + c * 0.01 * 0.01 Loop While t < z
第三步：将所得数据输出，利用Excel作出相应图像。其图像如下：
.word资料.
00.20.40.60.81012FTCS0.01--0.1FCTS0.01--0.100.20.40.60.81012
FTCS0.01--0.5
FCTS0.01--0.5
-6E+46-4E+46-2E+4602E+464E+466E+46012FTCS0.01--1FCTS0.01--100.050.10.150.20.250.30.350.4012
FTCS0.1--0.1
FCTS 0.1--0.1
00.10.20.30.4012FTCS0.1--0.5FCTS 0.1--0.500.000020.000040.00006012
FTCS1--0.1
FCTS 1--0.1
00.000010.000020.000030.000040.000050.00006012FTCS1--0.5FCTS 1--0.502E-444E-446E-448E-441E-431.2E-431.4E-43012
FTCS10--0.1
FCTS 10--0.1 .word资料.
（2）BTCS格式： 第一步：赋初始值，置n=0，用循环带入初始条件得到uk0的值 程序如下：
For i = 1 To 100 '赋初值 x(i) = x(i - 1) + 0.01 Next For i = 0 To 100 '求函数 If x(i) < 0.3 Then f(i) = 0 ElseIf x(i) >= 0.3 And x(i) <= 0.7 Then f(i) = 1 ElseIf x(i) > 0.7 And x(i) < 1 Then f(i) = (-10 / 3) * x(i) + 10 / 3 End If u(i, 0) = f(i) u(0, i) = 0 u(1, i) = 0 Next For i = 2 To 100 '为矩阵元素赋初值 b(i) = -1 - 2 * w c(i) = w a(i) = w Next For i = 1 To 99 '为矩阵元素中间变量赋初值 d(i, 0) = -u(i, 0) Next
第二步：同样采用循环迭代，并利用FTCS中变量替换的方式节省存,
02E-444E-446E-448E-441E-431.2E-431.4E-4300.20.40.60.811.2
FTCS 10--0.5
FCTS 10--0.5 .word资料. 同时以相应步长对t从零开始累加，直到大于等于输入t值时为止，
程序如下： Do '迭代循环 u(0, 1) = 0 u(100, 1) = 0 For i = 1 To 99 '迭代方程 v(i) = c(i) / (b(i) - v(i - 1) * a(i)) u(i, 0) = (d(i, 0) - u(i - 1, 0) * a(i)) / (b(i) - v(i - 1) * a(i)) Next For i = 99 To 1 Step -1 u(i, 1) = u(i, 0) - v(i) * u(i + 1, 1) Next For k = 0 To 100 '为矩阵元素赋初值 d(k, 0) = -u(k, 1) u(k, 0) = u(k, 1) Next k z = z + w * 0.01 * 0.01 Loop While z < t 第三步：将数据输出并作图表进行数据分析。图表如下：
-0.200.20.40.60.81012BTCS0.01--0.1BCTS 0.01--0.1-0.200.20.4
0.60.81
00.511.5BTCS0.01--0.5
BCTS0.01--0.5 .word资料.
-0.200.20.40.60.81012BTCS 0.01--1BCTS 0.01--1-0.0500.050.10.150.20.250.30.350.4012
BTCS0.1--0.1
BCTS 0.1--0.1
-0.100.10.20.30.400.511.5BTCS0.1--0.5BCTS 0.1--0.5-0.100.10.20.30.4012
BTCS0.1--1
BCTS 0.1--1
-0.0000100.000010.000020.000030.000040.00005012BTCS1--0.5BCTS 1--0.5-0.0000100.000010.000020.000030.000040.00005012
BTCS1--0.1
BCTS 1--0.1 .word资料.
四.总结： 计算参数的选择对结果的影响：向前差分法相对于向后差分法对数的要求较高，当σ=1时向前差分法只有在t=0.01时收敛，其他的t值均不收敛，而对于向后差分法则没有出现这种现象在给定的12组参数下均收敛。对于σ=0.1，,05时两种方法均可以。
-0.0000100.000010.000020.000030.000040.00005012BTCS1--1BCTS 1--1-2E-4202E-424E-426E-428E-421E-41012
BTCS10--0.1
BCTS 10--0.1
-5E-4305E-431E-421.5E-422E-4202
BTCS10--0.5
BCTS 10--0.5-2E-4302E-434E-436E-43012BTCS10--1BCTS 10--1