[meganrw]

Hi, I am trying to run a turbulent car simulation where a person coughs for 0.5 seconds with a wind coming in through the window. My friend ran the simulation in ansys and sent me the mesh file. However I keep getting errors - I get the error
" #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:? "
which I know is usually a mesh or boundary condition related error but I just can't figure out what the issue actually is. The mesh should be okay as the simulation runs in ansys so the issue is probably in the boundary. Thanks in advance !

p

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

internalField uniform 0;

boundaryField

{
inlet_wind
{
type zeroGradient;
}
outlet_wind
{
type fixedValue;
value 0;
}
inlet_mouth_front
{
type zeroGradient;
}
inlet_mouth_rear
{
type zeroGradient;
}
wall_left
{
type zeroGradient;
}
wall_right
{
type zeroGradient;
}
wall_bottom
{
type zeroGradient;
}
wall_top
{
type zeroGradient;
}
wall_front
{
type zeroGradient;
}
wall_rear
{
type zeroGradient;
}
wall_window_front_right_1
{
type zeroGradient;
}
wall_window_front_right_2
{
type zeroGradient;
}
wall_window_front_left_1
{
type zeroGradient;
}
wall_window_front_left_2
{
type zeroGradient;
}
wall_window_rear_right_1
{
type zeroGradient;
}
wall_window_rear_right_2
{
type zeroGradient;
}
wall_window_rear_left_1
{
type zeroGradient;
}
wall_window_rear_left_2
{
type zeroGradient;
}
wall_tunnel
{
type zeroGradient;
}
wall_head
{
type zeroGradient;
}
wall_seat_right
{
type zeroGradient;
}
wall_seat_left
{
type zeroGradient;
}
wall_trunk
{
type zeroGradient;
}
wall_seat_rear
{
type zeroGradient;
}
}




U

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

internalField uniform (0 0 0);

boundaryField
{
inlet_wind
{
type fixedValue;
value uniform (15 0 0);
}

outlet_wind
{
type inletOutlet;
inletValue uniform (0 0 0);
}

inlet_mouth_front
{
type uniformFixedValue;
uniformValue table
(
( 0 (0 10 0))
(0.5 (0 10 0 ))
(0.5001 (0 0 0 ))
);
}

inlet_mouth_rear
{
type noSlip;
}

wall_left
{
type noSlip;
}

wall_right
{
type noSlip;
}

wall_bottom
{
type noSlip;
}

wall_top
{
type noSlip;
}

wall_front
{
type noSlip;
}

wall_rear
{
type noSlip;
}

wall_window_front_right_1
{
type noSlip;
}

wall_window_front_right_2
{
type noSlip;
}

wall_window_front_left_1
{
type noSlip;
}

wall_window_front_left_2
{
type noSlip;
}

wall_window_rear_right_1
{
type noSlip;
}

wall_window_rear_right_2
{
type noSlip;
}

wall_window_rear_left_1
{
type noSlip;
}

wall_window_rear_left_2
{
type noSlip;
}

wall_tunnel
{
type noSlip;
}

wall_head
{
type noSlip;
}

wall_seat_right
{
type noSlip;
}

wall_seat_left
{
type noSlip;
}

wall_trunk
{
type noSlip;
}

wall_seat_rear
{
type noSlip;
}
}





nut

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

internalField uniform 0;

boundaryField

{
inlet_wind
{
type calculated;
value uniform 0;
}
outlet_wind
{
type calculated;
value uniform 0;
}
inlet_mouth_front
{
type calculated;
value uniform 0;
}
inlet_mouth_rear
{
type nutkWallFunction;
value uniform 0;
}
wall_left
{
type nutkWallFunction;
value uniform 0;
}
wall_right
{
type nutkWallFunction;
value uniform 0;
}
wall_bottom
{
type nutkWallFunction;
value uniform 0;
}
wall_top
{
type nutkWallFunction;
value uniform 0;
}
wall_front
{
type nutkWallFunction;
value uniform 0;
}
wall_rear
{
type nutkWallFunction;
value uniform 0;
}
wall_window_front_right_1
{
type nutkWallFunction;
value uniform 0;
}
wall_window_front_right_2
{
type nutkWallFunction;
value uniform 0;
}
wall_window_front_left_1
{
type nutkWallFunction;
value uniform 0;
}
wall_window_front_left_2
{
type nutkWallFunction;
value uniform 0;
}
wall_window_rear_right_1
{
type nutkWallFunction;
value uniform 0;
}
wall_window_rear_right_2
{
type nutkWallFunction;
value uniform 0;
}
wall_window_rear_left_1
{
type nutkWallFunction;
value uniform 0;
}
wall_window_rear_left_2
{
type nutkWallFunction;
value uniform 0;
}
wall_tunnel
{
type nutkWallFunction;
value uniform 0;
}
wall_head
{
type nutkWallFunction;
value uniform 0;
}
wall_seat_right
{
type nutkWallFunction;
value uniform 0;
}
wall_seat_left
{
type nutkWallFunction;
value uniform 0;
}
wall_trunk
{
type nutkWallFunction;
value uniform 0;
}
wall_seat_rear
{
type nutkWallFunction;
value uniform 0;
}
}





