| 查看: 2616 | 回复: 8 | |||
格物要致知木虫 (著名写手)
|
[交流]
求教fluent的激波捕捉格式
|
|
虫友好，投的一篇论文里面专家提出来激波捕捉格式问题，我原文里面只写了对流通量采用roe格式，方程离散采用二阶迎风格式。 fluent帮助文档里面说的激波捕捉格式是tvd类型的梯度限制器，参考文献是aiaa-89-0366，这个不是harten和yee的那篇报告，说实话不做cfd也看不懂这篇论文。 想请问一下各位激波捕捉格式能不能写tvd，另外迎风格式和roe格式以及激波捕捉之间有没有关系？ 发自小木虫Android客户端 |
回复此楼» 本帖已获得的红花（最新10朵）
dzprre044» 猜你喜欢
求调剂
已经有3人回复
-
一志愿北京化工大学，初试成绩350求调剂
已经有14人回复
-
复试调剂
已经有4人回复
-
323分（计算机视觉和大模型项目）能直接上手
已经有3人回复
-
311分 22408 求调剂
已经有3人回复
-
320分人工智能调剂
已经有7人回复
-
一志愿郑大0705求调剂
已经有4人回复
-
0703化学
已经有10人回复
-
301求调剂
已经有10人回复
-
306分材料与化工求调剂
已经有6人回复
 |
2楼2016-08-04 18:16:55
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
roe应该属于通量的分裂方法，是怎么用特征波理论来分解通量的。激波的捕捉格式，像有限体积法应指的是控制体面上的QL,QR是怎么通过控制体平均量的Q得到的，muscl是属于tvd(总残差减小)的，另外还有weno,eno是属于tvb（总残差有界）。有界和减小，通俗讲就是减小了震荡，震荡的幅值减小或有界。震荡的问题就是由Q到QR,QL可能产生的，如果Q是光滑的就不会有震荡，但是激波是间断的，数学上是不连续的。最简单的用q到qr可以是插值，但是插值的前提是函数连续，在不连续上进行插值就是震荡。像muscl和weno都叫重构（reconstruct）方法，不能理解为插值方法 发自小木虫Android客户端 |
» 本帖已获得的红花（最新10朵）
3楼2016-08-07 20:10:17
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
对流通量的离散格式一般主要有FVS和FDS两类，FVS如Steger-Warming, van Leer等；FDS格式包括Roe, HLL,HLLC 格式等。此外，另外基于van Leer格式并将压力项进行单独处理就形成所谓的AUSM类格式。主要为了提高空间精度，还有所谓的WENO等高分辨率格式。目前主流的流体商业软件大都基于有限体积法，空间计算精度最多两阶，fluent中涉及的roe格式是一阶,AUSM是二阶，AUSM类格式有多种，fluent给出的参考文献好像是95年Liu的，不太肯定。 为提高激波的捕捉精度，抑制激波前后的耗散和色散，通常会对格式中的修正项引入限制函数(限制器)，以满足TVD性质。引入的不同限制器决定了激波捕捉的精度，FLUENT中包括二阶迎风/二阶中心/三阶MUSCL等，三阶MUSCL不一定满足TVD性质。 |
» 本帖已获得的红花（最新10朵）
4楼2016-08-08 06:08:09
格物要致知
木虫 (著名写手)
- 应助: 96 (初中生)
- 金币: 5518.7
- 散金: 1206
- 红花: 28
- 帖子: 1817
- 在线: 357.8小时
- 虫号: 1569834
- 注册: 2012-01-09
- 专业: 内流流体力学
5楼2016-08-09 09:30:33
格物要致知
木虫 (著名写手)
- 应助: 96 (初中生)
- 金币: 5518.7
- 散金: 1206
- 红花: 28
- 帖子: 1817
- 在线: 357.8小时
- 虫号: 1569834
- 注册: 2012-01-09
- 专业: 内流流体力学
6楼2016-08-09 09:31:35
|
|
7楼2016-08-09 12:12:06
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
TVD只是提供了一种激波的捕捉方法，是一类方法的实现思路。就是不捕捉激波，时间推进方法的cfd也应该使用的是迎风格式。激波捕捉格式应该是某tvd格式或者tvb格式，格式的精度和分辨率高低，直接体现在激波的宽度上。激波是间断的，希望计算出来的也是这样，但是数值上很难做到，直接体体现在激波捕捉格式的精度和分辨率上，格式精度和分辨率如果高，激波可能就分辨在一个网格上，如果低的话，那可能就要几个网格上，激波的捕捉得就不好了。 发自小木虫Android客户端 |
8楼2016-08-09 20:50:22
9楼2016-12-16 23:30:44
