他可能要根据 1000-1500 确认 1000 没出错，能用:chigua2:

可以写个 bash 脚本，里面依序写上你要进行的任务：

计算 修改controlDict 计算 修改controlDict

可以写个循环，修改 controlDict 可以用 sed，查找 startTime 0， 替换。

@李东岳
好的，我试试。感谢回复:146:

@miss-zhou_dut 哎哟我勒个去 :143:

@李东岳 谢谢李老师😄

@李东岳 请问老师，可以实现计算一定步数后，时间步倒退，重新计算一遍嘛？我想尝试计算1000时间步，然后从第500步计算到1500步，再从第1000步计算到2000步，这样可以实现嘛？

@李东岳 李老师，再请教一个问题，在检查网格的时候，cell zones 和 Number of regions的数目必须一样吗？如果不一样的话，请问代表什么问题？

@callsy 关于这个问题你看东岳老师发的链接，我之前也试过这个链接，有时候能解决，有时候不行，不行的话就得重新画网格，还有就是切换版本，有时候我用2021R1版本会报错，但是用2019R2就没问题

求助各位老师，我试着在sixDoFRigidBodyMotion库中建立一个新的mylinearSpring文件夹（类似于原linearSpring），意图改动其中的锚点anchor的位置。先试着让它z坐标加1，在编译时出错，错误信息如下。

sixDoFRigidBodyMotion/restraints/mylinearSpring/mylinearSpring.C:error: assignment of read-only location ‘((const Foam::sixDoFRigidBodyMotionRestraints::mylinearSpring*)this)->Foam::sixDoFRigidBodyMotionRestraints::mylinearSpring::anchor_.Foam::Vector<Cmpt>::z<double>()’ anchor_.z()=anchor_.z()+1;

提示在给只读位置幅值，是不是意味着在这里anchor是不能变动的量，那要怎样才能改动它呢？

@李东岳 我试了一下，可以的，是我想要的结果。谢谢东岳老师！

@leeck Q562658790

@李东岳 懂了，谢谢李老师

multiphaseEulerfoam，低离散相分数的时候，没问题，离散相增加到0.35，就报错，尝试修改了边界条件、相间作用力模型、离散格式和网格都无果，差别就是跑的时间长短，最终都发散报错如下：

[8] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [8] #1 Foam::sigFpe::sigHandler(int) at ??:? [8] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" [8] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? [8] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/zhl/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" [8] #5 Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel>::nu() const at ??:? [8] #6 Foam::RASModels::kEpsilon<Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> >::correct() at ??:? [8] #7 Foam::phaseSystem::correctTurbulence() at ??:? [8] #8 ? in "/home/zhl/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" [8] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [8] #10 ?[ubuntu:09382] *** Process received signal *** 或者这样： [8] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [8] #1 Foam::sigFpe::sigHandler(int) at ??:? [8] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" [8] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? [8] #4 Foam::operator/(Foam::tmp<Foam::Field<double> > const&, Foam::UList<double> const&) at ??:? [8] #5 Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel>::nu(int) const at ??:? [8] #6 Foam::epsilonWallFunctionFvPatchScalarField::calculate(Foam::momentumTransportModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:? [8] #7 Foam::epsilonWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::momentumTransportModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:? [8] #8 Foam::epsilonWallFunctionFvPatchScalarField::updateCoeffs() at ??:? [8] #9 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::updateCoeffs() at ??:? [8] #10 Foam::RASModels::kEpsilon<Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> >::correct() at ??:? [8] #11 Foam::phaseSystem::correctTurbulence() at ??:? [8] #12 ? in "/home/zhl/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" [8] #13 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [8] #14 ?[ubuntu:07707] *** Process received signal *** 求大神指点一二！

@xpqiu 不是rpm包，是自己编译。

现在不管了，交给运维了，1000块，自己瞎搞，把同事的监测搞没了，被骂了一顿...

@冠竹 在 .h文件中参照nu的形式，把你新加的DT声明一下就好。建议去了解一下C++变量的声明和定义，，

@李东岳 ::146: 谢谢东岳大神

最近在入门流体力学中的线性稳定性分析相关工作，我的理解是通过CFD计算获得的流场作为base flow,对N-S方程进行线性化，获得小扰动特征信息。目前在全网找到了一篇入门相关的文献：https://www.sciencedirect.com/science/article/pii/S1000936116300024
我现在的问题是论文中提到A是右侧线性化的N-S方程的雅克比矩阵（如图1所示），那么如何在openfoam中获得该雅克比矩阵？
同时通过查阅另外一篇文献，文献2（https://arxiv.org/abs/2111.07356）给出了矩阵A的相关表达式，如图2所示。
PS1：通过查阅cfdonline，目前没有发现openfoam中附加线性稳定性分析相关的toolbox，而Nek5000中有相关的后处理模块，但是有限体积法目前刚入门算熟练，实在不想转去谱元法，所以想在openfoam中解决这个问题。
PS2:希望东岳老师能给点suggestion，最好能赐几行codes

图片1.png

@李东岳，另外看东岳老师最近公司在招人，说实话想问下收不收实习生呀，只要能学到技术就很满足！

图片2.png

大家好 想请教一个问题 超音速流计算机翼时，算了大概0.7秒发散了，显示好像是温度为负数了，用的rhoCentralFoam,这个有啥解决方法吗？谢谢大家:chouchou: Screenshot from 2022-03-10 21-44-24.png

@李东岳 好的，谢谢老师，麻烦您了:140:

@桥下卧龙 你好，请问一下，你的网格是用ICEM画的吗？

可以参考这篇文献：Determination of the Angle of Attack on Rotor Blades