rhoCentralFoam气体绝热膨胀算例求助
-
大家好。我做了一个rhoCentralFoam求解恒温、恒速、恒压入口气体膨胀的算例。密度和温度运算结果如下图所示：
解域为直径0.4 mm 半圆挖去直径0.1 mm 半圆；入射质量密度 0.76 mg/cc ；入射速度 120 km/s ；背景质量比 1000 :1；入射温度为 1E4 K；背景温度为 3E2 K 。
目前算例膨胀波后缘的温度过高了，与实验不符。实验用的是激光打固体靶，中心温度比边缘高。密度模拟的比较准确。我想将这个算例改为绝热膨胀的，但不知道该怎么修改入口边界条件。恳请大佬赐教。实验室有激光打靶的模拟程序，我个人只是想用OF重现一下。
PS：能否在计算0.2 ns 后将入口改为zerogradient边界，后面时间的模拟是绝热膨胀吗？
