[Juzzvy]

Good day. I am currently running a simulation in openFoam using the simple method, with the kEpsilon turbulence model. The run was going fine, but after I tried to experimentally change the wind direction, this error appears

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve() at ??:?
#9 ? at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? at ??:?
Floating point exception (core dumped)

I don't understand yet what could have gone wrong. the previous run it was ok, though I changed the inlet and outlet so as to suit with the wind direction I am aiming at. Thank you for the time and your suggestions are appreciated

[Joshua14]

Juzzvy,

Your "Osteam" error most likely is related to what you changed your boundary conditions to. Double check to see that they are correct or that you don't need to add additional changes based on the changes you made.

Feel free to post your old and new boundary conditions so that people can visually look at what you have.

Joshua

[Juzzvy]

Joshua,

Basically I only changed the wind direction, I attached the boundary condition for U which worked (2.83 -2.83 0) and the other that didn't (2.83 2.83 0). I wanted to use the latter, since it will exhibit the northwest direction which I was aiming for. I'll also put here the other boundary conditions/ Sorry as of now I can't attach it due to some reasons

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Uinlet (2.83 2.83 0);

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform $Uinlet;
}

outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value $internalField;
}

rooseveltGroup
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type slip;
}
lowerWall
{
type slip;
}
upperWall
{
type slip;
}
#includeEtc "caseDicts/setConstraintTypes"
}

// ************************************************** *********************** //



/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}

outlet
{
type totalPressure;
p0 uniform 0;
gamma 1.4;
value uniform 0;
}

rooseveltGroup
{
type zeroGradient;
}
frontAndBack
{
type slip;
}
lowerWall
{
type slip;
}
upperWall
{
type slip;
}
#includeEtc "caseDicts/setConstraintTypes"
}

// ************************************************** *********************** //

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}

outlet
{
type calculated;
value uniform 0;
}

rooseveltGroup
{
type nutkWallFunction;
value uniform 0;
}
frontAndBack
{
type slip;
}
lowerWall
{
type slip;
}
upperWall
{
type slip;
}
#includeEtc "caseDicts/setConstraintTypes"
}

// ************************************************** *********************** //


/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

kInlet 1.5; // approx k = 1.5*(I*U)^2 ; I = 0.1

dimensions [0 2 -2 0 0 0 0];

internalField uniform $kInlet;

boundaryField
{
inlet
{
type fixedValue;
value uniform $kInlet;
}

outlet
{
type inletOutlet;
inletValue uniform $kInlet;
value uniform $kInlet;
}

rooseveltGroup
{
type kqRWallFunction;
value uniform $kInlet;
}
frontAndBack
{
type slip;
}
lowerWall
{
type slip;
}
upperWall
{
type slip;
}
#includeEtc "caseDicts/setConstraintTypes"
}

// ************************************************** *********************** //



/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

epsilonInlet 0.03; // Cmu^0.75 * k^1.5 / L ; L =10

dimensions [0 2 -3 0 0 0 0];

internalField uniform $epsilonInlet;

boundaryField
{
inlet
{
type fixedValue;
value uniform $epsilonInlet;
}

outlet
{
type inletOutlet;
inletValue uniform $epsilonInlet;
value uniform $epsilonInlet;
}

rooseveltGroup
{
type epsilonWallFunction;
value uniform $epsilonInlet;
}
frontAndBack
{
type slip;
}
lowerWall
{
type slip;
}
upperWall
{
type slip;
}
#includeEtc "caseDicts/setConstraintTypes"
}

// ************************************************** *********************** //


/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices
(
(-200 -117 0)
(200 -117 0)
(200 177 0)
(-200 177 0)
(-200 -117 200)
(200 -117 200)
(200 177 200)
(-200 177 200)
);

blocks
(
hex (0 1 2 3 4 5 6 7) (14 77 14) simpleGrading (1 1 1)
);

edges
(
);

boundary
(

inlet
{
type patch;
faces
(
(0 4 7 3)
(1 5 4 0)
);
}
outlet
{
type patch;
faces
(
(2 6 5 1)
(3 7 6 2)
);
}
lowerWall
{
type slip;
faces
(
(0 3 2 1)
);
}
upperWall
{
type slip;
faces
(
(4 5 6 7)
);
}
);

// ************************************************** *********************** //