No MRF models present不是问题

你这个看起来是壁面函数修正的时候，有0存在。我看你是高超音速。是不是你算的是无粘流。

一楼那个图，如果不想让速度分散的话，感觉只能在上下两层加一层多孔介质，不知道其他大佬怎么看？

@zzcfd

这个边界条件我也没有用过。不过我刚才看了一下源代码。你看一下mapped和TimeVaryingMappedFixedValue的区别：

//- Name of the field data table, defaults to the name of the field word fieldTableName_; //- If true adjust the mapped field to maintain average value Switch setAverage_; //- Fraction of perturbation (fraction of bounding box) to add scalar perturb_; //- Interpolation scheme to use word mapMethod_; //- 2D interpolation (for 'planarInterpolation' mapMethod) autoPtr<pointToPointPlanarInterpolation> mapperPtr_; //- List of boundaryData time directories instantList sampleTimes_; //- Current starting index in sampleTimes label startSampleTime_; //- Interpolated values from startSampleTime Field<Type> startSampledValues_; //- If setAverage: starting average value Type startAverage_; //- Current end index in sampleTimes label endSampleTime_; //- Interpolated values from endSampleTime Field<Type> endSampledValues_; //- If setAverage: end average value Type endAverage_; //- Time varying offset values to interpolated data autoPtr<Function1<Type>> offset_; Example of the boundary condition specification: \verbatim <patchName> { type mapped; fieldName T; setAverage no; average 0; interpolationScheme cell; value uniform 0; } \endverbatim

TimeVaryingMappedFixedValue看起来就是附加了开始sample的时间以及结束sample时间的mapped边界条件。在时间之外，就是fixedValue。你核实下。

@zz给我趴下 在 OpenFOAM并行运算。 中说：

前几天我服务器里在跑着算例，但办公室停电了。来电之后我想接着跑算例，结果发现算例跑不了了

这是一个很诡异的问题。毫无理由啊

算是解决了吧，其实大可不必费力去猜，先写个Info看看这玩意执行的时间点合不合适，然后直接访问fvMesh对象，拿个引用出来就完事了:mihu:

codeCorrect #{ volVectorField &U = mesh().lookupObjectRef<volVectorField>("U"); #}

这个结果很正常，毕竟很难严格保界，但是这里的-1.1e-11不算越界，毕竟这个数基本都是0

如题，俺知道paraview可以.cgns格式读取并后处理，但是这只是一个时间步。如果想跟openfoam结果一样连续性的读取，该如何操作呢？诸位尝试过吗？

@楼道男孩 这就不太好弄了 :135: 要是openfoam的还能传上来让大家看看

大佬，我刚刚发了一个帖子，能不能帮我看一下，关于煤粉燃烧的。之前你帮我回答过相关问题，感谢了。:xinxin:

大佬，我刚刚发了一个帖子，能不能帮我看一下，关于煤粉燃烧的。之前你帮我回答过相关问题，感谢了。:xinxin:

@疏影横斜水清浅 我只能通过后处理来完成

@颗粒传书 或者你只显示lagrangian cloud

@李东岳 李老师，用的是reactingTwoPhaseEulerFoam里的class method

@小火人 帖子看起来比较老了，我想请教下 对于您这个case，周围空气的进出口边界条件如何设置呢？特别关心的是压力边界条件，totalPressure吗？因为我使用totalPressure设置后结果不准确也不稳定，谢谢

@wangfei9088 好的，谢谢你

@luca 类似下面的写法，参考写下。

top { type cyclicAMI; inGroups List<word> 1(cyclicAMI); nFaces 50000; startFace 21552500; matchTolerance 0.0001; neighbourPatch bottom; transformType translational; separation (0 0 2); //代表由bottom面平移 (0 0 2)得到top面 method faceAreaWeightAMI; } bottom { type cyclicAMI; inGroups List<word> 1(cyclicAMI); nFaces 50000; startFace 21602500; matchTolerance 0.0001; neighbourPatch top; transformType translational; separation (0 0 -2); //代表由top面平移 (0 0 -2)得到bottom面 method faceAreaWeightAMI; }

在官网上下载的paraview可以用了

@李东岳 收到，谢谢东岳老师

@wangfei9088 非常感谢

@李东岳 李老师我想请教一个问题，最近在多项欧拉Foam里植入了一个曳力模型，编译成功了但是求解动量方程的时候就会报错，我想请问这个报错的意思是我这pow函数没有用对吗？报错如下

PIMPLE: Iteration 1 MULES: Solving for alpha.air air fraction, min, max = 0.293333 0 1 MULES: Solving for alpha.air air fraction, min, max = 0.293333 0 1 Constructing momentum equations #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 ? in "/lib/x86_64-linux-gnu/libm.so.6" #4 powf64 in "/lib/x86_64-linux-gnu/libm.so.6" #5 Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:? #6 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) at ??:? #7 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, double const&) at ??:? #8 Foam::dragModels::bubbleDrag::CdRe() const at ??:? #9 Foam::dragModels::dispersedDragModel::Ki() const at ??:? #10 Foam::dragModels::dispersedDragModel::K() const at ??:? #11 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::BlendedInterfacialModel<Foam::dragModel>::evaluate<double, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > (Foam::dragModel::*)() const, Foam::word const&, Foam::dimensionSet const&, bool) const at ??:? #12 Foam::blendedDragModel::K() const at ??:? #13 Foam::MomentumTransferPhaseSystem<Foam::phaseSystem>::momentumTransfer() at ??:? #14 Foam::PhaseTransferPhaseSystem<Foam::OneResistanceHeatTransferPhaseSystem<Foam::MomentumTransferPhaseSystem<Foam::phaseSystem> > >::momentumTransfer() at ??:? #15 ? in "/home/gzy/OpenFOAM/OpenFOAM-10/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" #16 ? in "/lib/x86_64-linux-gnu/libc.so.6" #17 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #18 ? in "/home/gzy/OpenFOAM/OpenFOAM-10/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam" Floating point exception (core dumped)

我目前的想法是在主函数里输出一下这些场，然后逐个排查下哪里出了问题，请问这些场该怎么调用呢？直接info+场名行不通，我不大清楚这样的语句要怎么写，麻烦您指导一下，非常感谢，例如

Foam::dragModels::bubbleDrag::bubbleDrag ( const dictionary& dict, const phaseInterface& interface, const bool registerObject ) : dispersedDragModel(dict, interface, registerObject), residualRe_("residualRe", dimless, dict) {} // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * // Foam::dragModels::bubbleDrag::~bubbleDrag() {} // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // Foam::tmp<Foam::volScalarField> Foam::dragModels::bubbleDrag::CdRe() const { volScalarField TTT(interface_.continuous().thermo().T());:xiexie: TTT.dimensions().reset(dimless); const volScalarField g1(1-TTT/647); const volScalarField g2(max(g1, residualRe_)); const volScalarField g3(pow(g2, 0.98)); volScalarField ddd(interface_.dispersed().d()); ddd.dimensions().reset(dimless);

:xiexie: