CFD中文网

http://www.cfd-china.com/topic/3005 看看这个 我不太熟悉超临界流体

最近在入门流体力学中的线性稳定性分析相关工作，我的理解是通过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

@东岳 我现在接触的还很浅。。以后碰到这方面问题再来请教~~

@bestucan 好的，谢谢您。

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

缩小模型只能保证一个无量纲数相同，一般的流动考虑雷诺数就行，缩放会带来一些非真实的物理效应，但是计算量不够就不要讲究了

@李东岳 感谢大佬，我再去了解一下这个URANS

这更像是用代码格式化软件统一生成的。高级编辑器自带的代码格式化一般只改行前缩进，代码格式化软件是全部重新排版。

一般大型代码库都有非常细致的格式要求，代码格式化软件可以非常细致的定义每一个格式化要求：每行最多字符；算符前后加不加空格、换不换行，前换 后换还是都换……

然后用被处理过的代码非常满足每一个细致要求，所以某些地方就看着有一点点别扭了。 indent是GNU出的格式化软件。里面有一些常用的风格模板：GNU、K&R、linux。也可以跟“像一篇文章一样长的”参数来定义自己的代码排版风格。

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

大家好，估计你们都意识到广告条改了。嗯，因为几个月来，只收到几个3位数的广告，不足以支撑。索性自己全用了。

同行想做广告的可以撤了，没戏。我用心血做网站做这么几年，做到月几十万电机量，我不缺打广告那一顿饭钱，有的一顿饭钱都不想给

同行想做项目的可以试试，单打独斗年7位数没问题，前提是你混到CFD行业的塔尖。我个人每年只能做数量有限的项目，多了我也干不过来太占用时间

同行想做培训的也可以试试，请有点良心别瞎忽悠，好多机构只知道赚钱，就跟互联网公司跟风卖买肉夹馍

同行想做咨询的，不好说，这玩意不费事躺着就来钱。但做咨询的不是特别多。小客户这面不见得你竞争得过我。大客户那面我也拿不到。如果你还没开始做，那估计你也够呛

:papa:

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

@浪迹天大 在mathjax里面要设置跳过某些关键词，比如skipTags: ["<pre>","<span>","<code>"]

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 *** 求大神指点一二！

69cdd51f-33c7-4640-9840-dea844d49ba6-image.png
请问，想把这个设置成以红点为转轴，左右摆动的摆片，该怎么设置？

这个是用 fluent 的
https://grahamgroup.che.wisc.edu/pub/baier_thesis.pdf

这个是用 OpenFOAM 的
https://www.researchgate.net/publication/257986458_CFD_Simulation_of_Liquid-Liquid_Extraction_Columns_and_Visualization_of_Eulerian_Datasets

设初始压强和出口边界压强都为目标值，然后，内部气压低了，就会从出口回流进来空间。出口是恒压的，出口和内部是联通的，就能保证内部恒压。

请问各位有遇到过这个问题吗？谢谢各位答疑

Screenshot from 2022-11-11 10-57-06.png

@tidedrinker 感谢您的回复，确实如您所言，可以通过修改旋转中的速度矢量来改变旋转轴，而上面的那个参数确实是旋转中心位点