另外有没有可能是这类边界只适用于稳态求解器来忽略气体鼓入之后移动到液面这一过程的影响,对于multiphaseEulerFoam这种瞬态求解器而不适用呢?
青山
帖子
-
-
最近在用multiphaseEulerFoam做一个底部喷气的气液固三相搅拌釜,在设置顶部出口边界条件时发现了一个小问题,就是怎么为每个相设置合适的出口边界条件。对于液固体系或者没有气体存在的体系,将顶部出口设置为壁面即可。但是当存在气体鼓入时,必须将顶部设为patch并为气体提供出口,那么这时候应该对液相和固相的出口速度条件设为固定值(0 0 0)以防止流出。通过查阅文献大家也都是按照这个规则设置的。
以及
但是问题来了,当水和分散相都设置为无法流出时,由于初始场中并没有气体,那么在刚开始鼓气的过程中,气体还没接触到液面,也就是气体只进不出,液体和分散相也不流出,这显然是非物理的现象,并且当我使用这种类型的出口条件时,相分数之和会慢慢大于1,这显然也是不对的,所以文献中为什么可以设置成这种非物理的边界条件并且可以正常计算?望大佬解答,感谢! -
会不会是对于这种上层层流,下层湍流的存在较大湍动差异的案例,k-e模型不适用?k-e模型应该使用在完全发展的湍流案例中是吗?
-
目前在做一个气液固三相的搅拌模拟,我想在搅拌釜中先setFields一个有20%固含率和80%液体的液位,上部是空气,反应器底部有气体鼓入,使用openfoam9的multiphaseEulerFoam,气液湍流使用的是k-e模型,固体湍流使用kineticTheory模型。目前遇到的问题是:当我给上部气相设成纯气体的话,直接运算不起来,第一个时间步都跑不了就结束了。只有给上部气体设一个很小的液体含量才可以正常运行,但是这样的话,计算出来结果是液位上的k.water,k.gas,epsilon.water,epsilon.gas都会出现一块很大的值,计算过程中会出现k和epsilon的bounding,然后导致计算发散,如下图
alpha云图看着到是挺正常,
离散格式也都使用了迎风格式,
使用层流的话计算是没有任何问题的,而且如果不设置上层空气,使用液体填满的话,计算也是没问题的,一旦设置自由液位就有问题。而且为了简化网格中间搅拌桨我使用了一块平板代替,网格是很规整的六面体网格,所以应该也不是网格问题,请问各位大佬有没有遇到过这样的问题,望不吝赐教! -
@四季之中的隐星 谢谢,phaseproperties中好像只能开关相变模型和设置饱和蒸气压,具体的是从A相变为B相还是C相好像在alphat中指定,我看你之前做也是闪蒸,请问你的相变的具体相是在哪里设置的?
-
@四季之中的隐星 感谢,那么还有个问题,reactingMuletiPhaseEulerFoam中的发生相变的相和相变之后的相是需要人为指定的,通常在alphat文件中指定,在壁面沸腾案例中使用的是如下边界条件指定liquid相发生相变为vapor相:
那如果我现在不需要热源,我设置初始场温度就是沸点,相变只发生于液体中压力低于饱和蒸气压的部分,那么我的alphat的边界条件是否该设为calculated,那么发生相变的相和发生相变之后的相该在哪里指定? -
打算做一个汽液固三相的搅拌釜,体系温度为液体的沸点,由于搅拌过程中桨叶后方会出现低压区,桨叶后的高温液体会汽化出现气穴,类似空化但是由于温度是沸点所以产生的气泡不会溃灭,如何在openfoam中实现这种由于压力低于饱和蒸气压而发生的相变?因为体系还有固体颗粒,vof可能不太适用,所以打算用reactingMultiPhaseEluerFoam,看案例里都是壁面沸腾导致的液体气化跟我的不太符合,我这种的应该怎样实现?添加pvOption源项可以实现吗,望大佬解答,感谢!
-
@李东岳 好吧,谢谢大佬!
-
如图,我使用ICEM对网格进行了局部的细化,出现了这种一个网格相邻两个网格的现象,而且导入openfoam之后可以正常计算,那么请问像这种一个网格与两个网格相邻时,他们之间的面通量是怎么传输的,根据面积进行差值吗?
-
最近在模拟固体的溶解过程,由于我的固含量较高且需添加反应,因此使用reactingTwoPhaseEulerFoam,但是在openfoam中好像没有设置饱和溶解度来模拟液固传质的模型,只有饱和蒸气压这一说法。而我只有液固两相,如何才能实现饱和溶解度的设置或者说实现固体在液体中的溶解过程?还是只能自己写?望大佬解答!感谢!(openfoam自带的饱和模型在下图,都是计算饱和蒸气压的 )
-
目前正在使用reactingTwoPhaseEulerFoam做一个液固两相模拟,准备将化学反应放进去,已经在chemistryProperties中将化学反应开启,但是发现combustionProperties中必须选择一项combustionModel才能进行反应,如果在combustionModel中选择none,计算是可以计算,但是化学反应也进行不了。那怎么才能只进行化学反应而不使用燃烧模型呢?部分设置如下图,望各位大佬解答!
-
@田畔的风 感谢大佬,说的挺明白的,受教了
-
@李东岳 好的,感谢东岳老师!
-
@李东岳 而且固体会在高速碰撞时发生破碎,但不会聚并
-
@李东岳 感谢东岳老师回复,液-固不可压缩,没有超音速,可以近似认为是圆颗粒
-
最近想搞一个收缩管中的液固颗粒流动,想用欧拉-拉格朗日来模拟,看到openfoam中有DPMFoam和MPPICFoam可以使用,但是我看文献中很多都是OPENFOAM+LIGGGHTS+CFDEM耦合来模拟。请问他们各自的优势在哪里?如何选择适合自己的模拟方法?
-
@李东岳 感谢东岳老师,打算先找售后过来看看,如果他们说硬件没问题,换系统也没有改善的话,年后可能会麻烦东岳老师帮忙 先提前向东岳老师和CFD中文网的朋友们问声新年好!
-
@李东岳 感谢东岳老师,机子在杭州。超微的主板,买来都是配好的,一直没测试,最近测了几个案例才发现问题。标准算例里面8核的并行速度比16核还快我是真想不通
-
@星星星星晴 两个cpu槽插了一个cup,然后在cpu旁边的8个内存插槽都插了,搞不懂啥问题
-
@cccrrryyy 并行是没问题的,就是速度太慢,跟没并行差的不多
-
/*---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 7 \\/ M anipulation | \*---------------------------------------------------------------------------*/ Build : 7-3bcbaf946ae9 Exec : pisoFoam -parallel Date : Jan 18 2022 Time : 10:22:57 Host : "ps" PID : 4116 I/O : uncollated Case : /home/ps/OpenFoam/benchMark/core_16 nProcs : 16 Slaves : 15 ( "ps.4117" "ps.4118" "ps.4119" "ps.4120" "ps.4121" "ps.4122" "ps.4123" "ps.4124" "ps.4125" "ps.4126" "ps.4127" "ps.4128" "ps.4129" "ps.4130" "ps.4131" ) Pstream initialized with: floatTransfer : 0 nProcsSimpleSum : 0 commsType : nonBlocking polling iterations : 0 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No MRF models present Creating finite volume options from "constant/fvOptions" Selecting finite volume options model type meanVelocityForce Source: momentumSource - selecting all cells - selected 2097152 cell(s) with volume 10 Initial pressure gradient = 0 Starting time loop fieldAverage fieldAverage1: Starting averaging at time 0 Time = 0.1 Courant Number mean: 0.731999 max: 1.22083 GAMG: Solving for p, Initial residual = 0.629082, Final residual = 0.00572723, No Iterations 3 time step continuity errors : sum local = 1.06089e-05, global = -1.10517e-20, cumulative = -1.10517e-20 Pressure gradient source: uncorrected Ubar = 0.133489, pressure gradient = 0.000115089 GAMG: Solving for p, Initial residual = 0.0561029, Final residual = 0.000467237, No Iterations 4 time step continuity errors : sum local = 1.28615e-06, global = -3.71636e-20, cumulative = -4.82152e-20 Pressure gradient source: uncorrected Ubar = 0.133489, pressure gradient = 0.000117874 GAMG: Solving for p, Initial residual = 0.0152276, Final residual = 9.92056e-07, No Iterations 30 time step continuity errors : sum local = 2.62628e-09, global = -2.27344e-20, cumulative = -7.09497e-20 Pressure gradient source: uncorrected Ubar = 0.133489, pressure gradient = 0.000117425 ExecutionTime = 5.66 s ClockTime = 6 s Time = 0.2 Courant Number mean: 0.731698 max: 1.1641 GAMG: Solving for p, Initial residual = 0.758295, Final residual = 0.00701504, No Iterations 3 time step continuity errors : sum local = 9.27309e-06, global = -2.6639e-13, cumulative = -2.66391e-13 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119203 GAMG: Solving for p, Initial residual = 0.0736589, Final residual = 0.000696871, No Iterations 6 time step continuity errors : sum local = 6.00649e-07, global = -2.6639e-13, cumulative = -5.32781e-13 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000117205 GAMG: Solving for p, Initial residual = 0.0148478, Final residual = 9.08119e-07, No Iterations 43 time step continuity errors : sum local = 7.91859e-10, global = -2.66391e-13, cumulative = -7.99172e-13 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011676 ExecutionTime = 11.78 s ClockTime = 12 s Time = 0.3 Courant Number mean: 0.731607 max: 1.18248 GAMG: Solving for p, Initial residual = 0.476883, Final residual = 0.00315427, No Iterations 4 time step continuity errors : sum local = 2.37198e-06, global = -3.08998e-13, cumulative = -1.10817e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119187 GAMG: Solving for p, Initial residual = 0.0807666, Final residual = 0.000618685, No Iterations 5 time step continuity errors : sum local = 3.9245e-07, global = -3.08998e-13, cumulative = -1.41717e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000117229 GAMG: Solving for p, Initial residual = 0.018049, Final residual = 8.80851e-07, No Iterations 45 time step continuity errors : sum local = 5.48708e-10, global = -3.08998e-13, cumulative = -1.72617e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011678 ExecutionTime = 18.04 s ClockTime = 18 s Time = 0.4 Courant Number mean: 0.731599 max: 1.17717 GAMG: Solving for p, Initial residual = 0.1167, Final residual = 0.000969819, No Iterations 5 time step continuity errors : sum local = 6.01509e-07, global = -3.34128e-13, cumulative = -2.06029e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011923 GAMG: Solving for p, Initial residual = 0.0809347, Final residual = 0.000631617, No Iterations 5 time step continuity errors : sum local = 3.90613e-07, global = -3.34128e-13, cumulative = -2.39442e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.00011725 GAMG: Solving for p, Initial residual = 0.0184664, Final residual = 6.81898e-07, No Iterations 47 time step continuity errors : sum local = 4.15027e-10, global = -3.34128e-13, cumulative = -2.72855e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000116801 ExecutionTime = 24.52 s ClockTime = 25 s Time = 0.5 Courant Number mean: 0.731615 max: 1.19242 GAMG: Solving for p, Initial residual = 0.108148, Final residual = 0.00095674, No Iterations 5 time step continuity errors : sum local = 5.80503e-07, global = -3.47575e-13, cumulative = -3.07612e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000119271 GAMG: Solving for p, Initial residual = 0.0827827, Final residual = 0.000641763, No Iterations 5 time step continuity errors : sum local = 3.89751e-07, global = -3.47575e-13, cumulative = -3.4237e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000117272 GAMG: Solving for p, Initial residual = 0.0189672, Final residual = 9.67464e-07, No Iterations 42 time step continuity errors : sum local = 5.75885e-10, global = -3.47575e-13, cumulative = -3.77127e-12 Pressure gradient source: uncorrected Ubar = 0.1335, pressure gradient = 0.000116822 ExecutionTime = 30.63 s ClockTime = 31 s
-
------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 7 \\/ M anipulation | \*---------------------------------------------------------------------------*/ Build : 7-3bcbaf946ae9 Exec : reactingTwoPhaseEulerFoam -parallel Date : Jan 18 2022 Time : 13:20:12 Host : "ps" PID : 6349 I/O : uncollated Case : /home/ps/OpenFoam/zf/tube nProcs : 24 Slaves : 23 ( "ps.6350" "ps.6351" "ps.6352" "ps.6353" "ps.6354" "ps.6355" "ps.6356" "ps.6357" "ps.6358" "ps.6359" "ps.6360" "ps.6361" "ps.6362" "ps.6363" "ps.6364" "ps.6365" "ps.6366" "ps.6367" "ps.6368" "ps.6369" "ps.6370" "ps.6371" "ps.6372" ) Pstream initialized with: floatTransfer : 0 nProcsSimpleSum : 0 commsType : nonBlocking polling iterations : 0 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calculations will do 3 corrections PIMPLE: Operating solver in transient mode with 3 outer correctors Reading g Reading hRef Creating phaseSystem Selecting twoPhaseSystem populationBalanceTwoPhaseSystem Selecting phaseModel for air: purePhaseModel Selecting diameterModel for phase air: velocityGroup Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Calculating face flux field phi.air Selecting turbulence model type RAS Selecting RAS turbulence model mixtureKEpsilon RAS { RASModel mixtureKEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 1.92; Cp 0.25; sigmak 1; sigmaEps 1.3; } Selecting phaseModel for water: purePhaseModel Selecting diameterModel for phase water: constant Selecting thermodynamics package { type heRhoThermo; mixture pureMixture; transport const; thermo eConst; equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } Calculating face flux field phi.water Selecting turbulence model type RAS Selecting RAS turbulence model mixtureKEpsilon RAS { RASModel mixtureKEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 1.92; Cp 0.25; sigmak 1; sigmaEps 1.3; } No MRF models present Selecting default blending method: linear Selecting drag blending method: linear Selecting surfaceTensionModel for (air and water): constant Selecting aspectRatioModel for (air in water): constant Selecting aspectRatioModel for (water in air): constant Selecting dragModel for (air and water): segregated Selecting swarmCorrection for (air and water): none Selecting dragModel for (air in water): IshiiZuber Selecting swarmCorrection for (air in water): Tomiyama Selecting dragModel for (water in air): IshiiZuber Selecting swarmCorrection for (water in air): Tomiyama Selecting virtualMassModel for (air in water): constantCoefficient Selecting virtualMassModel for (water in air): constantCoefficient Selecting heatTransferModel for (air in water): RanzMarshall Selecting heatTransferModel for (water in air): RanzMarshall Calculating field g.h Reading field p_rgh Courant Number mean: 4.48994e-05 max: 0.00609375 Starting time loop Courant Number mean: 4.48994e-05 max: 0.00609375 Max Ur Courant Number = 0 Time = 5e-05 PIMPLE: Iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 8.97988e-05 Min(alpha1) = 0 Max(alpha1) = 0.25 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.990339, Final residual = 7.12114e-09, No Iterations 43 PIMPLE: Iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 8.97988e-05 Min(alpha1) = 0 Max(alpha1) = 0.25 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.00172072, Final residual = 8.24851e-09, No Iterations 14 PIMPLE: Iteration 3 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 8.97988e-05 Min(alpha1) = 0 Max(alpha1) = 0.25 populationBalance bubbles: Iteration 1 DILUPBiCGStab: Solving for f0.air.bubbles, Initial residual = 1, Final residual = 1.07449e-10, No Iterations 1 DILUPBiCGStab: Solving for f1.air.bubbles, Initial residual = 1, Final residual = 1.0744e-10, No Iterations 1 DILUPBiCGStab: Solving for f2.air.bubbles, Initial residual = 1, Final residual = 1.07565e-10, No Iterations 1 DILUPBiCGStab: Solving for f3.air.bubbles, Initial residual = 1, Final residual = 1.07601e-10, No Iterations 1 DILUPBiCGStab: Solving for f4.air.bubbles, Initial residual = 1, Final residual = 1.07472e-10, No Iterations 1 DILUPBiCGStab: Solving for f5.air.bubbles, Initial residual = 1, Final residual = 1.07428e-10, No Iterations 1 DILUPBiCGStab: Solving for f6.air.bubbles, Initial residual = 1, Final residual = 1.07231e-10, No Iterations 1 DILUPBiCGStab: Solving for f7.air.bubbles, Initial residual = 1, Final residual = 1.07023e-10, No Iterations 1 DILUPBiCGStab: Solving for f8.air.bubbles, Initial residual = 1, Final residual = 1.06289e-10, No Iterations 1 DILUPBiCGStab: Solving for f9.air.bubbles, Initial residual = 1, Final residual = 1.04711e-10, No Iterations 1 DILUPBiCGStab: Solving for f10.air.bubbles, Initial residual = 1, Final residual = 1.00915e-10, No Iterations 1 DILUPBiCGStab: Solving for f11.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for f12.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for f13.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for f14.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 air Sauter mean diameter, min, max = 4.65318e-05 4.65237e-05 4.65318e-05 air sizeGroups-sum volume fraction, min, max = 1 1 1 bubbles sizeGroup phase fraction first, last = 8.85697e-08 0 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.000995987, Final residual = 6.44985e-09, No Iterations 11 smoothSolver: Solving for epsilonm, Initial residual = 0.360485, Final residual = 5.98015e-08, No Iterations 6 smoothSolver: Solving for km, Initial residual = 1, Final residual = 4.14411e-08, No Iterations 6 ExecutionTime = 0.14 s Courant Number mean: 0.0537623 max: 1.16043 Max Ur Courant Number = 0.0897572 Time = 0.0001 PIMPLE: Iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.000179598 Min(alpha1) = 0 Max(alpha1) = 0.25 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.00793744, Final residual = 8.90127e-09, No Iterations 30 PIMPLE: Iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.000179598 Min(alpha1) = 0 Max(alpha1) = 0.25 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.0765251, Final residual = 7.71052e-09, No Iterations 16 PIMPLE: Iteration 3 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.000179598 Min(alpha1) = 0 Max(alpha1) = 0.25 populationBalance bubbles: Iteration 1 DILUPBiCGStab: Solving for f0.air.bubbles, Initial residual = 0.147551, Final residual = 4.86871e-11, No Iterations 1 DILUPBiCGStab: Solving for f1.air.bubbles, Initial residual = 0.0394195, Final residual = 1.02867e-11, No Iterations 1 DILUPBiCGStab: Solving for f2.air.bubbles, Initial residual = 0.736843, Final residual = 2.25619e-10, No Iterations 1 DILUPBiCGStab: Solving for f3.air.bubbles, Initial residual = 0.793873, Final residual = 2.42534e-10, No Iterations 1 DILUPBiCGStab: Solving for f4.air.bubbles, Initial residual = 0.857026, Final residual = 2.6112e-10, No Iterations 1 DILUPBiCGStab: Solving for f5.air.bubbles, Initial residual = 0.608674, Final residual = 1.86077e-10, No Iterations 1 DILUPBiCGStab: Solving for f6.air.bubbles, Initial residual = 0.196465, Final residual = 5.34906e-11, No Iterations 1 DILUPBiCGStab: Solving for f7.air.bubbles, Initial residual = 0.596202, Final residual = 1.77238e-10, No Iterations 1 DILUPBiCGStab: Solving for f8.air.bubbles, Initial residual = 0.32325, Final residual = 9.33873e-11, No Iterations 1 DILUPBiCGStab: Solving for f9.air.bubbles, Initial residual = 0.785256, Final residual = 1.98355e-10, No Iterations 1 DILUPBiCGStab: Solving for f10.air.bubbles, Initial residual = 0.727823, Final residual = 1.4766e-10, No Iterations 1 DILUPBiCGStab: Solving for f11.air.bubbles, Initial residual = 1, Final residual = 1.35917e-10, No Iterations 1 DILUPBiCGStab: Solving for f12.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for f13.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for f14.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 air Sauter mean diameter, min, max = 4.6532e-05 4.65318e-05 4.65537e-05 air sizeGroups-sum volume fraction, min, max = 1 0.999972 1 bubbles sizeGroup phase fraction first, last = 1.77403e-07 0 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.0649563, Final residual = 6.62621e-09, No Iterations 17 smoothSolver: Solving for epsilonm, Initial residual = 0.24891, Final residual = 1.7812e-08, No Iterations 6 smoothSolver: Solving for km, Initial residual = 0.212581, Final residual = 2.92726e-08, No Iterations 6 ExecutionTime = 0.27 s Courant Number mean: 0.0538033 max: 1.18845 Max Ur Courant Number = 0.12208 Time = 0.00015 PIMPLE: Iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.000269397 Min(alpha1) = 0 Max(alpha1) = 0.25 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.148212, Final residual = 7.88969e-09, No Iterations 17 PIMPLE: Iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.000269397 Min(alpha1) = 0 Max(alpha1) = 0.25 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.0744514, Final residual = 5.55115e-09, No Iterations 18 PIMPLE: Iteration 3 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.000269397 Min(alpha1) = 0 Max(alpha1) = 0.25 populationBalance bubbles: Iteration 1 DILUPBiCGStab: Solving for f0.air.bubbles, Initial residual = 0.0782433, Final residual = 4.79582e-11, No Iterations 1 DILUPBiCGStab: Solving for f1.air.bubbles, Initial residual = 0.163708, Final residual = 8.79559e-11, No Iterations 1 DILUPBiCGStab: Solving for f2.air.bubbles, Initial residual = 0.412078, Final residual = 2.36687e-10, No Iterations 1 DILUPBiCGStab: Solving for f3.air.bubbles, Initial residual = 0.425663, Final residual = 2.43918e-10, No Iterations 1 DILUPBiCGStab: Solving for f4.air.bubbles, Initial residual = 0.43729, Final residual = 2.49922e-10, No Iterations 1 DILUPBiCGStab: Solving for f5.air.bubbles, Initial residual = 0.376317, Final residual = 2.14072e-10, No Iterations 1 DILUPBiCGStab: Solving for f6.air.bubbles, Initial residual = 0.974911, Final residual = 5.33381e-10, No Iterations 1 DILUPBiCGStab: Solving for f7.air.bubbles, Initial residual = 0.367985, Final residual = 1.98277e-10, No Iterations 1 DILUPBiCGStab: Solving for f8.air.bubbles, Initial residual = 0.238507, Final residual = 1.18735e-10, No Iterations 1 DILUPBiCGStab: Solving for f9.air.bubbles, Initial residual = 0.396462, Final residual = 1.61137e-10, No Iterations 1 DILUPBiCGStab: Solving for f10.air.bubbles, Initial residual = 0.365836, Final residual = 1.05167e-10, No Iterations 1 DILUPBiCGStab: Solving for f11.air.bubbles, Initial residual = 0.798764, Final residual = 1.27828e-10, No Iterations 1 DILUPBiCGStab: Solving for f12.air.bubbles, Initial residual = 0.999974, Final residual = 7.83893e-11, No Iterations 1 DILUPBiCGStab: Solving for f13.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for f14.air.bubbles, Initial residual = 0, Final residual = 0, No Iterations 0 air Sauter mean diameter, min, max = 4.65323e-05 4.65318e-05 4.65886e-05 air sizeGroups-sum volume fraction, min, max = 0.999999 0.999917 1 bubbles sizeGroup phase fraction first, last = 2.66368e-07 0 Constructing momentum equations GAMG: Solving for p_rgh, Initial residual = 0.0789003, Final residual = 8.16298e-09, No Iterations 16 smoothSolver: Solving for epsilonm, Initial residual = 0.217737, Final residual = 1.58357e-08, No Iterations 6 smoothSolver: Solving for km, Initial residual = 0.191623, Final residual = 2.29973e-08, No Iterations 8 ExecutionTime = 0.38 s
-
@星星星星晴 好的吧
-
@星星星星晴 这里有问题吗?
-
@星星星星晴 来了,大佬帮我看看,内存是16*8的
*---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 7 \\/ M anipulation | \*---------------------------------------------------------------------------*/ Build : 7-3bcbaf946ae9 Exec : reactingTwoPhaseEulerFoam -parallel Date : Jan 18 2022 Time : 13:20:12 Host : "ps" PID : 6349 I/O : uncollated Case : /home/ps/OpenFoam/zf/tube nProcs : 24 Slaves : 23 ( "ps.6350" "ps.6351" "ps.6352" "ps.6353" "ps.6354" "ps.6355" "ps.6356" "ps.6357" "ps.6358" "ps.6359" "ps.6360" "ps.6361" "ps.6362" "ps.6363" "ps.6364" "ps.6365" "ps.6366" "ps.6367" "ps.6368" "ps.6369" "ps.6370" "ps.6371" "ps.6372" ) Pstream initialized with: floatTransfer : 0 nProcsSimpleSum : 0 commsType : nonBlocking polling iterations : 0 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10) allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calculations will do 3 corrections PIMPLE: Operating solver in transient mode with 3 outer correctors
-
@cccrrryyy 我看东岳老师测试的,7502跑32核60多秒就跑完了,我比他慢了十倍 不知道是硬件的问题还是设置的问题
-
@李东岳 大佬跑完了,这7502感觉有点拉胯啊
-
@星星星星晴 没问题,明天算一下再回来反馈
-
@星星星星晴 是的,原生的20.04
-
@cccrrryyy 是的,跑起来的时候24个cup的占用都起来了,就是速度没块,我再研究研究,感谢!
-
@cccrrryyy mpirun -np 24 twoPhaseEulerFoam -parallel ,命令应该没问题,我在其他电脑上并行的话都挺快的,就这机子不行。没有超线程,32核64线程我只用了24个核。
-
组里刚弄了台主机,32核64线程的AMD EPYC 7502,ubuntu18.04上跑OF7,但是测试了一下,12核并行、24核并行和不并行的速度都是一样的,几十万网格和几万网格案例都试了,都是并行速度和没并行一样,这是什么原因?跑的时候我也看了系统监视器,的确并行的核的占用都起来了,风扇声音也贼大,但是速度一点没提升是哪里出了问题?
-
@李东岳 好的,感谢东岳老师!这是我们自己组内写的小文档,看来还有很多要完善的地方!
-
新手一枚!请问在openfoam中,时间项的欧拉隐形离散是否真的不受限于库朗数?为什么我使用twoPhaseEulerFoam时对时间项进行欧拉离散,但是库朗数一大还是会发散?这对我的案例造成了严重限制,案例是一个管道中孔径突然缩小,所以小孔中流速骤增,而且网格很小很密,这样库朗数自然会很大,自动调节步长会到1e-7s,计算的就太慢了,有什么办法可以解决这类问题吗?感谢大佬!
搅拌罐的出口边界设置问题
搅拌罐的出口边界设置问题
气液固三相模拟k-epsilon发散
气液固三相模拟k-epsilon发散
reactingMultiPhaseEluerFoam如何实现压力引起的相变?
reactingMultiPhaseEluerFoam如何实现压力引起的相变?
reactingMultiPhaseEluerFoam如何实现压力引起的相变?
网格加密后的通量传输问题
网格加密后的通量传输问题
固体饱和溶解度的设置问题
如何关闭燃烧模型,只开启化学反应
DPMFoam、MPPICFoam以及CFDEM怎么选择?
DPMFoam、MPPICFoam以及CFDEM怎么选择?
DPMFoam、MPPICFoam以及CFDEM怎么选择?
DPMFoam、MPPICFoam以及CFDEM怎么选择?
DPMFoam、MPPICFoam以及CFDEM怎么选择?
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
主机并行运算问题
时间项离散问题
时间项离散问题