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

Multispecies jet flow in a garage, HELP !!

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 25, 2010, 04:36
Question Multispecies jet flow in a garage, HELP !!
  #1
New Member
 
Join Date: Apr 2010
Posts: 5
Rep Power: 7
Pierre Mathieu is on a distinguished road
Hi everyone,

I want to simulate the injection of Helium into air in a confined space, and I am quite new with OpenFoam.

I used the reactingFoam solver adding the gravity effect (rho*g) to the velocity equation (UEqn.H) but my results are bad : my species are mixing like if there was no buoyancy effects between them...

I want to figure the density effects between the species during the injection.

If someone have some advices to give, I would be grateful.
Thank you !!
Pierre Mathieu is offline   Reply With Quote

Old   May 27, 2010, 11:04
Default
  #2
Senior Member
 
Eelco van Vliet
Join Date: Mar 2009
Location: The Netherlands
Posts: 122
Rep Power: 9
eelcovv is on a distinguished road
Hi Pierre,

Perhaps you can use something simular as buoyantBoussinesqPisoFoam, which means you only solve one phase (air) and solve helium as a passive scalar A mixing into the room with A=1 at inlet and A=0 at t=0 in the room. Solve the transport of A and add a momentum source to your U equation based on the fraction of helium in air (just a the boussinesq equation). Something like below, but I haven't checked it, so probably needs some adjustments. Goodluck

To solve transport of A (helium fraction)
Code:
{
    volScalarField kappaEff
    (
        "kappaEff",
        turbulence->nu()/Pr + turbulence->nut()/Prt
    );

    fvScalarMatrix AEqn
    (
        fvm::ddt(A)
      + fvm::div(phi, A)
      - fvm::laplacian(kappaEff, A)
    );

    AEqn.relax();

    AEqn.solve();

    rhok = (rho_air-rho_hel)/rho_air*A;
and to solve U

Code:
    fvVectorMatrix UEqn
    (
        fvm::ddt(U)
      + fvm::div(phi, U)
      + turbulence->divDevReff(U)
    );

    UEqn.relax();

    if (momentumPredictor)
    {
        solve
        (
            UEqn
         ==
            fvc::reconstruct
            (
                (
                    fvc::interpolate(rhok)*(g & mesh.Sf())
                  - fvc::snGrad(p)*mesh.magSf()
                )
            )
        );
    }
eelcovv is offline   Reply With Quote

Old   May 27, 2010, 14:08
Default
  #3
Senior Member
 
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 271
Rep Power: 9
phsieh2005 is on a distinguished road
Hi,

Have you looked into twoLiquidMixingFoam?

Pei
phsieh2005 is offline   Reply With Quote

Old   May 28, 2010, 10:33
Default
  #4
New Member
 
Join Date: Apr 2010
Posts: 5
Rep Power: 7
Pierre Mathieu is on a distinguished road
Hi,

I am actually trying with twoLiquidMixingFoam, I am waiting for the results.

Thanks for the advices Eelco, I wil try this way too!
Pierre
Pierre Mathieu 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
how to simulate jet flow? Benny FLUENT 18 May 12, 2014 10:53
Approximate Numerical modelling of Jet Flow jphastings Main CFD Forum 0 March 22, 2010 10:10
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Slot jet in cross flow. Hamdy FLUENT 0 March 5, 2006 05:04
2D Jet Flow Soon CFX 2 August 16, 2005 11:46


All times are GMT -4. The time now is 13:09.