|
[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: 12
Rep Power: 3 |
Hello, I am trying to run a 0.5 Mach flow around a double-wedge 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/t-yGfIZvyy96 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 1e-5; relTol 0; } pFinal { $p; tolerance 1e-6; relTol 0; } "pcorr.*" { $p; tolerance 1e-2; relTol 0; } "(U|k|omega|e)" { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e-6; relTol 0; } "(U|k|omega|e)Final" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-6; relTol 0; } } SIMPLE { //correctPhi yes; //nOuterCorrectors 50; //nCorrectors 2; nNonOrthogonalCorrectors 1; outerCorrectorResidualControl { "(U|k|omega|p)" { relTol 0; tolerance 1e-6; } } } 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((phi|interpolate(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.782625e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.002013097, Final residual = 3.645934e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.07622737, Final residual = 2.000609e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001246553, Final residual = 2.294851e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.06899724, Final residual = 4.537767e-06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.0006720902, Final residual = 3.72458e-06, No Iterations 2 time step continuity errors : sum local = 2.503897e-05, global = 7.908043e-06, cumulative = 2.704127e-05 smoothSolver: Solving for omega, Initial residual = 3.659904e-06, Final residual = 2.418697e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.225559e-05, Final residual = 1.641878e-07, No Iterations 4 ExecutionTime = 37.19 s ClockTime = 52 s Time = 1212 smoothSolver: Solving for Ux, Initial residual = 0.0002468802, Final residual = 4.731498e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003135094, Final residual = 5.480615e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.1135995, Final residual = 2.676881e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001361711, Final residual = 2.963579e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0379519, Final residual = 4.56387e-06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0006912517, Final residual = 4.099163e-06, No Iterations 2 time step continuity errors : sum local = 2.212178e-05, global = 2.862645e-06, cumulative = 2.990392e-05 smoothSolver: Solving for omega, Initial residual = 3.675242e-06, Final residual = 2.458315e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.180378e-05, Final residual = 1.710549e-07, No Iterations 4 ExecutionTime = 40.28 s ClockTime = 56 s Time = 1213 smoothSolver: Solving for Ux, Initial residual = 0.0002241639, Final residual = 4.268223e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003180212, Final residual = 5.78672e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.113732, Final residual = 2.98224e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001200583, Final residual = 2.528982e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.03340179, Final residual = 5.555115e-06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.003993981, Final residual = 8.071736e-06, No Iterations 3 time step continuity errors : sum local = 1.510507e-05, global = -7.377566e-06, cumulative = 2.252635e-05 smoothSolver: Solving for omega, Initial residual = 3.655159e-06, Final residual = 2.476639e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.145421e-05, Final residual = 1.733164e-07, No Iterations 4 ExecutionTime = 43.28 s ClockTime = 60 s Time = 1214 smoothSolver: Solving for Ux, Initial residual = 0.0002224467, Final residual = 3.594115e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.002194249, Final residual = 4.026306e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.08134264, Final residual = 2.214839e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001169864, Final residual = 1.975844e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.04619271, Final residual = 7.885469e-06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.003035966, Final residual = 7.081895e-06, No Iterations 3 time step continuity errors : sum local = 1.514269e-05, global = -7.148483e-06, cumulative = 1.537787e-05 smoothSolver: Solving for omega, Initial residual = 3.718127e-06, Final residual = 2.549116e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.172644e-05, Final residual = 1.798811e-07, No Iterations 4 ExecutionTime = 46.49 s ClockTime = 64 s Time = 1215 smoothSolver: Solving for Ux, Initial residual = 0.0002521356, Final residual = 4.383906e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.002684685, Final residual = 4.70982e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.0989898, Final residual = 2.311732e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001274286, Final residual = 2.504468e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.09298018, Final residual = 5.140296e-06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.001378192, Final residual = 8.791779e-06, No Iterations 2 time step continuity errors : sum local = 1.658667e-05, global = -4.826956e-06, cumulative = 1.055091e-05 smoothSolver: Solving for omega, Initial residual = 3.621697e-06, Final residual = 2.507616e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.015491e-05, Final residual = 1.793056e-07, No Iterations 4 ExecutionTime = 49.8 s ClockTime = 68 s Time = 1216 smoothSolver: Solving for Ux, Initial residual = 0.0002364796, Final residual = 4.730065e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003417261, Final residual = 6.382248e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.1162457, Final residual = 2.941088e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001207433, Final residual = 2.65458e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.2301868, Final residual = 5.026898e-06, No Iterations 10 GAMG: Solving for p, Initial residual = 0.003008972, Final residual = 5.638641e-06, No Iterations 3 time step continuity errors : sum local = 1.209115e-05, global = 4.109161e-06, cumulative = 1.466007e-05 smoothSolver: Solving for omega, Initial residual = 3.612333e-06, Final residual = 2.524231e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.040957e-05, Final residual = 1.829401e-07, No Iterations 4 ExecutionTime = 52.89 s ClockTime = 73 s Time = 1217 smoothSolver: Solving for Ux, Initial residual = 0.0002197738, Final residual = 3.65705e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.00269095, Final residual = 4.907718e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.0876302, Final residual = 2.269867e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001141772, Final residual = 1.895016e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.1576677, Final residual = 4.94643e-06, No Iterations 10 GAMG: Solving for p, Initial residual = 0.001269915, Final residual = 6.67411e-06, No Iterations 2 time step continuity errors : sum local = 3.63759e-05, global = 9.716975e-06, cumulative = 2.437705e-05 smoothSolver: Solving for omega, Initial residual = 3.644851e-06, Final residual = 2.566505e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.173449e-05, Final residual = 1.898109e-07, No Iterations 4 ExecutionTime = 55.96 s ClockTime = 77 s Time = 1218 smoothSolver: Solving for Ux, Initial residual = 0.0003425709, Final residual = 4.074055e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.003214182, Final residual = 3.976894e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.09876538, Final residual = 1.654451e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.002131877, Final residual = 2.5353e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.05914744, Final residual = 9.914611e-06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0005479791, Final residual = 3.606309e-06, No Iterations 2 time step continuity errors : sum local = 2.396532e-05, global = 8.534779e-06, cumulative = 3.291183e-05 smoothSolver: Solving for omega, Initial residual = 3.649319e-06, Final residual = 2.579692e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.25176e-05, Final residual = 1.950485e-07, No Iterations 4 ExecutionTime = 59.08 s ClockTime = 81 s Time = 1219 smoothSolver: Solving for Ux, Initial residual = 0.0003253982, Final residual = 4.879725e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.008147821, Final residual = 5.986769e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.1219774, Final residual = 2.230793e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001623936, Final residual = 3.040002e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.04126586, Final residual = 6.367415e-06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.0009362383, Final residual = 5.65737e-06, No Iterations 2 time step continuity errors : sum local = 2.725185e-05, global = -4.739751e-07, cumulative = 3.243785e-05 smoothSolver: Solving for omega, Initial residual = 3.688261e-06, Final residual = 2.611056e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.274386e-05, Final residual = 1.96346e-07, No Iterations 4 ExecutionTime = 62.36 s ClockTime = 86 s Time = 1220 smoothSolver: Solving for Ux, Initial residual = 0.0003331128, Final residual = 4.089699e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.004665316, Final residual = 5.607162e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2331707, Final residual = 1.439117e-07, No Iterations 12 smoothSolver: Solving for e, Initial residual = 0.001604169, Final residual = 2.4175e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.07735929, Final residual = 8.633926e-06, No Iterations 8 GAMG: Solving for p, Initial residual = 0.005453768, Final residual = 4.782639e-06, No Iterations 4 time step continuity errors : sum local = 6.636625e-06, global = -2.473585e-06, cumulative = 2.996427e-05 smoothSolver: Solving for omega, Initial residual = 3.761971e-06, Final residual = 2.660331e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.476273e-05, Final residual = 1.965373e-07, No Iterations 4 ExecutionTime = 65.74 s ClockTime = 90 s Time = 1221 smoothSolver: Solving for Ux, Initial residual = 0.0005524315, Final residual = 3.806285e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.006251112, Final residual = 3.782886e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.09837613, Final residual = 9.192444e-07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.002134693, Final residual = 2.034312e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.1362519, Final residual = 4.802647e-06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.002245415, Final residual = 5.309531e-06, No Iterations 3 time step continuity errors : sum local = 1.33171e-05, global = -5.968186e-06, cumulative = 2.399608e-05 smoothSolver: Solving for omega, Initial residual = 3.82601e-06, Final residual = 2.698813e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 3.587721e-05, Final residual = 1.993957e-07, No Iterations 4 ExecutionTime = 68.83 s ClockTime = 94 s Time = 1222 smoothSolver: Solving for Ux, Initial residual = 0.002155755, Final residual = 6.634064e-07, No Iterations 6 smoothSolver: Solving for Uy, Initial residual = 0.02980976, Final residual = 5.285516e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2594953, Final residual = 5.583388e-07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.01653358, Final residual = 2.847587e-07, No Iterations 8 GAMG: Solving for p, Initial residual = 0.2467346, Final residual = 6.36241e-06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.001604496, Final residual = 3.598005e-06, No Iterations 3 time step continuity errors : sum local = 6.503624e-06, global = -1.469501e-06, cumulative = 2.252658e-05 smoothSolver: Solving for omega, Initial residual = 3.694153e-06, Final residual = 2.595912e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 4.202257e-05, Final residual = 1.993223e-07, No Iterations 4 ExecutionTime = 72.13 s ClockTime = 98 s Time = 1223 smoothSolver: Solving for Ux, Initial residual = 0.002054502, Final residual = 9.479456e-08, No Iterations 8 smoothSolver: Solving for Uy, Initial residual = 0.03078286, Final residual = 5.86827e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2124491, Final residual = 5.636619e-07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.4303646, Final residual = 2.595784e-07, No Iterations 10 GAMG: Solving for p, Initial residual = 0.6008846, Final residual = 4.838204e-06, No Iterations 9 GAMG: Solving for p, Initial residual = 0.001194825, Final residual = 7.152647e-06, No Iterations 2 time step continuity errors : sum local = 4.058805e-05, global = 1.127372e-05, cumulative = 3.38003e-05 smoothSolver: Solving for omega, Initial residual = 3.593287e-06, Final residual = 2.512341e-07, No Iterations 2 smoothSolver: Solving for k, Initial residual = 6.240293e-05, Final residual = 2.134858e-07, No Iterations 4 ExecutionTime = 75.12 s ClockTime = 103 s Time = 1224 smoothSolver: Solving for Ux, Initial residual = 0.02408419, Final residual = 9.202684e-07, No Iterations 8 smoothSolver: Solving for Uy, Initial residual = 0.06743824, Final residual = 6.133449e-07, No Iterations 8 smoothSolver: Solving for Uz, Initial residual = 0.2085229, Final residual = 2.407561e-07, No Iterations 10 smoothSolver: Solving for e, Initial residual = 0.6109017, Final residual = 6.899753e-08, 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.208029e-07, Final residual = 5.208029e-07, No Iterations 0 smoothSolver: Solving for k, Initial residual = 0.05769598, Final residual = 9.754113e-18, No Iterations 2 ExecutionTime = 83.78 s ClockTime = 115 s Time = 1225 -------------------------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, 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 2-phase flow modeling (long) | DS & HB | Main CFD Forum | 0 | January 8, 2000 15:00 |