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

Two phases flow with implementation of chemical reactions

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

Like Tree3Likes
  • 1 Post By juho
  • 1 Post By RamCmd
  • 1 Post By A.jafarinia

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 19, 2016, 03:46
Post Two phases flow with implementation of chemical reactions
  #1
New Member
 
M.H. Radanielina
Join Date: Oct 2016
Location: Madagascar
Posts: 4
Rep Power: 5
RamCmd is on a distinguished road
Hi everyone,

I am a new member. I am very worry because, I work on my topic til several months but I have no idea how to find solution on this.

I am currently working on gasification process and using OpenFOAM to solve it. The mathematical model I would like to implement consists of two phases (gas and solid) and chemical reactions.

I tried to use reatingFoam solver but it only considers single phase, whereas multiple/twoPhaseEuleurFoam or reactingTwo/multiplePhaseEulerFoam solvers allow to describe two phases but no chemical reaction inside it. I thought in first on parallelisation calculation with two solvers but I found it is difficult or impossible to implement above all I have not a lot experience about OpenFoam manipulation. I can mesh, define my boundary conditions, but I could not directly enter inside .C or .H file to modify it.

Could someone help me how I can find solution on it? and if there are other solutions.

Thanks in advance for all your replies.

Best Regards,

RamCmd
RamCmd is offline   Reply With Quote

Old   October 19, 2016, 04:14
Default
  #2
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 89
Rep Power: 12
juho is on a distinguished road
I think that the reactingEulerFoam phaseModel does support reactions:

https://github.com/OpenFOAM/OpenFOAM...ntPhaseModel.C
Code:
    return
    (
        fvm::ddt(alpha, rho, Yi)
      + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)")
      - fvm::Sp(this->continuityError(), Yi)

      - fvm::laplacian
        (
            fvc::interpolate(alpha)
           *fvc::interpolate(this->turbulence().nut()*rho/Sc_),
            Yi
        )
     ==
        this->R(Yi)

      + fvc::ddt(residualAlpha_*rho, Yi)
      - fvm::ddt(residualAlpha_*rho, Yi)
    );
https://github.com/OpenFOAM/OpenFOAM...ngPhaseModel.C

Code:
Foam::ReactingPhaseModel<BasePhaseModel, ReactionType>::R
(
    volScalarField& Yi
) const
{
    return reaction_->R(Yi);
}
A tutorial case with reactions:

https://github.com/OpenFOAM/OpenFOAM...ratingReacting

Are there some components still missing?
RamCmd likes this.
juho is offline   Reply With Quote

Old   October 19, 2016, 05:11
Default
  #3
New Member
 
M.H. Radanielina
Join Date: Oct 2016
Location: Madagascar
Posts: 4
Rep Power: 5
RamCmd is on a distinguished road
Quote:
Originally Posted by juho View Post
I think that the reactingEulerFoam phaseModel does support reactions:

https://github.com/OpenFOAM/OpenFOAM...ntPhaseModel.C
Code:
    return
    (
        fvm::ddt(alpha, rho, Yi)
      + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)")
      - fvm::Sp(this->continuityError(), Yi)

      - fvm::laplacian
        (
            fvc::interpolate(alpha)
           *fvc::interpolate(this->turbulence().nut()*rho/Sc_),
            Yi
        )
     ==
        this->R(Yi)

      + fvc::ddt(residualAlpha_*rho, Yi)
      - fvm::ddt(residualAlpha_*rho, Yi)
    );
https://github.com/OpenFOAM/OpenFOAM...ngPhaseModel.C

Code:
Foam::ReactingPhaseModel<BasePhaseModel, ReactionType>::R
(
    volScalarField& Yi
) const
{
    return reaction_->R(Yi);
}
A tutorial case with reactions:

https://github.com/OpenFOAM/OpenFOAM...ratingReacting

Are there some components still missing?
Juho, thanks a lot for your reply, it's very kind of you

I think there are all components I need to solve my case study. I 'll try with it.

I'll keep you posted.
RamCmd is offline   Reply With Quote

Old   October 21, 2016, 04:37
Default
  #4
New Member
 
M.H. Radanielina
Join Date: Oct 2016
Location: Madagascar
Posts: 4
Rep Power: 5
RamCmd is on a distinguished road
Hi Juho,

I saw the tuto case bubbleColumnEvaporatingReacting you told me, and have some questions about this, I post you only two for the moment because it is the base of my study:

1. in constant/reactions.gas file, there are species
species
(
CO
CO2
H2
H2O
AIR
);
So if I want to introduce other species like in my case, CHAR, it should be a SOLID particle interpreted as a CONTINUUM phase, where can I define its properties (thermo and physical)? is that in the thermo.gas file or in another one because it is not a gas?

