CFD中文网

@xpqiu 感谢老师详细指点，试了下这个方法确实简洁:xiexie:

@mohui 挺好挺好！燥动起来！ :high:

各位老师好，我想在OF里植入一个与流动直径有关的边界条件。通过sHM后得到的网格，其各流动通道直径不一且比较复杂，能否获取各网格所在流动通道的当量直径呢？在网上查阅wallDist(mesh).y()可以计算网格到壁面的最近的距离，但不知道能否用来计算出直径。

@bestucan 非常感谢老师的详细建议。因为集群之前就已经安装过OF，我是在个人账户里安装OF7版本的，所以可能就差这一个包。我已经使用了跑了计算，目前来看没什么问题:xiexie:

@东岳 谢谢东岳老师！:xiexie:

dimensionedScalar nu ( "nu", dimViscosity, transportProperties ); wordList patchTypes(mesh.boundaryMesh().size(),"fixedValue"); volScalarField nu0 ( IOobject ( "nu0", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh, nu, patchTypes ); const scalar nuWanted(1000.0); labelUList neiCells(nu0.boundaryFieldRef()[0].patch().faceCells()); forAll(neiCells, cellI) { nu0.primitiveFieldRef()[neiCells[cellI]]=nuWanted; }

@东岳 老师，这个公式（1）是他们参考另一篇文章
Wu, K., Chen, Z., Li, J., Li, X., Xu, J., & Dong, X. (2017). Wettability effect on nanoconfined water flow. Proceedings of the National Academy of Sciences, 114(13), 3358-3363. https://www.pnas.org/content/114/13/3358.short
我的理解是因为在纳米管内靠近壁面的位置，粘度变化没有统一的解析解，所以看到一些文章就这么分成“多层”粘度，算是粗略简化的方法吧 :xiezuoye:

@麦迪文 那您觉得应该从哪方面入手，小白我还没有改过边界

@东岳 非常感谢东岳老师回复。请问不成立的原因是什么呢？