CFD中文网

@李东岳 是的，之前搞过飞机结冰的模拟。现在已远离学术圈哈哈

哈哈的确是的，当年实习时候做的项目。希望各位老铁/老师 多多引用

@李东岳 东岳前辈觉得硕士/博士应届毕业，平均有30万么？有没有对行业选择的一些建议？感谢感谢！

@liujunCFD 都是海淀区的985

大家好，

想开一个就业贴，大家来讨论一下国内CFD就业情况，互相学习~

以下是个人的一些认知：国内CFD就业市场像科研一样也分为两类：第一类是使用CFD技术为其他产品开发做支撑的，这类就是CFD的应用。比如系统热分析，化工，暖通，整车等等。往往CFD不是主业，并且主要使用商软，Fluent，Star-ccm等。一般title是 XX工程师。第二类是关于CFD算法的开发，对学历要求高一些，要会写code，有算法植入的经历。

个人感觉对于应届毕业生或者3年经验内的朋友，第一类偏应用的工作，占大多数，但薪资较低。第二类门槛高一些，坑位较少，但薪资高一点。但也有可能，未来5-10年，国家会更重视cfd软件开发，从而导致第二类工作蓬勃发展？？

总体来说，对国内CFD就业环境，机会多少，城市，以及工资待遇方面不是很清楚，毕竟不像码工市场大，大家信息流通快。所以开贴希望大家讨论，也希望论坛里的前辈给新人分享经验，建议。

先说一些我自己及朋友的信息来抛砖引玉：
杭州，汽车行业，ansys工程师，25k-12期, 硕士应届;
成都，ANSYS，应用工程师，16k-12期, 硕士应届;
长春，汽车行业，CFD工程师，20万年薪，硕士应届；
上海，上汽，研发工程师，25-30万年薪，海外硕士应届；

@bestucan 感谢回复！加thermo在incomprehensible里面，理论上可以。不过今天试了一下，报错挺多好多需要改。。。并且主要是我们想做一个普适性的BC, 因为会用到多个不可压缩solvers，每个都改太麻烦了。所以还是想从改bc上入手。

大家好，

最近在用externalCoupledTemperature 边界条件，这个是用来把OpenFOAM solver和其他cfd联系起来。但是目前这个BC只能用在compressible solver里面，原因可见下面的code. 基本上就是，它需要‘he'和’alpha‘，但这些传热量只能在thermo库里找到。因为incomopressible solver是分开求解温度的，只有transport库，所以会报错。

if (db().foundObject<cmpTurbModelType>(turbName))
    {
        const cmpTurbModelType& turbModel =
            db().lookupObject<cmpTurbModelType>(turbName);

        const basicThermo& thermo = turbModel.transport();

        const fvPatchScalarField& hep = thermo.he().boundaryField()[patchi];

        qDot = turbModel.alphaEff(patchi)*hep.snGrad();
    }
    else if (db().foundObject<basicThermo>(thermoName))
    {
        const basicThermo& thermo = db().lookupObject<basicThermo>(thermoName);

        const fvPatchScalarField& hep = thermo.he().boundaryField()[patchi];

        qDot = thermo.alpha().boundaryField()[patchi]*hep.snGrad();
    }
    else
    {
        FatalErrorInFunction
            << "Condition requires either compressible turbulence and/or "
            << "thermo model to be available" << exit(FatalError);
    }

现在想修改一下这个BC, 让它可以在incompressible solver中也得到’he', 'alpha', 还请大家给点建议！十分感谢！

@东岳 几年前在学校用student license跑燃烧，连100万都不能算的。并行的话，还有核数限制。忘记具体多少了。。

@东岳 还没毕业。。都是泪 结冰也在搞，也和别人合作搞搞别的。希望早日上岸。

感谢东岳前辈回复！
刚比较了一下OF5和7的结果，几乎是一样的。就是感觉这个不完全符合一般意义的Boussinesq approximation了。

大家好，之前用的OF5中，heat transfer有4个solvers，OF7把buoyantBoussinesqPimple[Sp]Foam 删掉了，通过把Boussinesq approximation作为 equation of state，使得之前的 buoyantPimpleFoam也可以使用Boussinesq approximation。
但是一般来讲，Boussinesq approximation 只对buoyancy term 进行处理，就像之前版本的buoyantBoussinesqPimple[Sp]Foam里，基于不可压缩，只有动量方程的buoyancy term 有修改过的rho项。现在，在of7里面，buoyantPimpleFoam完全是可压缩的方程，当我们使用Boussinesq approximation作为 equation of state时，就把 rho = rho0*(1 - beta*(T - T0)) 应用到了所有的rho项上。。。
请各位前辈讲解一下这是什么操作？？多谢！

@RebelYoung
建议先查boundary conditions；还有mesh，你的checkMesh结果呢

@RebelYoung
你这个log是rhoC的吧？你的case是inviscid还是viscous flow? rhoCentralFoam对待inviscid flow是用的explicit方法求解，所以假如你没有viscosity，用diagonal solver时，residual就是0.

可能的原因太多了，先检查boundary conditions，你可以把case上传来看一下。
假如试了不同的BCs, 还是有问题。可以先做一些简单尝试。
rhoC 在做高mach number时候，出现这个错误，可以尝试把第一步dt降低很多。
rhoS 可能和simple算法有关，你可以在fvOption里面加个限制温度的：

limitTemperature
{
    type	limitTemperature;
    active	true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        Tmin	300;
        Tmax	400;
    }
}

感觉你的thermo设定可能也有问题，你可以把file放上来看看。
还有mesh，checkMesh结果怎么样？

哈哈专家不敢当。。现在没有现成的求解器，最接近的可能是reactingParcelFilmFoam, 可以做water film在机翼上的模拟，但是没有结冰过程。
Messinger和myers model都不难植入，可以直接做成平衡方程的那种形式，但是宏观的机翼结冰模拟已经做的太多了。。不推荐搞这个。

@东岳 请问东岳老师，可以用速度基求解器比如rhoCentralFoam求解这个算例么？还得加个温度去解压力？
没算过DNS，望指点！谢谢！

@东岳 是啊我也看到了，改成压力基说是在low mach时更好。感觉各种拼凑。。

@东岳 感谢前辈回复！

看来如果要做全隐式，还需要一些其他方法捕获间断。比如改成压力基然后用pimple？

rhoCentralFoam里面隐式的求解应该对速度修正吧？加入粘性项对速度的影响。

这里附上@东岳 前辈的 rhoCentralFoam解析
以及rhoCentralFoam作者的 文章

大家好，

rhoCentralFoam 在处理inviscid flow时完全显式，对于粘性项做隐式处理。

好奇为什么不直接做成 fully implicit? 因为作者在文中也提到考虑到步长问题，将粘性项做了隐式处理。还是说 central-upwind 的 KT/KNP 在全隐式时不可行？

欢迎大家讨论，谢谢！