[omid20110]

Hello everybody
I want to simulate sandia Spray H heptane combustion via sprayFoam
Please look at this: http://www.ca.sandia.gov/ecn/cvdata/...h/frameset.php
First I simulated the nonreactive case with O2 set to zero. And after some difficulties I concluded that I should run my case in 3D form (because 2D wedge gave wrong results for LPL) in order to obtain right liquid penetration length (LPL) and vapor penetration length (VPL) so I used a 30 degree sector of a cylinder (the chamber in spray H is a cylinder 108*108 mm) after that I wanted to simulate the reactive case with O2 set to 21% but my case diverges and gives unusual temperatures sometimes it gets negative but most often it rises up to more than 6000 celsius degrees, that is out of Janaf thermo range.
I hope I described everything clearly.
Now can anyone help me to solve these problems and simulate spray H combustion case?

[arvindpj]

Try using the temperature limiters and it has helped me in the past.

/*--------------------------------*- 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;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

temperatureLimit
{
	type limitTemperature;


	active true;
	limitTemperatureCoeffs
	{
		selectionMode all;
		max 5000; //Tmax
		min 250.0;
	}
}/*--------------------------------*- 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;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

temperatureLimit
{
	type limitTemperature;


	active true;
	limitTemperatureCoeffs
	{
		selectionMode all;
		max 5000; //Tmax
		min 250.0;
	}
}

Cheers:
-J

[omid20110]

Originally Posted by arvindpj

Try using the temperature limiters and it has helped me in the past.

/*--------------------------------*- 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;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

temperatureLimit
{
    type limitTemperature;


    active true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        max 5000; //Tmax
        min 250.0;
    }
}/*--------------------------------*- 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;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

temperatureLimit
{
    type limitTemperature;


    active true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        max 5000; //Tmax
        min 250.0;
    }
}

Cheers:
-J

Thanks
How should I do that?
Should I set the min and max correspondence to that of thermo.dat file?
My openFoam version is 4.1 does it have difference?

[omid20110]

Originally Posted by arvindpj

Try using the temperature limiters and it has helped me in the past.

/*--------------------------------*- 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;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

temperatureLimit
{
    type limitTemperature;


    active true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        max 5000; //Tmax
        min 250.0;
    }
}/*--------------------------------*- 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;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

temperatureLimit
{
    type limitTemperature;


    active true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        max 5000; //Tmax
        min 250.0;
    }
}

Cheers:
-J

I tried temperature limiter but it gives error!

[arvindpj]

Could you provide additional details?
-J