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

Solver for Diffusion of Species

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 16, 2018, 05:15
Default Solver for Diffusion of Species
  #1
New Member
 
manu ebn
Join Date: Aug 2015
Location: Switzerland
Posts: 18
Rep Power: 11
Triggin is on a distinguished road
Hi to everyone

Could you please give me an advice to choose the right solver?

Case:
I have a domain which is filled with oxygen (< 200ppm) and nitrogen (>99.9%). Furthermore there is an outlet and an inlet.
Through the inlet pure nitrogen flows into the domain, at the outlet nitrogen flows out of the domain (with Reynoldsnumber lower than 50). But, and here comes the tricky part, oxygen flows from the environment into the domain (probably by diffusion) - at the end i would like to have a steady state solution (laminar).

I could simulate this with a scalarTransportFoam, but there wouldn't be any interaction between the different species and in a steady-state solution everything would be just nitrogen.

So which solver would you advice?

Thanks for every single answer!
Triggin
Triggin is offline   Reply With Quote

Old   January 19, 2018, 15:28
Default
  #2
New Member
 
manu ebn
Join Date: Aug 2015
Location: Switzerland
Posts: 18
Rep Power: 11
Triggin is on a distinguished road
If someone could help me with this one, i would be very thankful.

Like i sad, scalarTransport could be an option, but i would like to know if it's possible to solve the interaction of a binary mixture? I think it could be tricky, because to simulate the interaction of particels in flow is already demanding.

This information would already be helpful.

Best wishes
Triggin
Triggin is offline   Reply With Quote

Old   January 22, 2018, 04:42
Default
  #3
New Member
 
Sebastian
Join Date: Feb 2017
Posts: 22
Rep Power: 9
sepp.zell is on a distinguished road
Hi Triggin,

I am a little bit confused. You are writing that you want to use scalarTransportFoam, however, later you write that you have interaction of particles. Usually scalarTransportFoam works with the concentration of the quantity of interest and not with single particles. Sticking to concentration, you might have to look for a multiphase flow solver. I don't know if interFoam could be a suitable solver for your case.
Concerning the inflow oxygen: If you can not model it as a classical inflow, you might add a source term, possibly using fvOptions.

Regards

Sebastian
sepp.zell is offline   Reply With Quote

Old   January 22, 2018, 08:58
Default
  #4
New Member
 
manu ebn
Join Date: Aug 2015
Location: Switzerland
Posts: 18
Rep Power: 11
Triggin is on a distinguished road
Hi Sebastian

Thank you for your comments.

Quote:
You are writing that you want to use scalarTransportFoam, however, later you write that you have interaction of particles. Usually scalarTransportFoam works with the concentration of the quantity of interest and not with single particles.
You are completely right. I a modified scalarTransport solver does it.

I am not very familiar with multi-gas simulations (therefore i also thought about reactingFoam - without reactions). I thought the different gases (nitrogen and oxygen) might interact and turbulences even at low reynoldsnumbers can occur.

But assuming that i have only a very low oxygen concentration and the flow is mainly influenced by nitrogen, i can use a modified scalarTransportFoam for steady state cases.

Quote:
Concerning the inflow oxygen: If you can not model it as a classical inflow, you might add a source term, possibly using fvOptions.
I used the following conditions, U-file:
Code:
    outlet
    {
      type              pressureInletOutletVelocity;
      value             $internalField;
    }
C-file for concentrations:
Code:
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0.23;
        value           $internalField;
     }
Therefore: If flow goes out of the domain, i will have zeroGradient, if air enters the domain it will enter with a concentration of 0.77 nitrogen (an consequently 0.23 oxygen).

Thank you for your comments!
Triggin is offline   Reply With Quote

Old   December 28, 2022, 03:45
Default
  #5
New Member
 
孙睿
Join Date: Dec 2022
Posts: 2
Rep Power: 0
R.Sun is on a distinguished road
Quote:
Originally Posted by Triggin View Post
Hi Sebastian

Thank you for your comments.

You are completely right. I a modified scalarTransport solver does it.

I am not very familiar with multi-gas simulations (therefore i also thought about reactingFoam - without reactions). I thought the different gases (nitrogen and oxygen) might interact and turbulences even at low reynoldsnumbers can occur.

But assuming that i have only a very low oxygen concentration and the flow is mainly influenced by nitrogen, i can use a modified scalarTransportFoam for steady state cases.



I used the following conditions, U-file:
Code:
    outlet
    {
      type              pressureInletOutletVelocity;
      value             $internalField;
    }
C-file for concentrations:
Code:
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0.23;
        value           $internalField;
     }
Therefore: If flow goes out of the domain, i will have zeroGradient, if air enters the domain it will enter with a concentration of 0.77 nitrogen (an consequently 0.23 oxygen).

Thank you for your comments!
Hi, I really want to know that have you managed to get a steady state flow field using scalarTransportFoam? I tried to set the ddtSchemes to steadyState and change the other settings accordingly, but the results were just weird...
R.Sun is offline   Reply With Quote

Reply

Tags
diffusion, laminar, nitrogen, solver, steady-state

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
Hybrid discretisation - blend factor gcoopermax CFX 5 September 23, 2016 08:05
Creating New Solver: For particle-laden compressible jets sankarv OpenFOAM Running, Solving & CFD 17 December 3, 2014 19:41
Working directory via command line Luiz CFX 4 March 6, 2011 20:02
Creating New Solver: For particle-laden compressible jets sankarv OpenFOAM 0 April 4, 2010 18:06
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


All times are GMT -4. The time now is 20:00.