CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

rhoSimpleFoam High Pressure Cell Crashes Simulation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 3, 2023, 05:02
Default rhoSimpleFoam High Pressure Cell Crashes Simulation
  #1
New Member
 
Kuzey Can Derman
Join Date: Oct 2022
Posts: 12
Rep Power: 3
NorthCFD is on a distinguished road
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;
    }
}
p:
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;
    }
}
T:
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;
    }
}
Solver properties:
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);
}
fvSchemes:
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;
}
fvOptions:
Code:
limitT
{
    type       limitTemperature;
    min        100;
    max        350;
    selectionMode all;
}
Last entries to terminal:
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 1223 showing pressure:


Photo from step 1100 showing pressure:


Do not hesitate to ask anything. Thanks for any help
NorthCFD is offline   Reply With Quote

Reply

Tags
crash, openfoam, rhosimplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 17:01.