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

Numerical problem in species Transport in InterFoam at the interface

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

Like Tree1Likes
  • 1 Post By Cyp

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 11, 2018, 06:15
Default Numerical problem in species Transport in InterFoam at the interface
  #1
Member
 
Mehdi Aminyavari
Join Date: Feb 2016
Location: Milan
Posts: 35
Rep Power: 10
Mehdi3031 is on a distinguished road
Hello Forum;
I am trying to introduce a species transport equation in interfoam for the liquid phase (Alpha1 = 1) as bellow:


Code:
Info << "solve ylEqn.H"<<endl;
    surfaceScalarField rho1Phialpha1 = rho1*phi*fvc::interpolate(alpha1);
    fvScalarMatrix ylEqn
    (
        fvm::ddt(alpha1*rho1, yl)
      + fvm::div(rho1Phialpha1, yl, "div(rhoPhiAlpha,yl)")
      - fvm::laplacian(alpha1*rho1*Dl, yl)
    );
This equation has been cited 100s of times and thats how people cancell out the gas domain(Alpha1=0) when solving the species balanse for the liquid phase and just solves it where ever Alpha1 is 1 or between 0 and 1.

But the problem is that I get crazy yl values at the interface where 0<Alpha1<1 as you can see in the attached photo! and if I don't bound the yl between 0 and 1 by force, it will go to infinity.
any Idea how to solve this numerical problem? did anybody have done it before?

Thanks a lot in advance
Attached Images
File Type: jpg Alpha1.jpg (80.0 KB, 55 views)
File Type: jpg yl.jpg (79.2 KB, 56 views)
Mehdi3031 is offline   Reply With Quote

Old   October 11, 2018, 10:20
Default
  #2
Member
 
Mehdi Aminyavari
Join Date: Feb 2016
Location: Milan
Posts: 35
Rep Power: 10
Mehdi3031 is on a distinguished road
why my posts have always minimum views among the list??????!!!!
is there a trick or something to increase the view and reply??
Mehdi3031 is offline   Reply With Quote

Old   October 12, 2018, 03:45
Default
  #3
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
Hi,

Please, have a look at this paper :

https://www.sciencedirect.com/science/article/pii/S0009250918303944


Cheers,
Cyprien
Cyp is offline   Reply With Quote

Old   October 12, 2018, 08:01
Default
  #4
Member
 
Mehdi Aminyavari
Join Date: Feb 2016
Location: Milan
Posts: 35
Rep Power: 10
Mehdi3031 is on a distinguished road
Dear Cyprien,
Thank you so much for your response, I have seen the article, but the point is that you could use the Henry coefficient in this work and make a CST equation because you had a dilute solution which in my case I don't and cant do that. I have to solve the species equation separately for both phases like for phase 1:

Code:
fvScalarMatrix ylEqn
    (
        fvm::ddt(alpha1*rho1, yl)
      + fvm::div(rho1Phialpha1, yl, "div(rhoPhiAlpha,yl)")
      - fvm::laplacian(alpha1*rho1*Dl, yl)
    );
and the same for Phase 2; and then couple them with a source term. What I can't solve is that as soon as Alpha comes around in the species transfer equations, All goes wrong due to numerical problems introduced at the interface:

Any suggestion on that?
Mehdi3031 is offline   Reply With Quote

Old   October 12, 2018, 08:03
Default
  #5
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
You can if you put H=0.
Mehdi3031 likes this.
Cyp is offline   Reply With Quote

Old   October 12, 2018, 08:28
Default
  #6
Member
 
Mehdi Aminyavari
Join Date: Feb 2016
Location: Milan
Posts: 35
Rep Power: 10
Mehdi3031 is on a distinguished road
I may explained bad,
(Using the notion in your paper): Cj,2 = Hj*Cj,1
so if Hj=0, means Cj,2 = 0 which means in phase two I have no Cj!!
In another word phase 2 (gas phase) is pure species i and liquid phase is binary mixture of Species j and i with Cj,1!
My problem is Ammonia-water mixture which means in both phases I have both species and it is non-diluted (Cj,2 = 0.95-0.995) and (Cj,1=0.1-0.6) which is as far as I know, not the case to use Henry coefficient.

Isn't that right?
Mehdi3031 is offline   Reply With Quote

Reply

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
UDF for Sorption in Species Transport through Porous Medium Kosch Fluent UDF and Scheme Programming 6 October 21, 2020 18:46
Species transport model problem argeus FLUENT 8 September 20, 2019 08:23
InterFoam and species transport PKeller OpenFOAM Running, Solving & CFD 7 October 7, 2014 03:40
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 18:44
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 21:09


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