CFD中文网

好的好的谢谢大家的意见~！！！！非常感谢大家，根据大家的意见已经调整好了，最近再准备毕业事宜所以没有及时感谢大家！

在breaking wave领域，这种海浪结构，会将空气卷吸并形成很多小气泡。最近提交一个稿件，审稿人推荐我们讨论一下这方面的研究。我详细看了一下，确实有点意思。

看下图，一些实验研究已经证明，这种海浪拍打过程卷吸进的气泡具有尺寸分布。见过很多文献，都是用的DNS直接模拟研究的，结果也都能证实。

但目前我想的是，能否用PBE模型，或者其他类似模型做相关的研究。

捕获.JPG

气泡直径和分布的关系
捕获.JPG

一些参考文献

Deike, L., Melville, W. K. & Popinet, S. 2016. Air entrainment and bubble statistics in breaking waves. J. Fluid Mech. 801, 91–129.
Wang, Z., Yang, J., Stern, F., 2016. High-fidelity simulations of bubble, droplet and spray formation in breaking waves. J. Fluid Mech. 792, 307–327.
G. Soligo, A. Roccon, and A. Soldati, Breakage, coalescence and size distribution of surfactant-laden droplets in turbulent flow, J. Fluid Mech. 881, 244 (2019).
Ahmed, Z., Izbassarov, D., Costa, P., Muradoglu, M., Tammisola, O. 2020. Turbulent bubbly channel flows: Effects of soluble surfactant and viscoelasticity, Computers & Fluids 212, 104717.

Why the two-fluid model fails to predict the bed expansion characteristics of Geldart A particles in gas-fluidized beds: A tentative answer

http://www.sciencedirect.com/science/article/pii/S0009250908005137

short communication类的，找了几个都是做的生物、实验什么的，看的不明不白，打算研究研究文体和思路。谢谢

太厉害了。。。。。。又做实验又写代码， 全才啊

https://www.esi-group.com/sites/default/files/resource/other/6395/abstract-roe-university-michigan-upwinding-past-present-future.pdf

We are pleased to announce that Prof. Philip Roe, best known for his work on high-resolution compressible flows based on approximate Riemann solvers and Godunov methods will deliver the keynote at this year’s conference.

@veen 大佬搞好没:140:

那就在做O形网格的时候，“O穿”整个区域。就是假设那个圆柱两端无限延伸，对这个无限长的圆柱做O形网格。

如果不知道空心实心怎么选。做的时候；part里有个part，是专门放已经被删除的block；把这个part打开。就当是个全实心的圆柱。完了再关了那个part。

明白了，感谢各位的回复

休斯顿市区内150平米loft卖15万美元：https://www.zillow.com/homes/for_sale/Houston-TX/pmf,pf_pt/39051_rid/globalrelevanceex_sort/29.893043,-95.053368,29.656432,-95.539513_rect/10_zm/

:136: 好吧找到原因了，是这个版本的link edge factor默认是0

组分运输问题，组分守恒方程里除对流项之外，还含有一项是质量分数对x方向的导数乘上一个系数，这一项能否看作是质量源项处理？可以的话怎么编写这个udf？
除此之外还有别的处理方法吗？这是一个对流扩散问题，控制方程为质量守恒和动量守恒。
谢谢大家

@winsway_zero 前辈，请问gmesh怎么划分block

因为分散力模型无法找到，https://zhuanlan.zhihu.com/p/90057423

1.JPG
src文件夹下没有fvoptions，但想利用其中的CellSetOption功能，有人知道吗？

@羽之下 thank you，谢谢你的回复。我记得好像在哪看过说这里引用了一个假设，我决定把这部分系统的看一看，有什么新发现我会继续在这里交流~

@lixin19981013 可以联系我，我教你。qq，1259451345

在DPM中粒子所受力例如重力、虚拟质量力等都是通过forceSuSp中的vector Su以及scalar Sp来决定的，最后通过KinematicParcel.C来计算Fcp以及Fncp(这个地方是每个力单独耦合求解速度还是把所有的力都考虑进来一起耦合),不知道理解的对不对，现在我想要把每个粒子的单独受力比如重力、虚拟质量力等都单独写出来，不知道各位老师有没有什么办法去实现？Ps:我看另一个帖子中写到F=Su*(U-Ud)+Sp.这个是在哪里出现的我也没找到，Ud是啥，可以直接用吗？

这是单独一个力的以重力为例：

template<class CloudType> Foam::forceSuSp Foam::GravityForce<CloudType>::calcNonCoupled ( const typename CloudType::parcelType& p, const scalar dt, const scalar mass, const scalar Re, const scalar muc ) const { forceSuSp value(Zero, 0.0); value.Su() = mass*g_*(1.0 - p.rhoc()/p.rho());//这里求出来Su return value; }

这是KinematicParcel.C出现通过力算速度的：

const forceType& forces = td.cloud().forces(); //这是最后粒子受力耦合求解速度 // Momentum source due to particle forces const parcelType& p = static_cast<const parcelType&>(*this); //这个地方是每个力单独耦合求解速度还是把所有的力都考虑进来一起耦合？？ const forceSuSp Fcp = forces.calcCoupled(p, dt, mass, Re, mu); const forceSuSp Fncp = forces.calcNonCoupled(p, dt, mass, Re, mu); const forceSuSp Feff = Fcp + Fncp; const scalar massEff = forces.massEff(p, mass); // New particle velocity //~~~~~~~~~~~~~~~~~~~~~~ // Update velocity - treat as 3-D const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/massEff; const scalar bp = Feff.Sp()/massEff; Spu = dt*Feff.Sp(); IntegrationScheme<vector>::integrationResult Ures = td.cloud().UIntegrator().integrate(U_, dt, abp, bp); vector Unew = Ures.value();

glgoo，我一直在用的

请教各位老师，外力驱动的流动，常用哪些算例做benchmark，比如Poiseuille流动，非定常的Womersley流动，还有吗？