非定常流体的瞬态动力学CFD分析

工程上,有大量的非定常流体问题无法简单地通过稳态方法来求解。非定常流动主要由两种因素产生,一种是由于流体内部不稳定因素或初始流体状态的非平衡状态,如各种尺度的湍流旋涡,激波,对流等。另一种是由于变换的边界条件或者源项,如脉动流,旋转机械的定转子转动。对于这些非定常流动,就必须通过瞬态的分析方法来了解流体及其固体接触表面的状态。

瞬态分析由于计算多个步长下的流场,计算量会随着模拟的时间长短线性增长。数值上,常用的方法是在每一个时间步长下,计算一个短时间内的稳态,并将这个稳态的计算结果,作为初始条件带入到下一个时间步长的计算中。根据不同的时间求解器,时间步长的选择略有不同,通常隐式求解器的时间步长可以设置的比显式求解器的时间步大。著名的开源多物理场求解器SU2已经被证明可以很好地求解瞬态CFD问题,WELSIM也在2023R3中增加了SU2的支持,详情参见《使用WELSIM生成SU2求解器文件》一文。

瞬态CFD分析的步骤

下面通过实例,演示如何进行瞬态CFD分析。

1. 以二维模型为例,打开WELSIM后,新建一个项目,并将模型设置为2维的流体瞬态模型。

2. 导入几何模型。

3. 划分网格,将最大单元尺寸设置为0.03m。

4. 设置求解的时间步为0.0005秒,总共运行0.6秒。

5. 求解器使用SU2。

6. 控制方程使用可压缩流体的RANS,湍流模型选用Spalart-Allmaras。

7. 设置求解器的相关参数。

8. 设置自由流场区域条件。其中,马赫数为0.3,攻角17度。常温常压。雷诺数为1000。

9. 设置远场边界条件。

10. 设置热流边界条件,数值为零,无热对流。

点击计算按钮,由于是瞬态计算,根据网格密度和时长,需要较长的物理计算时间。计算完成后添加马赫数结果节点,并显示结果云图。下图分别显示流场在0.027秒和0.597秒时的马赫数。

本算例的计算结果视频如下。

SU2+WELSIM可压缩流的湍流计算分析

同时,本算例加入了WELSIM的自动化回归测试库,能够有益于求解器和前端软件的长期维护。测试文件已经开源并分享在GitHub上。地址如下:

https://github.com/WelSimLLC/WelSimAutoTests

WELSIM对SU2命令的支持

SU2的控制文件是主要的求解器输入文件,目前WELSIM作为前处理器已经支持了大量核心命令。SU2的7.5.1版本共含有671个控制命令,WELSIM已经支持了134个命令,占全部命令的20%。已经支持的命令列举如下。

SOLVER

MATH_PROBLEM

KIND_TURB_MODEL

KIND_TRANS_MODEL

BODY_FORCE

BODY_FORCE_VECTOR

RESTART_SOL

FLUID_MODEL

SPECIFIC_HEAT_CP

VISCOSITY_MODEL

MU_CONSTANT

CONDUCTIVITY_MODEL

THERMAL_CONDUCTIVITY_CONSTANT

REYNOLDS_NUMBER

REYNOLDS_LENGTH

PRANDTL_LAM

PRANDTL_TURB

MACH_NUMBER

INIT_OPTION

FREESTREAM_OPTION

FREESTREAM_PRESSURE

FREESTREAM_DENSITY

FREESTREAM_TEMPERATURE

FREESTREAM_TEMPERATURE_VE

INC_DENSITY_MODEL

INC_ENERGY_EQUATION

INC_DENSITY_INIT

INC_VELOCITY_INIT

INC_TEMPERATURE_INIT

FREESTREAM_VELOCITY

FREESTREAM_VISCOSITY

FREESTREAM_INTERMITTENCY

FREESTREAM_TURBULENCEINTENSITY

FREESTREAM_NU_FACTOR