2. in constant/phaseProperties file, could you explain me how the type of phase is define? and how I can define the type of my phase properties like in this :
gas
{
type reactingPhaseModel;//
...
}

liquid // change into dispersed solid phase
{
type purePhaseModel; //
...
}
Is it already in OpenFoan, and how I can use it?

Thanks for time and consideration.

Regards,

RamCmd
Kummi likes this.
RamCmd is offline   Reply With Quote

Old   June 19, 2018, 05:11
Default
  #5
New Member
 
Alireza Jafarinia
Join Date: Jun 2018
Location: Austria
Posts: 6
Rep Power: 3
A.jafarinia is on a distinguished road
Quote:
Originally Posted by RamCmd View Post
Hi Juho,

I saw the tuto case bubbleColumnEvaporatingReacting you told me, and have some questions about this, I post you only two for the moment because it is the base of my study:

1. in constant/reactions.gas file, there are species
species
(
CO
CO2
H2
H2O
AIR
);
So if I want to introduce other species like in my case, CHAR, it should be a SOLID particle interpreted as a CONTINUUM phase, where can I define its properties (thermo and physical)? is that in the thermo.gas file or in another one because it is not a gas?

2. in constant/phaseProperties file, could you explain me how the type of phase is define? and how I can define the type of my phase properties like in this :
gas
{
type reactingPhaseModel;//
...
}

liquid // change into dispersed solid phase
{
type purePhaseModel; //
...
}
Is it already in OpenFoan, and how I can use it?

Thanks for time and consideration.

Regards,

RamCmd
Hi. did you manage to do it? I have a similar problem.
I am going to solve convection-diffusion-reaction equations in multiphase flow using OpenFOAM 5. As far as I know it is possible to use reactingMultiphaseEulerFoam solver, but the point is that there is no solved tutorial similar to my case. I would like to know your recommendations about how to start. Which tutorial do you suggest to modify?
There is another solver reactingTwophaseEulerFoam which solves convection-diffusion-reaction in two phase. Currently I think it is possible to start with this two phase flow problem or a tutorial in multiphase flow without solving reaction like multiphaseEulerFoam. I am really confused how to start and I appreciate any kind of help. My equation are related to biochemical reactions in multiphase flow with several species like thrombus formation in blood flow. The problem is the formation of thrombus in blood. It means that at first we have only fluid which is blood; then after biochemical reactions a solid phase might form near the wall which is the thrombus.
Kummi likes this.
A.jafarinia is offline   Reply With Quote

Old   June 18, 2019, 10:07
Default
  #6
New Member
 
Sepehr
Join Date: May 2019
Posts: 1
Rep Power: 0
sepehrmk is on a distinguished road
Quote:
Originally Posted by A.jafarinia View Post
Hi. did you manage to do it? I have a similar problem.
I am going to solve convection-diffusion-reaction equations in multiphase flow using OpenFOAM 5. As far as I know it is possible to use reactingMultiphaseEulerFoam solver, but the point is that there is no solved tutorial similar to my case. I would like to know your recommendations about how to start. Which tutorial do you suggest to modify?
There is another solver reactingTwophaseEulerFoam which solves convection-diffusion-reaction in two phase. Currently I think it is possible to start with this two phase flow problem or a tutorial in multiphase flow without solving reaction like multiphaseEulerFoam. I am really confused how to start and I appreciate any kind of help. My equation are related to biochemical reactions in multiphase flow with several species like thrombus formation in blood flow. The problem is the formation of thrombus in blood. It means that at first we have only fluid which is blood; then after biochemical reactions a solid phase might form near the wall which is the thrombus.
Hi Alireza
Is there any update/solution for your problem?
sepehrmk is offline   Reply With Quote

Old   June 19, 2019, 06:07
Default
  #7
New Member
 
Alireza Jafarinia
Join Date: Jun 2018
Location: Austria
Posts: 6
Rep Power: 3
A.jafarinia is on a distinguished road
I haven't tried tried the multiphase solvers yet, currently I am solving only passive scalar transport equations for the change in species concentrations.
A.jafarinia is offline   Reply With Quote

Reply

Tags
gasification cfd, parallel calculation, reactingfoam, twophaseeulerfoam

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
About Some Concepts:Laminar flow, turbulent flow, steady flow and time-dependent flow Jing Main CFD Forum 8 October 5, 2018 18:02
Review: Reversed flow CRT FLUENT 1 May 7, 2018 06:36
Can 'shock waves' occur in viscous fluid flows? diaw Main CFD Forum 104 February 16, 2006 06:44
How to specify chemical reactions in certain areas Esti Siemens 0 September 24, 2004 07:03
Finite rate chemical reactions Marc Segovia Siemens 3 November 28, 2001 05:21


All times are GMT -4. The time now is 08:48.