MoeCFD |
July 7, 2017 12:55 |
Problem with SnappyHexMesh Cell Removal
Hello All,
I have to admit that my level in OpenFoam is basic. I am going through the learning process. So, without further ado, I am trying to run simpleFoam on a fibrous material (think about it as a porous media). I created the mesh using SnappyHexMesh and and an stl file that I have. It is 3D problem with one inlet and one outlet, below are the boundary conditions that I am using:
For pressure:
Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
fixedWalls_patch20083
{
type zeroGradient;
}
left
{
type fixedValue;
value uniform 1e-3;
}
right
{
type fixedValue;
value uniform 0;
}
bottom
{
type zeroGradient;
}
top
{
type zeroGradient;
}
frontAndBack
{
type zeroGradient;
}
}
For velocity
Code:
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
fixedWalls_patch20083
{
type fixedValue;
value uniform (0 0 0);
}
left
{
type zeroGradient;
}
right
{
type zeroGradient;
}
bottom
{
type fixedValue;
value uniform (0 0 0);
}
top
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type fixedValue;
value uniform (0 0 0);
}
}
Here is my fvSolution:
Code:
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.1;
}
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0;
}
}
/*
U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-12;
relTol 0;
nSweeps 1;
}
p
{
solver GAMG;
tolerance 1e-9;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration on;
agglomerator faceAreaPair;
nCellsInCoarsestLevel 10;
mergeLevels 1;
}
*/
SIMPLE
{
nNonOrthogonalCorrectors 0;
residualControl
{
p 1e-6;
U 1e-6;
}
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
p 0.15;
U 0.9;
}
// ************************************************************************* //
When I run the command SimpleFoam to start the simulation based on the inputs shown above, I get the following (see below)
which doesn't show me what went wrong with the simulation. Is it something with the way I defined the boundary conditions or simulation settings? I appreciate your help in advance. Also, I am running the simulation in serial for now. However, I created the Mesh by running SnappyHexMesh in parallel then used reconstruct to put it together.
Code:
Build : 4.1
Exec : simpleFoam
Date : Jul 07 2017
Time : 21:57:43
Host :
PID : 592
Case :
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
SIMPLE: convergence criteria
field p tolerance 1e-06
field U tolerance 1e-06
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
No MRF models present
No finite volume options present
Starting time loop
Time = 1
DILUPBiCG: Solving for Ux, Initial residual = 0.219962, Final residual = 2.60365e-08, No Iterations 10
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.89906e-08, No Iterations 11
DILUPBiCG: Solving for Uz, Initial residual = 0.143777, Final residual = 9.97743e-08, No Iterations 9
moe@DESKTOP:/example$
|