# P Divergence for k-Epsilon Model (Simple Foam)

 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 6, 2021, 21:05 P Divergence for k-Epsilon Model (Simple Foam) #1 New Member     UncleRodger Join Date: Jan 2021 Location: Calgary, AB Posts: 1 Rep Power: 0 Hello, I am trying to simulate water flow inside an emitter (drip irrigation) using RANS and LES (Samagorin–Lilly), I have constructed an structured quad mesh in ANSA (https://imgur.com/a/4tp2tA5) it is a very coarse mesh of 25600 volume elements. I am just trying to perform RANS as first step using k-Epsilon model. I reused the pitzDaily simpleFoam tutorial to try and run it using my mesh, I changed everything that was needed (PIDs match to those ones in the tutorial). But for some reason my simulation explodes and diverges, specially p. I left the initial conditions as they are, this is a preliminary simulation, will change them later. Here is the last iteration: Time = 5 smoothSolver: Solving for Ux, Initial residual = 0.579153, Final residual = 0.0324501, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.530541, Final residual = 0.0388022, No Iterations 2 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0258436, No Iterations 2 time step continuity errors : sum local = 9.49971e+40, global = -4.51565e+40, cumulative = -4.51565e+40 smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0233028, No Iterations 2 bounding epsilon, min: -9.60072e+58 max: 2.86832e+60 average: 1.98086e+57 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0699013, No Iterations 1 bounding k, min: -5.66487e+50 max: 1.4514e+52 average: 3.68843e+48 ExecutionTime = 0.47 s ClockTime = 0 s Time = 6 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.086456, No Iterations 8 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0864664, No Iterations 8 *Code Stops* What could be the possible causes of my divergence? I have a suspicion it could be my mesh but not 100% sure, I tried another mesh that was for another school assignment in the same file and it did converge! I am very confused at the moment Any comment would be helpful, thanks!

 Tags cfd, k-epsilon, openfoam, rans