CFD中文网

@李东岳 非常感谢，最近一直卡在这个问题上

最近想在动量方程里面植入一个比较复杂的模型，导过来导过去的，量纲就很复杂，常数项的量纲不好给，想直接关上量纲对齐。

最近要做一些跨临界的问题，但是在跨临界中，如下图，由于伪沸腾的影响，流体的热力学性质会在伪沸腾区呈现剧烈的变化，形成类似于激波的区域。但是在亚临界区域，压力对密度的变化极为敏感，所以传统的全守恒格式会引起压力震荡。很多研究者针对这个问题做了很多修正。比如：求解压力传播方程并加入人工扩散、计算左右通量求黎曼解并对最后的结果熵修正。但是目前的算法总体来说很繁琐，且由于通过压力反推能量会引起能量不守恒。主要有两个问题，非常感谢您的解答：

1. 求解压力传播方程的那个算法很类似simple算法的pressure based方法，为什么现在几乎所有的计算空气动力学论文里面几乎没有使用过simple算法，更倾向于TVD RK这种显示迭代？
2. 显式格式很适合GPU加速计算，想开发一套GPU的代码用来计算。但是GPU对双精度计算的支持很差。当在激波附近使用高阶离散格式的时候如果使用单精度浮点数计算，截断误差是否会对计算结果产生比较大的影响？
   

@bestucan 时间太短了，烧不起来，只能强算，往后算

目前在做微尺度燃烧的一些机理性的研究，需要计算详细的机理，但是只要选择chemkin，计算的过程中甲烷和空气就不会燃烧。同样的条件下使用单步的反应机理就能燃烧，同样的条件下加大反应步长和空间尺度，也能烧的起来。想请教一下大佬们有没有遇到过相似的问题，如能解答，甚是感激。已经在这个问题上卡了两周了。
这个是我的controlDict文件

  application     reactingFoam;
  
  startFrom       latestTime;
  
  startTime       0;
  
  stopAt          endTime;
  
  endTime         0.01;
  
  deltaT          2e-7;
  
  writeControl    timeStep;
  
  writeInterval   1000;
  
  purgeWrite      0;
  
  writeFormat     ascii;
  
  writePrecision  6;
  
  writeCompression off;```