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

compressibleInterFoam - problem with boundary conditions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 25, 2014, 11:38
Default compressibleInterFoam - problem with boundary conditions
  #1
Member
 
Join Date: Feb 2013
Posts: 30
Rep Power: 5
Natalie2210 is on a distinguished road
Hi everybody,

I want to simulate pressure driven flow through a pipe. Up to now I have used the solver interFoam successfully, now, however, I want to test the effect of compressibility. Thus, I'd like to switch to the solver compressibleInterFoam.

I have found the tutorial depthCharge in the tutorials section and found out how to set psi and rho0, at least in the transportProperties-file. However, when I try to run the case I run into troubles. I somehow suspect the boundary conditions. I'm using the boundary conditions working well with interFoam (creating a pressure drop which makes the fluid move through the pipe):

alpha1

HTML Code:
*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      alpha;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 0 0 0 0 0 0];
internalField   uniform 0;
boundaryField
{
    wall
    {
        type            zeroGradient;
    }
    wall2
    {
        type            zeroGradient;
    }
    outflow 
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
   inflow
   {
      type    fixedValue;
      value   1;  
   }
}
// ************************************************************************* //
p_rgh:

HTML Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [1 -1 -2 0 0 0 0];
internalField   uniform 0;
boundaryField
{
    wall
    {
        type           zeroGradient;
    }
    wall2
    {
        type           zeroGradient;
    }
    inflow
    {
        type            fixedValue;
        value           3.5e4;
    }
    outflow 
    {
        type            fixedValue;
        value           $internalField;
    }
}
// ************************************************************************* //
and U:

HTML Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 1 -1 0 0 0 0];
internalField   uniform (0 0 0);
boundaryField
{
    wall
    {
        type            slip;
    }
    wall2
    {
        type            slip;
    }
    inflow 
    {
        type            pressureInletUniformVelocity;
        phi             phi; 
        rho             rho; 
        value           uniform (0 0 0);
    }
    outflow 
    {
        type            pressureInletOutletVelocity;
        phi             phi; 
        value           uniform (0 0 0);
    }
}
// ************************************************************************* //
Do I need to specify something else? Psi for example?

Thanks a lot,
Natalie
Natalie2210 is offline   Reply With Quote

Old   March 26, 2014, 05:19
Default
  #2
Member
 
Join Date: Feb 2013
Posts: 30
Rep Power: 5
Natalie2210 is on a distinguished road
Ok, I think I got it - I just used a time increment way to big in order to get a stable solution.
Natalie2210 is offline   Reply With Quote

Old   April 21, 2016, 12:14
Default
  #3
New Member
 
Paul
Join Date: May 2012
Posts: 23
Rep Power: 6
pmdelgado2 is on a distinguished road
Could you share all of your case files?
pmdelgado2 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
Problem with Periodic Boundary Conditions Help!!! otsigun FLUENT 0 July 11, 2013 03:20
Boundary Conditions problem o_mars_2010 Main CFD Forum 2 July 8, 2013 02:10
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
Problem with Boundary conditions Mahiboobswamulu Main CFD Forum 10 August 26, 2003 13:24


All times are GMT -4. The time now is 15:35.