k

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

internalField uniform 0.375;

boundaryField

{
inlet_wind
{
type fixedValue;
value uniform 0.375;
}
outlet_wind
{
type zeroGradient;
}
inlet_mouth_front
{
type fixedValue;
value uniform 0.375;
}
inlet_mouth_rear
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_left
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_right
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_bottom
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_top
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_front
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_rear
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_front_right_1
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_front_right_2
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_front_left_1
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_front_left_2
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_rear_right_1
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_rear_right_2
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_rear_left_1
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_window_rear_left_2
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_tunnel
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_head
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_seat_right
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_seat_left
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_trunk
{
type kqRWallFunction;
value uniform 0.00375;
}
wall_seat_rear
{
type kqRWallFunction;
value uniform 0.00375;
}
}




omega

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

internalField uniform 440.15;

boundaryField

{
inlet_wind
{
type fixedValue;
value $internalField;
}
outlet_wind
{
type zeroGradient;
}
inlet_mouth_front
{
type fixedValue;
value $internalField;
}
inlet_mouth_rear
{
type omegaWallFunction;
value uniform 22.4;
}
wall_left
{
type omegaWallFunction;
value uniform 22.4;
}
wall_right
{
type omegaWallFunction;
value uniform 22.4;
}
wall_bottom
{
type omegaWallFunction;
value uniform 22.4;
}
wall_top
{
type omegaWallFunction;
value uniform 22.4;
}
wall_front
{
type omegaWallFunction;
value uniform 22.4;
}
wall_rear
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_front_right_1
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_front_right_2
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_front_left_1
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_front_left_2
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_rear_right_1
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_rear_right_2
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_rear_left_1
{
type omegaWallFunction;
value uniform 22.4;
}
wall_window_rear_left_2
{
type omegaWallFunction;
value uniform 22.4;
}
wall_tunnel
{
type omegaWallFunction;
value uniform 22.4;
}
wall_head
{
type omegaWallFunction;
value uniform 22.4;
}
wall_seat_right
{
type omegaWallFunction;
value uniform 22.4;
}
wall_seat_left
{
type omegaWallFunction;
value uniform 22.4;
}
wall_trunk
{
type omegaWallFunction;
value uniform 22.4;
}
wall_seat_rear
{
type omegaWallFunction;
value uniform 22.4;
}
}





epsilon

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

internalField uniform 14.855;

boundaryField

{
inlet_wind
{
type fixedValue;
value uniform 14.855;
}
outlet_wind
{
type zeroGradient;
value uniform 0;
}
inlet_mouth_front
{
type fixedValue;
value uniform 14.855;
}
inlet_mouth_rear
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_left
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_right
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_bottom
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_top
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_front
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_rear
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_front_right_1
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_front_right_2
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_front_left_1
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_front_left_2
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_rear_right_1
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_rear_right_2
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_rear_left_1
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_window_rear_left_2
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_tunnel
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_head
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_seat_right
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_seat_left
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_trunk
{
type epsilonWallFunction;
value uniform 0.00754;
}
wall_seat_rear
{
type epsilonWallFunction;
value uniform 0.00754;
}
}




SYSTEM



controlDict

application simpleFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 10;

deltaT 0.0005;

writeControl timeStep;

writeInterval 100;

purgeWrite 0;

writeFormat binary;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;





fvSchemes

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(U) cellLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwindV grad(U);
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}





fvSolution

solvers
{
p
{
solver GAMG;
smoother GaussSeidel;
tolerance 1e-7;
relTol 0.01;
}

Phi
{
$p;
}

U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}

k
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}

omega
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-8;
relTol 0.1;
nSweeps 1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
consistent yes;
}

potentialFlow
{
nNonOrthogonalCorrectors 10;
}

relaxationFactors
{
equations
{
U 0.9;
k 0.7;
omega 0.7;
}
}

cache
{
grad(U);
}