CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Problems when change the boundary conditions in rhoSimplicFoam tutorial example

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 16, 2014, 09:52
Default Problems when change the boundary conditions in rhoSimplicFoam tutorial example
  #1
New Member
 
Join Date: Dec 2014
Posts: 10
Rep Power: 3
SylvainWang is on a distinguished road
Hello, Foamers;

I just begin to use OpenFOAM to do my thesis. Since my problem is a compressible case, so I'd like to choose the rhoSimplecFoam solver to solve it. In the tutorial example offered in the rhoSimplecFoam solver, it gives a mass flowrate inlet condition and pressure outlet condition. However, in my case, it should be a total pressure inlet condition and pressure outlet condition. Problem comes as I changed the boundary conditions of U and P after 3 iterations. The error messages are as following:
Quote:
#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::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception (core dumped)
and the new boundary conditions are
P
Code:
dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 100000;

boundaryField
{
    Default_Boundary_Region
    {
        type            zeroGradient;
    }
    inlet
    {
        type            totalPressure;
        U               U;
        phi             phi;
        rho             none;
        psi             none;
        gamma           1.4;
        p0              uniform 200000;
    }
    outlet
    {
        type            inletOutlet;
        phi             phi;
        inletValue      uniform 100000;
        value           uniform 100000;
    }
}
U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    Default_Boundary_Region
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    inlet
    {
        type            pressureInletVelocity;
        phi             phi;
        rho             rho;
        value           uniform (0 0 0);
    }
    outlet
    {
        type            pressureInletOutletVelocity;
        phi             phi;
        value           uniform (0 0 0);
    }
}
The version I used is 2.3.0, from the error message, it seems that the errors comes from the thermal properties. Nevertheless, we should not encounter any errors when just change the velocity inlet boundary condition to a total pressure boundary condition. Is there any problems in the modified boundary conditions?

Endless appreciation if any foamers can help. Sylvain
SylvainWang is offline   Reply With Quote

Old   December 17, 2014, 16:43
Default
  #2
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 328
Rep Power: 12
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi Sylvain,

I think you should change your pressure boundary condition for the outlet to either a fixed static pressure or to outletInlet instead of inletOutlet. Look up their differences in the sourcecode!

Some hints:

Code:
$> cd $FOAM_SRC
$> find . -iname *inletoutlet*
If you have a supersonic outlet, zeroGradient or waveTransmissive may be better options.

You may also need to ramp your total pressure at the inlet over several iterations. I think some examples may be found on the forum.

Regards,
Tom
tomf is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Implementation of boundary conditions for FVM Tom Main CFD Forum 7 August 26, 2014 05:58
Boundary conditions for hassan hemida interFoam tutorial vishal_s OpenFOAM Pre-Processing 0 August 21, 2013 00:21
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
[ICEM] Problems with coedge curves and surfaces tommymoose ANSYS Meshing & Geometry 0 August 5, 2011 16:02


All times are GMT -4. The time now is 23:57.