CFD中文网

你好，请问这个问题你解决了吗？我也遇到类似问题。如果A程序是稳态问题，B程序是瞬态问题。而且AB还要循环多次求解。请问这个问题该如何解决呢？

@bestucan openfoam 5.0和6.0都测了··一样的崩

各位好，我想把拓扑优化与openfoam中自适应网格相结合。但在mesh.update()运行时报错了。

• 第一次优化后，根据密度场加密网格。此时程序正常运行
  

• 但第二步优化时，加密网格时出错

感觉像是点在映射时出了问题。我把源程序和算例文件都放在这里了，希望各位大佬不吝赐教。谢谢！solver.zip case.zip （mesh.update()在sensitivity.H里）

谢谢两位大佬@史浩 @东岳

大佬，我也遇到和你一样的疑问，为啥G前面多了个rho，不吝赐教万分感谢！！！

问题：在scalarTransportFoam中，目前是把rhoc当做常数，并除到K下面，得到DT

请问如果rho随空间变化（比如说现有 volScalarfield rho），该如何实现该方程的求解呢？

最近在学习solidDisplacementFoam时遇到了一个问题，在计算热应力是，openfoam是这样写的：

       if (thermalStress)
        {
            const volScalarField& T = Tptr();
            sigma = sigma - I*(rho*threeKalpha*T);
        }

可弹性力学中，应该是用温度变化计算热应力，即：

 sigma = sigma - I*(rho*threeKalpha*(T-T0));

其中T0是初始的参考温度。请问有人了解这个问题吗？请不吝赐教

谢谢东岳老师

最近看了东岳老师的pimpleDyMFoam解析，想请问下如果只通过AMR加密，那ub=0吗？

本人做流体优化，其中需要求解伴随控制方程，以及伴随边界条件。现有边界条件如下： 其中u为速度，v为伴随速度。等号右边不用管。
文献中把该边界条件分解边界处为切向与法向两个式子分别添加。我在做的时候切向分解下列一项
该项变为 其中Vn即为伴随速度的法向分量，Xt为切向，请问这个在openfoam中如何实现呢？
补充：文献中对 在切向分解时得到 这一项可以用面心和体心处的V值进行插值（由于是法向Xn），但切向该如何做呢？？请各位赐教。