《cfx培训教材06瞬态模拟》由会员分享,可在线阅读,更多相关《cfx培训教材06瞬态模拟(20页珍藏版)》请在金锄头文库上搜索。
1、ansys cfx 培训教材第六节:瞬态模拟安世亚太科技(北京)有限公司原因n自然界几乎所有流动都是瞬态流动!n在以下假设下,可以认为流动是稳态的:n忽略非稳态波动n采用总/时均方法,忽略流动的波动性 (这也是为什么采用湍流模型模化湍流的原因)n在cfd,首选稳态的计算方法n更小的计算代价n更易处理和分析n很多情况下,要求瞬态求解:n空气动力学(aerodynamics), 如飞行器、汽车等 涡脱落(vortex shedding)n旋转机械rotating machinery 转子/静子相互作用, 停车, 飞逸n多相流multiphase flows 自由液面, 空泡动力学n变形计算域def
2、orming domains 缸内燃烧n非稳态热传递unsteady heat transfer 瞬态加热和冷却n更多瞬态流动的起源n自然的非稳定性n流动的不稳定性或非稳定的初始流动状态导致非稳定流动n例如: 自然对流, 所有尺度的湍流涡, 流动波(重力波, 冲击波)n强迫非稳定性n时间相关的边界条件, 源项驱动的非稳定流场n例如: 喷嘴处的脉冲流, 转子-静子节间的相互影响kelvin-helmholtz cloud instabilityrotor-stator interaction in an axial compressor瞬态cfd 模拟n瞬态流域的模拟都是基于一个指定的时间周期n
3、可能的求解方法:n稳态求解 流动变量不随时间改变n时间周期求解 流动变量以某种反复的模式波动n也可以通过指定时间间隔的方法简单的进行流动分析.n自由面流动n冲击波的运动n等.n详细分析关心的量n固有频率 ,如.斯特劳哈尔数(strouhal number)n时均(time-averaged)和/或 均方根(rms)值n与时间相关的参数 (如. 冷却一个热的固体的时间要求, 污染物的扩散时间)n谱数据 快速傅里叶变换(fft)20timestep = 2 sinitial time = 0 stotal time = 20 scoefficient loops = 52468 10 12 14
4、16 18time (seconds)5 coefficient loopsn通过计算不同离散时间点的方法完成瞬态模拟n在每个时间点,都需要进行计算迭代如何求解瞬态问题n对稳态计算也有相似的设置n通常的工作流程:1.设置分析类型为transient2.指定瞬态求解时间和时间步长3.设置物理模型和边界条件边界条件可以随时间变化4.指定初始条件最好采用符合实际的物理条件, 比如稳态计算结果5.指定求解器设置6.设定瞬态结果文件(results files), 瞬态统计(transient statistics), 监测点(monitors points)7.运行求解器如何求解瞬态问题1. 分析类型
5、在目录树上编辑在目录树上编辑analysis type,设置,设置option为为transient设置设置time duration控制求解结束时间控制求解结束时间选项选项:total time求解的总时间求解的总时间time per run不计以前计算的所有时间不计以前计算的所有时间maximum number of timesteps 包括以前所完成的计算时间步包括以前所完成的计算时间步number of timesteps per run对本次求解有效对本次求解有效. 不计以前完成的时间步不计以前完成的时间步2. 时间长度和步长2.时间长度和步长设置设置time step size控制
6、两计算时间点之间的间隔控制两计算时间点之间的间隔选项选项:timesteps / timesteps for the run可采用不同的方式可采用不同的方式, 例如例如0.0010.001, 0.002, 0.002, 0.0035*0.001, 10*0.05, 20*0.06adaptive通过设定收敛标准或通过设定收敛标准或courant数的方法,在规定的数的方法,在规定的范围内,时间步可以动态的改变范围内,时间步可以动态的改变2.时间长度和步长n在瞬态模拟中,时间步尺寸是一个很重要的参数n求解与时间相关的特征的时候,需要时间步较小 实际值实际值时间时间关心的变量关心的变量d dt时间时
7、间关心的变量关心的变量d dt时间步较大,时间步较大, 求解点通常不能反映真实的现象求解点通常不能反映真实的现象小的时间步,小的时间步, 求解点可以反映真实的现象求解点可以反映真实的现象2.时间长度和步长n 为了维持求解器的稳定性,也需要较小的时间步n关心的量可能变化缓慢 (例如. 固体内的稳定变化), 如果其它量(如, 速度)具有较小的时间尺度,就不能采用较大的时间步ncourant数常用于评估时间步:n在一个时间步内,流体通过的网格的数目ncourant数通常介于2 10之间, 但在一些情况下该值稍高也可以接受n平均的courant数和最大courant数在每个时间步求解器.out文件中会
8、显示出来n小的时间步会加速收敛必要的情况下必要的情况下, 边界条件可以是时间的函数边界条件可以是时间的函数速度, 质量流, 压力条件, 温度, 等. 都能作为函数表达式在cel表达式中使用 “t”或者“time”3. 边界条件最好采用与物理实际符合的初值条件最好采用与物理实际符合的初值条件经常采用稳态计算的收敛解作为初值经常采用稳态计算的收敛解作为初值瞬态模拟的初值越准确,瞬态计算得到的初瞬态模拟的初值越准确,瞬态计算得到的初值就越精确值就越精确前面几个时间步可能不收敛前面几个时间步可能不收敛为了保证求解器的稳定性,瞬态计算的最初为了保证求解器的稳定性,瞬态计算的最初阶段需要较小的时间步阶段需
9、要较小的时间步计算获得稳定计算获得稳定(周期模式周期模式)后,可以忽略开始后,可以忽略开始的几个循环的几个循环2468 10 12 14 16time (seconds)4. 初始化residualsntransient scheme定义瞬态项的数值算法n两个隐式时间步格式可选:nfirst order backward euler (更稳定)nsecond order backward euler (更精确)n对大多数瞬态模拟,推荐采用默认的second order backward euler格式ntimestep initialisation控制前一时间步作为后一时间步计算起点的方式n可
10、以与前一求解的时间步初值(previous timestep)保持一致n或求解器可以通过对前一求解的外推(extrapolate ), 以获得一更好的计算起点n在高courant数的情况下,不推荐使用nautomatic (默认)模式 ,根据courant数的大小在两种模式之间进行切换5. 求解控制min. 和和 max. coeff. loops设置每个时间步设置每个时间步内的迭代次数的上下限内的迭代次数的上下限争取每个时间步可以在争取每个时间步可以在3-5个循环内达到收敛个循环内达到收敛complex physics may need more loops如果在最大循环内没有达到收敛,通常
11、来说减如果在最大循环内没有达到收敛,通常来说减小时间步是较好的选择,而不是增加循环次数小时间步是较好的选择,而不是增加循环次数不管收敛与否,求解将转到下一个时间步不管收敛与否,求解将转到下一个时间步监测求解很重要监测求解很重要5.求解控制瞬态结果瞬态结果(transient results)默认只写出最后的默认只写出最后的res文件文件不包含瞬态求解信息不包含瞬态求解信息需要在需要在output控制下定义控制下定义transient resultstransient results选项选项standard写出包括全部结果的文件写出包括全部结果的文件将占据大量的磁盘空间将占据大量的磁盘空间sma
12、llest写出最小文件,该文件依然可以作为重新计算写出最小文件,该文件依然可以作为重新计算的初值文件的初值文件 (文件依然较大文件依然较大)selected variables写出选取的关心的变量,文件将较小写出选取的关心的变量,文件将较小output frequency控制结果文件写出的频率控制结果文件写出的频率6. 输出控制瞬态统计瞬态统计(transient statistics)对求解变量生成求解统计对求解变量生成求解统计option:算法平均:算法平均(arithmetic average), rms, minimum, maximum, standard deviation和和fu
13、ll (everything)拾取关心的变量拾取关心的变量start和和stop iteration list:定义收集统计:定义收集统计开始和结束的时间开始和结束的时间6.输出控制6.输出控制nmonitor points通常用于稳态模拟nmonitor coefficient loop convergence在时间步内为每个迭代创建监测记录n用于观察关心量在时间步内是否在收敛n默认的仅仅显示时间步末的检测值ntip: 监测一个表达式,在求解器管理器中将创建一个瞬态记录图表. 与从瞬态结果文件创建图表,这个方法更加容易,而且不需要写出瞬态结果文件n与稳态输出不同的是,瞬态输出在每个时间步都包括ncourant数信息显示在每个时间步的开始处n通过检测rms和max残差,确保在时间步结束时达到收敛7.输出谢谢a pera global company pera china