
[Sponsors] 
rhoSimpleFoam High Pressure Cell Crashes Simulation 

LinkBack  Thread Tools  Search this Thread  Display Modes 
March 3, 2023, 05:02 
rhoSimpleFoam High Pressure Cell Crashes Simulation

#1 
New Member
Kuzey Can Derman
Join Date: Oct 2022
Posts: 4
Rep Power: 2 
Hello, I am trying to run a 0.5 Mach flow around a doublewedge on rhoSimpleFoam. After around 1100 steps, random high pressure cell(s) appears on a single cell at the corner. Only one cell is visible and there is no occurence of symmetry of this cell. And after 1200s steps, errors start to rise up and at step 1225 simulation crashes.
Note: Values of some boundary conditions may not correspond to a real condition as this is just a study for me to learn compressible flow simulations. You can download the case from this link if you want (Sorry I could only put one timestep) https://we.tl/tyGfIZvyy96 Boundary Conditions are like this for U,p and T U: Code:
dimensions [0 1 1 0 0 0 0]; internalField uniform (150 0 0); boundaryField { inlet { type fixedValue; value $internalField; } wedge { type noSlip; } atmosphere { type inletOutlet; inletValue $internalField; value $internalField; } } Code:
dimensions [1 1 2 0 0 0 0]; internalField uniform 1.0e5; boundaryField { inlet { type zeroGradient; } wedge { type zeroGradient; } atmosphere { type fixedValue; value uniform 1.0e5; } } Code:
dimensions [0 0 0 1 0 0 0]; internalField uniform 250; boundaryField { inlet { type zeroGradient; //value uniform 250; } wedge { type zeroGradient; } atmosphere { type inletOutlet; inletValue $internalField; value $internalField; } } fvSolution: Code:
solvers { p { solver GAMG; smoother DICGaussSeidel; cacheAgglomeration no; maxIter 50; tolerance 1e5; relTol 0; } pFinal { $p; tolerance 1e6; relTol 0; } "pcorr.*" { $p; tolerance 1e2; relTol 0; } "(Ukomegae)" { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e6; relTol 0; } "(Ukomegae)Final" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e6; relTol 0; } } SIMPLE { //correctPhi yes; //nOuterCorrectors 50; //nCorrectors 2; nNonOrthogonalCorrectors 1; outerCorrectorResidualControl { "(Ukomegap)" { relTol 0; tolerance 1e6; } } } relaxationFactors { fields { p 0.4; rho 0.05; } equations { U 0.5; e 0.5; k 0.4; omega 0.4; } cache { grad(U); } Code:
ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; grad(U) cellLimited Gauss linear 1; grad(k) cellLimited Gauss linear 1; grad(omega) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) bounded Gauss linearUpwind grad(U); div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div(phi,e) bounded Gauss linearUpwind limited; div(phi,K) bounded Gauss linearUpwind limited; div(phi,Ekp) bounded Gauss linearUpwind limited; div(phid,p) Gauss upwind; div((phiinterpolate(rho)),p) bounded Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.33; } interpolationSchemes { default linear; } snGradSchemes { default limited corrected 0.33; } wallDist { method meshWave; } Code:
limitT { type limitTemperature; min 100; max 350; selectionMode all; } Code:
Time = 1211 smoothSolver: Solving for Ux, Initial residual = 0.0002279731, Final residual = 3.782625e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.002013097, Final residual = 3.645934e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.07622737, Final residual = 2.000609e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001246553, Final residual = 2.294851e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.06899724, Final residual = 4.537767e06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.0006720902, Final residual = 3.72458e06, No Iterations 2 time step continuity errors : sum local = 2.503897e05, global = 7.908043e06, cumulative = 2.704127e05 smoothSolver: Solving for omega, Initial residual = 3.659904e06, Final residual = 2.418697e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.225559e05, Final residual = 1.641878e07, No Iterations 4 ExecutionTime = 37.19 s ClockTime = 52 s Time = 1212 smoothSolver: Solving for Ux, Initial residual = 0.0002468802, Final residual = 4.731498e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003135094, Final residual = 5.480615e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.1135995, Final residual = 2.676881e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001361711, Final residual = 2.963579e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0379519, Final residual = 4.56387e06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0006912517, Final residual = 4.099163e06, No Iterations 2 time step continuity errors : sum local = 2.212178e05, global = 2.862645e06, cumulative = 2.990392e05 smoothSolver: Solving for omega, Initial residual = 3.675242e06, Final residual = 2.458315e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.180378e05, Final residual = 1.710549e07, No Iterations 4 ExecutionTime = 40.28 s ClockTime = 56 s Time = 1213 smoothSolver: Solving for Ux, Initial residual = 0.0002241639, Final residual = 4.268223e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003180212, Final residual = 5.78672e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.113732, Final residual = 2.98224e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001200583, Final residual = 2.528982e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.03340179, Final residual = 5.555115e06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.003993981, Final residual = 8.071736e06, No Iterations 3 time step continuity errors : sum local = 1.510507e05, global = 7.377566e06, cumulative = 2.252635e05 smoothSolver: Solving for omega, Initial residual = 3.655159e06, Final residual = 2.476639e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.145421e05, Final residual = 1.733164e07, No Iterations 4 ExecutionTime = 43.28 s ClockTime = 60 s Time = 1214 smoothSolver: Solving for Ux, Initial residual = 0.0002224467, Final residual = 3.594115e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.002194249, Final residual = 4.026306e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.08134264, Final residual = 2.214839e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001169864, Final residual = 1.975844e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.04619271, Final residual = 7.885469e06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.003035966, Final residual = 7.081895e06, No Iterations 3 time step continuity errors : sum local = 1.514269e05, global = 7.148483e06, cumulative = 1.537787e05 smoothSolver: Solving for omega, Initial residual = 3.718127e06, Final residual = 2.549116e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.172644e05, Final residual = 1.798811e07, No Iterations 4 ExecutionTime = 46.49 s ClockTime = 64 s Time = 1215 smoothSolver: Solving for Ux, Initial residual = 0.0002521356, Final residual = 4.383906e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.002684685, Final residual = 4.70982e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.0989898, Final residual = 2.311732e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001274286, Final residual = 2.504468e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.09298018, Final residual = 5.140296e06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.001378192, Final residual = 8.791779e06, No Iterations 2 time step continuity errors : sum local = 1.658667e05, global = 4.826956e06, cumulative = 1.055091e05 smoothSolver: Solving for omega, Initial residual = 3.621697e06, Final residual = 2.507616e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.015491e05, Final residual = 1.793056e07, No Iterations 4 ExecutionTime = 49.8 s ClockTime = 68 s Time = 1216 smoothSolver: Solving for Ux, Initial residual = 0.0002364796, Final residual = 4.730065e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003417261, Final residual = 6.382248e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.1162457, Final residual = 2.941088e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001207433, Final residual = 2.65458e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.2301868, Final residual = 5.026898e06, No Iterations 10 GAMG: Solving for p, Initial residual = 0.003008972, Final residual = 5.638641e06, No Iterations 3 time step continuity errors : sum local = 1.209115e05, global = 4.109161e06, cumulative = 1.466007e05 smoothSolver: Solving for omega, Initial residual = 3.612333e06, Final residual = 2.524231e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.040957e05, Final residual = 1.829401e07, No Iterations 4 ExecutionTime = 52.89 s ClockTime = 73 s Time = 1217 smoothSolver: Solving for Ux, Initial residual = 0.0002197738, Final residual = 3.65705e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.00269095, Final residual = 4.907718e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.0876302, Final residual = 2.269867e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001141772, Final residual = 1.895016e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.1576677, Final residual = 4.94643e06, No Iterations 10 GAMG: Solving for p, Initial residual = 0.001269915, Final residual = 6.67411e06, No Iterations 2 time step continuity errors : sum local = 3.63759e05, global = 9.716975e06, cumulative = 2.437705e05 smoothSolver: Solving for omega, Initial residual = 3.644851e06, Final residual = 2.566505e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.173449e05, Final residual = 1.898109e07, No Iterations 4 ExecutionTime = 55.96 s ClockTime = 77 s Time = 1218 smoothSolver: Solving for Ux, Initial residual = 0.0003425709, Final residual = 4.074055e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003214182, Final residual = 3.976894e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.09876538, Final residual = 1.654451e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.002131877, Final residual = 2.5353e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.05914744, Final residual = 9.914611e06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0005479791, Final residual = 3.606309e06, No Iterations 2 time step continuity errors : sum local = 2.396532e05, global = 8.534779e06, cumulative = 3.291183e05 smoothSolver: Solving for omega, Initial residual = 3.649319e06, Final residual = 2.579692e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.25176e05, Final residual = 1.950485e07, No Iterations 4 ExecutionTime = 59.08 s ClockTime = 81 s Time = 1219 smoothSolver: Solving for Ux, Initial residual = 0.0003253982, Final residual = 4.879725e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.008147821, Final residual = 5.986769e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.1219774, Final residual = 2.230793e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001623936, Final residual = 3.040002e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.04126586, Final residual = 6.367415e06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0009362383, Final residual = 5.65737e06, No Iterations 2 time step continuity errors : sum local = 2.725185e05, global = 4.739751e07, cumulative = 3.243785e05 smoothSolver: Solving for omega, Initial residual = 3.688261e06, Final residual = 2.611056e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.274386e05, Final residual = 1.96346e07, No Iterations 4 ExecutionTime = 62.36 s ClockTime = 86 s Time = 1220 smoothSolver: Solving for Ux, Initial residual = 0.0003331128, Final residual = 4.089699e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.004665316, Final residual = 5.607162e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2331707, Final residual = 1.439117e07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001604169, Final residual = 2.4175e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.07735929, Final residual = 8.633926e06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.005453768, Final residual = 4.782639e06, No Iterations 4 time step continuity errors : sum local = 6.636625e06, global = 2.473585e06, cumulative = 2.996427e05 smoothSolver: Solving for omega, Initial residual = 3.761971e06, Final residual = 2.660331e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.476273e05, Final residual = 1.965373e07, No Iterations 4 ExecutionTime = 65.74 s ClockTime = 90 s Time = 1221 smoothSolver: Solving for Ux, Initial residual = 0.0005524315, Final residual = 3.806285e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.006251112, Final residual = 3.782886e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.09837613, Final residual = 9.192444e07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.002134693, Final residual = 2.034312e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.1362519, Final residual = 4.802647e06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.002245415, Final residual = 5.309531e06, No Iterations 3 time step continuity errors : sum local = 1.33171e05, global = 5.968186e06, cumulative = 2.399608e05 smoothSolver: Solving for omega, Initial residual = 3.82601e06, Final residual = 2.698813e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.587721e05, Final residual = 1.993957e07, No Iterations 4 ExecutionTime = 68.83 s ClockTime = 94 s Time = 1222 smoothSolver: Solving for Ux, Initial residual = 0.002155755, Final residual = 6.634064e07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.02980976, Final residual = 5.285516e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2594953, Final residual = 5.583388e07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.01653358, Final residual = 2.847587e07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.2467346, Final residual = 6.36241e06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.001604496, Final residual = 3.598005e06, No Iterations 3 time step continuity errors : sum local = 6.503624e06, global = 1.469501e06, cumulative = 2.252658e05 smoothSolver: Solving for omega, Initial residual = 3.694153e06, Final residual = 2.595912e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 4.202257e05, Final residual = 1.993223e07, No Iterations 4 ExecutionTime = 72.13 s ClockTime = 98 s Time = 1223 smoothSolver: Solving for Ux, Initial residual = 0.002054502, Final residual = 9.479456e08, No Iterations 8 smoothSolver: Solving for Uy, Initial residual = 0.03078286, Final residual = 5.86827e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2124491, Final residual = 5.636619e07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.4303646, Final residual = 2.595784e07, No Iterations 10 GAMG: Solving for p, Initial residual = 0.6008846, Final residual = 4.838204e06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.001194825, Final residual = 7.152647e06, No Iterations 2 time step continuity errors : sum local = 4.058805e05, global = 1.127372e05, cumulative = 3.38003e05 smoothSolver: Solving for omega, Initial residual = 3.593287e06, Final residual = 2.512341e07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 6.240293e05, Final residual = 2.134858e07, No Iterations 4 ExecutionTime = 75.12 s ClockTime = 103 s Time = 1224 smoothSolver: Solving for Ux, Initial residual = 0.02408419, Final residual = 9.202684e07, No Iterations 8 smoothSolver: Solving for Uy, Initial residual = 0.06743824, Final residual = 6.133449e07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2085229, Final residual = 2.407561e07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.6109017, Final residual = 6.899753e08, No Iterations 12 GAMG: Solving for p, Initial residual = 0.9618585, Final residual = 1.220895e+53, No Iterations 50 GAMG: Solving for p, Initial residual = 0.9792507, Final residual = 6.461137e+53, No Iterations 50 time step continuity errors : sum local = 3.807355e+108, global = 3.241565e+108, cumulative = 3.241565e+108 smoothSolver: Solving for omega, Initial residual = 5.208029e07, Final residual = 5.208029e07, No Iterations 0 smoothSolver: Solving for k, Initial residual = 0.05769598, Final residual = 9.754113e18, No Iterations 2 ExecutionTime = 83.78 s ClockTime = 115 s Time = 1225  Primary job terminated normally, but 1 process returned a nonzero exit code. Per userdirection, the job has been aborted.  Photo from step 1100 showing pressure: Do not hesitate to ask anything. Thanks for any help 

Tags 
crash, openfoam, rhosimplefoam 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
viscosity udf don't use correct temperature and strain rate  rezvani  Fluent UDF and Scheme Programming  8  May 27, 2021 05:40 
High negative pressure in incompressible flow simulation  user_HH  OpenFOAM Programming & Development  0  February 4, 2021 14:40 
Simulation crashes at high tolerance  nandiganavishal  OpenFOAM  0  February 7, 2011 10:40 
Neumann pressure BC and velocity field  Antech  Main CFD Forum  0  April 25, 2006 02:15 
Hydrostatic pressure in 2phase flow modeling (long)  DS & HB  Main CFD Forum  0  January 8, 2000 15:00 