CFD中文网

激光固体靶相互作用在激光脉冲时间内是等温膨胀，激光结束后是绝热膨胀。两个过程我都想用OF试一下

大家好。我做了一个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边界，后面时间的模拟是绝热膨胀吗？

@李东岳 老师我OF-4.1版本算例constant文件夹下得有turbulenceProperties字典，里面写的是laminar，不涉及湍流。好像OF-8以后都不需要turbulenceProperties字典了。

@尚善若水 谢谢老师。我觉得rCF这点不太好，，有没有更新的算法解决这个问题？现在有些既是高阶精度稳定性又好的格式，有能移植到OF的相关文献吗？

@李东岳 李老师好。对于shockTube算例，OF-4.1和OF-6是100mesh的;OF-8是1000mesh。我用OF-4.1算了OF-8的算例，只添加了turbulenceProperties字典，计算结果如图：
原始算例maxCo=0.2时Ux参数仍有震荡，调小maxCo后震荡被抑制。那位老师能给一下OF-8的计算结果吗？可能是版本的问题吗？

两位老师这是设置不同maxCo的计算结果。800网格

调小时间步长后，我觉得振荡被抑制的很好，但没有完全消去。似乎rhoCentralFoam的鲁棒性不是很好？我印象中有些求解器的CFL数设置在0.98都没有数值振荡的。不知道我的理解对不对

感谢李老师。这是我原始算例计算的结果。绿线是100网格./Allrun计算结果，蓝线是将网格改为800然后./Allrun的计算结果。其它参数都没有动。 这是算例shockTube_org.rar 我用的OF-4.1

大家好。我发现用rhoCentralFoam求解shockTube算例时x方向速度会有一些数值震荡。我调了网格数和maxCo数还是消不掉它。计算结果如下面两个图所示：

根据这个网站[Shocktube – rhoCentralFoam TVD Schemes Test – curiosityFluids](Shocktube – rhoCentralFoam TVD Schemes Test – curiosityFluids)我将限制器改为vanAlbada可是仍有数值震荡。
我将shockTube算例也上传了shockTube_test.rar 。我觉的应该哪里的参数配置需要调整一下，可以消除这个数值震荡。希望大佬们指点迷津。