[mohammad morovati]

Hi,


I am trying to simulate a fan inside a canal by rhoPimpleFoam solver.



Fan is rotating by w=104, and I face to negative initial temperature in 3rth part of first iteration.


It would be great appreciated if someone could guide me.

all files are attached

[Fouch]

Hi,

For compressible cases you have to limit T and U in fvOption (constant directory) :

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

/*viscousDissipation
{
    type            viscousDissipation;
    enabled         true;
}*/

limitT
{
    type       limitTemperature;
    min        230;
    max        400;
    selectionMode all;
}
limitU
{
    type            limitVelocity;
    selectionMode   all;
    max            400;
}
// ************************************************************************* //

This is an example above but you have to adapt the values to your case.

Best regards,