@李东岳 不开网格自适应的正在尝试，晚点看看效果。
U
unravelkaka
@unravelkaka
帖子
-
单个液滴在横流中的破碎 -
单个液滴在横流中的破碎@李东岳 是fvSolution中的通量修正吗？这个是开启了的。
-
单个液滴在横流中的破碎@李东岳 老师AMI是什么，是指AMR吗？另外dynamicMeshDict中的phi的修正是否也要开启？
correctFluxes ( (phi none) (rhoPhi none) (alphaPhi10 none) (rhoCpLPhi none) (rhoCpVPhi none) (alphaPhi_ none) (alphaPhi.water none) (U none) (interfaceNormal.water none) (interfaceCentre.water none) (nHat none) (nHatf none) (nHatf_ none) (dVf_ none) (Kf_ none) (ghf none) (deltaFunction_ none) (talphaPhiCorr0_ none) ); -
单个液滴在横流中的破碎@李东岳 是的李老师，的确开启了correctPhi。不过我对correctPhi的认识也不够深，记得好像是修正通量确保质量守恒的功能，所以我一般都默认开启了。是开启这一项导致的原因吗？
-
单个液滴在横流中的破碎@李东岳 这次将速度改回30m/s，重新计算了一段时间。
计算1 D/U，2 D/U，3 D/U后，相分布如下。背风侧出现凸起，迎风侧未出现凹陷，与预期的袋装破碎不符。其实不论何种破碎模式，首先会出现一个变形阶段。液滴在前后滞止压的作用下，会逐渐在流向上变薄，向椭球形转变。但计算结果中没有看到这样的趋势。
最大压力出现在背风侧凸起处，而非迎风侧滞止点。
参考前人工作，当前算例应处于袋装破碎。
-
单个液滴在横流中的破碎液滴为水滴，直径2mm，横流流速60m/s时，We大约是103；之前也试过30m/s，We大约是25，是更典型的袋装破碎区，但实际效果也和上图差不多
-
单个液滴在横流中的破碎@李东岳 不好意思李老师，忘记贴出结果了。
这是计算大约3D/U后的相分布图，想不明白背风侧为什么会是这样的形态？按道理这个时候应该已经有鼓袋发展了，或者在受到迎风滞止点的压力，液滴应该在Y方向上延申Z方向上缩小，呈现沿流向被压缩的形态。
-
单个液滴在横流中的破碎@李东岳 李老师，打扰了，能否抽空帮我看看这个小case
-
单个液滴在横流中的破碎实在想不明白问题在哪，各位同学老师能否抽空帮忙看一下
-
单个液滴在横流中的破碎目前使用interFoam想要模拟一下单个液滴的袋装破碎，但总是达不到预期效果，不知道问题出在哪。
以下是边界条件设置，/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: plus | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 60); boundaryField { inlet { type fixedValue; value uniform (0 0 60); } outlet { type pressureInletOutletVelocity; value uniform (0 0 0); } walls { type slip; } } /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: plus | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedFluxPressure; gradient uniform 0; value uniform 0; } outlet { type prghPressure; p uniform 0; value uniform 0; } walls { type zeroGradient; } } /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: plus | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object alpha.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value uniform 0; } outlet { type inletOutlet; inletValue uniform 0; value uniform 0; } walls { type zeroGradient; } } // ************************************************************************* //以下是相场初始化情况，
-
怎样添加新的equationOfState？另外我在youtube上找到一个较为详细的教程，需要的朋友可以看看。
https://www.youtube.com/watch?v=jJHOeDZJUrU -
怎样添加新的equationOfState？@李东岳 感谢老师指点，确实是个很巧妙的方法。
-
怎样添加新的equationOfState？@李东岳 老师，实在打扰，是否有时间帮我看看这个小问题？
-
怎样添加新的equationOfState？前期在of v1812版本中开发了某求解器，现在发现该版本的EOS中并没有rPolynominal，只有以下几类（看到的rPolynominal是我后期添加的），
~/OpenFOAM/v1812/OpenFOAM-v1812/src/thermophysicalModels/specie/equationOfState$ ls Boussinesq adiabaticPerfectFluid incompressiblePerfectGas perfectFluid rPolynomial PengRobinsonGas icoPolynomial linear perfectGas rhoConst于是我从v2206中复制了rPolynominal，当我准备编译时，却没有找到Make文件，在上一级文件夹~/OpenFOAM/v1812/OpenFOAM-v1812/src/thermophysicalModels/specie中有Make/file如下，
atomicWeights/atomicWeights.C specie/specie.C reaction/reactions/makeReactions.C reaction/reactions/makeLangmuirHinshelwoodReactions.C LIB = $(FOAM_LIBBIN)/libspecie这些好像并不涉及EOS的编译，所以我不明白要添加一个EOS应该怎么操作呢？有没有大佬指教一下
-
compressibleInterFoam求解器怎么加入空化模型前辈你好，我最经也遇到这个问题，请问有时间交流一下吗？
-
有关多孔介质求解器中阻滞系数不均匀该怎么处理@李东岳 好的，谢谢老师！
-
有关多孔介质求解器中阻滞系数不均匀该怎么处理@李东岳 老师，我接触openfoam没多久，理解还没有那么深入。我初步是这样想：首先去读取csv文件，生成一个全域的foamfield场文件；然后在createPorousZones.H文件里加载阻滞系数场，并赋到对应的区块上；最后动量方程里有关多孔介质的源项也改成那个全域的阻滞系数场。老师你看看按这样的思路做有没有问题
-
有关多孔介质求解器中阻滞系数不均匀该怎么处理@李东岳 老师，是否有时间给点建议？或者openfoam现在有自带的赋值非均匀多孔介质的方法吗？
-
有关多孔介质求解器中阻滞系数不均匀该怎么处理这是porosityProperties文件和toposet文件有关区块porosity1的内容，
porosity1 //烟丝 { type DarcyForchheimer; cellZone porousz1; D1 1.17e9; d ($D1 $D1 $D1); f (0 0 0); coordinateSystem { type cartesian; origin (0 0 0); rotation { type axes; e1 (1 0 0); e2 (0 1 0); } } } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // actions ( { name porous1; type cellSet; action new; source cylinderToCell; sourceInfo { p1 ( 0 0 0); p2 ( 0 0 0.062); radius 0.00385; } } { name porousz1; type cellZoneSet; action new; source setToCellZone; sourceInfo { set porous1; } }可见区块1是一个半径0.00385，长0.062的圆柱体，整个区域的阻滞系数都是1.17e9，现在我想将阻滞系数设置为沿轴向z非均匀分布的形式，可以通过外部文件（如txt，csv等）读取阻滞系数参数，然后赋值在轴向的不同位置上。原来的DarcyForchheimer是不是并不支持这种赋值方式，是不是需要修改源代码才能实现呢？