似乎第一步就发散了，是我的初边条设置的有问题吗？
-
是基于tutorial里面的waterChannel改的算例，模型是二维的一个溢洪道，边界条件就在原先的基础上改了一点：
//alpha.water dimensions [0 0 0 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value uniform 1; } walls { type zeroGradient; } outlet { type zeroGradient; value uniform 0; } atmosphere { type inletOutlet; inletValue uniform 0; value uniform 0; } BaseAndTop { type empty; } }//k dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.0001; boundaryField { inlet { type fixedValue; intensity 0.05; value $internalField; } walls { type kqRWallFunction; value $internalField; } outlet { type inletOutlet; value $internalField; inletValue uniform 0; } BaseAndTop { type empty; } ".*" { type inletOutlet; inletValue $internalField; value $internalField; } }//nut dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { walls { type nutkWallFunction; value uniform 0; } BaseAndTop { type empty; } ".*" { type calculated; value uniform 0; } }//omega dimensions [0 0 -1 0 0 0 0]; internalField uniform 0.003; boundaryField { inlet { type fixedValue; value $internalField; } walls { type omegaWallFunction; value $internalField; } outlet { type inletOutlet; value $internalField; inletValue uniform 0; } BaseAndTop { type empty; } ".*" { type inletOutlet; inletValue $internalField; value $internalField; } }//p_rgh dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { atmosphere { type totalPressure; p0 uniform 0; } inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } BaseAndTop { type empty; } walls { type fixedFluxPressure; value uniform 0; } }//s dimensions [0 0 0 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value uniform 0; } walls { type zeroGradient; } outlet { type zeroGradient; value uniform 0; } atmosphere { type inletOutlet; inletValue uniform 0; value uniform 0; } BaseAndTop { type empty; } }//U dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type flowRateInletVelocity; volumetricFlowRate constant 50; } walls { type noSlip; } atmosphere { type pressureInletOutletVelocity; value uniform (0 0 0); } outlet { type inletOutlet; inletValue uniform (0 0 0); value $internalField; } BaseAndTop { type empty; } }第一个迭代是个情况：
Starting time loop Courant Number mean: 0.000528922 max: 0.136511 Interface Courant Number mean: 0 max: 0 deltaT = 0.0120048 Time = 0.0120048 PIMPLE: iteration 1 smoothSolver: Solving for alpha.water, Initial residual = 1, Final residual = 0.000281706, No Iterations 1000 Phase-1 volume fraction = 0.00125164 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water Phase-1 volume fraction = 0.00125164 Min(alpha.water) = -2.54212e-24 Max(alpha.water) = 1 GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 2.7816e+83, No Iterations 50 time step continuity errors : sum local = 3.5015e+80, global = 2.12903e+77, cumulative = 2.12903e+77请问大佬们，这是边界条件的问题吗？