SIDESLIP_ANGLE

AOA

REF_ORIGIN_MOMENT_X

REF_ORIGIN_MOMENT_Y

REF_ORIGIN_MOMENT_Z

REF_AREA

REF_LENGTH

REF_DIMENSIONALIZATION

MARKER_PLOTTING

MARKER_MONITORING

MARKER_ANALYZE

MARKER_DESIGNING

MARKER_EULER

MARKER_FAR

MARKER_SYM

MARKER_NEARFIELD

INLET_TYPE

INC_INLET_TYPE

MARKER_INLET

MARKER_INLET_SPECIES

MARKER_INLET_TURBULENT

MARKER_SUPERSONIC_INLET

MARKER_SUPERSONIC_OUTLET

MARKER_OUTLET

INC_OUTLET_TYPE

MARKER_ISOTHERMAL

MARKER_HEATFLUX

MARKER_HEATTRANSFER

MARKER_PRESSURE

MARKER_DAMPER

TIME_MARCHING

CFL_NUMBER

CFL_ADAPT

CFL_ADAPT_PARAM

RK_ALPHA_COEFF

TIME_DISCRE_FLOW

TIME_DISCRE_FEM_FLOW

TIME_DISCRE_ADJFLOW

TIME_DISCRE_TURB

LINEAR_SOLVER

LINEAR_SOLVER_PREC

LINEAR_SOLVER_ERROR

LINEAR_SOLVER_ITER

CONV_RESIDUAL_MINVAL

CONV_STARTITER

CONV_CAUCHY_ELEMS

CONV_CAUCHY_EPS

CONV_FIELD

MGLEVEL

MGCYCLE

MG_PRE_SMOOTH

MG_POST_SMOOTH

MG_CORRECTION_SMOOTH

MG_DAMP_RESTRICTION

MG_DAMP_PROLONGATION

NUM_METHOD_GRAD

NUM_METHOD_GRAD_RECON

VENKAT_LIMITER_COEFF

ADJ_SHARP_LIMITER_COEFF

CONV_NUM_METHOD_FLOW

MUSCL_FLOW

SLOPE_LIMITER_FLOW

JST_SENSOR_COEFF

LAX_SENSOR_COEFF

CONV_NUM_METHOD_ADJFLOW

MUSCL_ADJFLOW

SLOPE_LIMITER_ADJFLOW

MESH_FORMAT

MESH_FILENAME

MESH_OUT_FILENAME

CONV_FILENAME

SOLUTION_FILENAME

SOLUTION_ADJ_FILENAME

RESTART_FILENAME

RESTART_ADJ_FILENAME

VOLUME_FILENAME

VOLUME_ADJ_FILENAME

GRAD_OBJFUNC_FILENAME

VALUE_OBJFUNC_FILENAME

SURFACE_FILENAME

SURFACE_ADJ_FILENAME

SURFACE_SENS_FILENAME

VOLUME_SENS_FILENAME

TIME_DOMAIN

TIME_ITER

ITER

RESTART_ITER

TIME_STEP

SCREEN_OUTPUT

HISTORY_OUTPUT

VOLUME_OUTPUT

OUTPUT_WRT_FREQ

OUTPUT_FILES

总结

SU2是一款性能优异且协议友好的开源CFD求解器,可以快速的计算含有湍流的瞬态流体问题。同时,得益与友好的图形化界面,WELSIM令SU2的使用变得简单。WELSIM既可以无缝调用SU2进行求解并显示结果,也可以只生成用户所需的SU2计算输入文件。当前,WELSIM是世界上对SU2支持最好的前后处理软件。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 213,864评论 6 494
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,175评论 3 387
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,401评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,170评论 1 286
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,276评论 6 385
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,364评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,401评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,179评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,604评论 1 306
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,902评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,070评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,751评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,380评论 3 319
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,077评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,312评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,924评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,957评论 2 351

推荐阅读更多精彩内容