CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

rhoPimpleFoam hardship

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 7, 2016, 02:41
Default rhoPimpleFoam hardship
  #1
New Member
 
Pedro
Join Date: Feb 2013
Location: Beijing (CHINA)
Posts: 3
Rep Power: 13
petrus is on a distinguished road
Send a message via Skype™ to petrus
https://www.dropbox.com/home?preview...%3A24%3A15.png

https://www.dropbox.com/home?preview...%3A29%3A07.png


Hi.

My first post, I hope give all necessary information to get a pleasant, good and quickly answer.
Thanks in advance.

So, trying simulate a jet injection inside a combustion chamber (without combustion of course). Nitrogen gas.
The mesh is axisymmetric, he nozzle of injector and the combustion chamber.
Split in 3 blocks.
blocks
(
hex ( 0 2 2 0 1 3 9 8) (20 1 8) simpleGrading (1.00 1.00 0.3) // block number0
hex ( 2 5 5 2 3 6 11 9) (100 1 8) simpleGrading (10.00 1.00 0.3) // block number1
hex ( 3 6 11 9 4 7 12 10) (100 1 25) simpleGrading (10.00 1.00 35.00) // block number2
);

It run correctly, but the animation is not as expected, better said it is very strange.

p;
}

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

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

internalField uniform 2.14e6;

boundaryField
{

inlet
{
type zeroGradient;
}

outlet
{
type fixedValue;
value uniform 2.14e6;
}

wall
{
type zeroGradient;
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}

}
---------------------------------------------------------------------------------------------
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 298;

boundaryField
{

inlet
{
type fixedValue;
value uniform 298;
}

outlet
{
type zeroGradient;
}

wall
{
type zeroGradient;
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * /

object U;
}

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

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

internalField uniform (0 0 0);

boundaryField
{

inlet
{
type fixedValue;
value uniform (400 0 0);
}

outlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}


wall
{
type fixedValue;
value uniform (0 0 0);
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

object sig;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{

inlet
{
type fixedValue;
value uniform 0;
}

outlet
{
type zeroGradient;
}

wall
{
type fixedValue;
value uniform 0;
}

top
{
type zeroGradient;
}

frontface
{
type wedge;
}

backface
{
type wedge;
}

axis
{
type empty;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
top
{
type calculated;
value uniform 0;
}
wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 14903;

boundaryField
{
inlet
{
type fixedValue;
value uniform 14903;
}
outlet
{
type zeroGradient;
}
top
{
type zeroGradient;
}
wall
{
type compressible::epsilonWallFunction;
value uniform 14903;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 14903;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 37.5;

boundaryField
{
inlet
{
type fixedValue;
value uniform 37.5;
}
outlet
{
type zeroGradient;
}
top
{
type zeroGradient;
}
wall
{
type compressible::kqRWallFunction;
value uniform 37.5;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object mut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
top
{
type calculated;
value uniform 0;
}
wall
{
type mutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

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;
}
top
{
type calculated;
value uniform 0;
}
wall
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 1118000;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1118000;
}
outlet
{
type zeroGradient;
}
top
{
type zeroGradient;
}
wall
{
type omegaWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
beta1 0.075;
value uniform 1118000;
}
frontface
{
type wedge;
}
backface
{
type wedge;
}
axis
{
type empty;
}
}
// ************************************************** *********************** //

blockMesh is ok, checkMesh also is ok.

I think the problem is in the boundary conditions, but I dunno where.



Thanking again any suggestion.
petrus is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
rhoPimpleFoam Boundary Condition Problem dancfd OpenFOAM Pre-Processing 18 September 16, 2021 07:43
rhoPimpleFoam solver adrieno OpenFOAM Running, Solving & CFD 11 April 6, 2016 11:01
Pressure stair-step behaviour using rhopimplefoam joegi.geo OpenFOAM Running, Solving & CFD 3 December 12, 2014 12:10
rhoPimpleFoam floating point error dancfd OpenFOAM Running, Solving & CFD 6 January 5, 2014 20:57
Questions about buoyantPimpleFoam and rhoPimpleFoam Mojtaba.a OpenFOAM Running, Solving & CFD 6 August 1, 2012 04:50


All times are GMT -4. The time now is 11:12.