采用大涡模型进行风工程计算及方法1.docVIP

前言 风工程计算中流场尺度大，建筑风工程的流场计算域可达百米甚至数百米，模型单元数量多，收敛过程需要的迭代步数多。由此导致的计算量大以及实际存在的收敛困难是必须要克服的问题。 在ADINA中，全耦合流场的控制体积或者有限元计算方法尽管精度非常高，但由于风工程模型网格数量多（10~1000万），流场自由度全耦合方法会耗费大量时间。我们意见一般采用FCBI-C算法，即SimpleC算法，在空间离散和时间积分上选择二阶精度算法，还包括其它一些关键建模环节的合理处理，可实现高效，高精度的风工程计算。 以下我们用一个方块障碍的风洞试验模拟来说明主要问题。 示例模型说明 方块障碍物尺寸为1*1*1米，风洞高度2米。计算模型初步尺寸设定为，上游取3米，下流6米，侧面3米。 网格划分尺度参加附带模型。初步计算的网格单元数量为57408。 分析一：常规瞬态风场计算（入口速度不变） 第一步 采用K-E湍流模型,获得平均速度下的稳态流场流态。只计算一步，输出map文件给后续计算使用，作为后续计算的初始流场状态。主要参数设置如下： ********************************************************** * SOLUTION PARAMETERS-Step 1 ********************************************************** *总体参数设置：分析类型STEADY-STATE，不可压缩，采用K-E湍流模型，输出map文 *件MAP-OUTP=YES *MAP-FORM=YES；采用FCBI-C算法。 MASTER ANALYSIS=STEADY-STATE MODEX=EXECUTE TSTART=0.00000000000000, IDOF=1 TURBULEN=KE HYDRO=YES STREAM=YES TRACTB=DEFAULT, IRINT=DEFAULT AUTOMATI=NO SOLVER=AMG1 COMPRESS=NO FSINTERA=NO, NMASS=0 MASSCOUP=NO MAP-OUTP=YES MAP-FORM=YES NONDIMEN=YES, MAXSOLME=0 MTOTM=2 RECL=3000 ALE=NO THERMAL-=NO, UPWINDIN=CONTROL-VOLUME MESHUPDA=CURRENT MESHADAP=NO, COUPLING=ITERATIVE POROUS-C=NO CELL-BCD=YES VOF=NO FCBI=CENTER, TURB-ITE=COUPLED EM-MODEL=NO *无量纲化：注意长度单位为流场核心区域单元平均边长，速度为入口平均速度。 NONDIMENSION INPUT=YES OUTPUT=NO L=.05, XREF=.0 YREF=.0 ZREF=.0, V=20.0 D=1.0 CP=1.0, T=1.0 TREF=.0 *方程残差为0.0001，方程类型为质量方程；参数残差0.0001，参数类型为所有参数；速度 *压力场耦合算法SIMPLEC；空间离散精度一阶First；最大迭代次数999。 OUTER-ITERAT EPSER=1.E-4 ERTYPE=MASS, EPSVR=1.E-4 VRTYPE=ALL MAXITE=999 MINITE=5, PVCOUP=SIMPLEC ORDER=FIRST FTSTYPE=NO FTS=1.E+21, VPT-MAXI=1 VP-MAXIT=1 KE-MAXIT=1 VIDPRT=DEFAULT ITSMAX=0, MS-SOLVE=DEFAULT *松弛因子：速度0.35,压力0.85，湍流动能 0.4,耗散率0.4 。 OVAR-CONTROL RELAX-PRI-SS VELOCITY=0.35, PRESSURE=0.85 TEMPERAT=0.95, TURB-K=0.4 TURB-E=0.4 * *************************************************************** * INNER ITERATION SCALES-内迭代过程设置 *************************************************************** *最少迭代次数：速度4，压力25，湍流5。 IVAR-