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

interFoam and alphaEqn.H

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 25, 2012, 04:50
Default interFoam and alphaEqn.H
  #1
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 9
pablodecastillo is on a distinguished road
Hello,

Why in alphaEqn.H is it not added alpha1 = max(min(alpha1, scalar(1)), scalar(0)); i mean this artificial bounded?

Thanks
pablodecastillo is offline   Reply With Quote

Old   April 25, 2012, 05:55
Default
  #2
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Singapore
Posts: 306
Rep Power: 8
Phicau is on a distinguished road
Hi

MULES solver ensures boundedness, that is why the last 2 inputs are 0 and 1, the bounds of alpha1. This is achieved by the use of limiters.

The final result may show values greater than 1 or lower than 0, but with differences of the order of 10^-5, which is negligible.

Regards
Phicau is offline   Reply With Quote

Old   May 1, 2012, 08:52
Default
  #3
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 430
Rep Power: 16
santiagomarquezd will become famous soon enough
Using an artificial bounding of this kinds leads to a non-phase-preserving algorithm due it 'cuts' with having in account the conservative properties of the alpha equation. In other words you will lose or gain mass.

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   May 7, 2012, 13:57
Default
  #4
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 9
pablodecastillo is on a distinguished road
But looking for steady state loss or gain mass is not going to be relevant or yes?

Thanks
pablodecastillo is offline   Reply With Quote

Old   May 8, 2012, 12:59
Default
  #5
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 430
Rep Power: 16
santiagomarquezd will become famous soon enough
Hhmm, I don't know, but since the solver is unsteady probably it will blow up before reach the steady state due mass preservation problems.

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   March 28, 2016, 10:47
Default
  #6
Senior Member
 
Saideep
Join Date: Apr 2015
Location: INDIA
Posts: 144
Rep Power: 3
Saideep is on a distinguished road
Hi guys,

I guess this is the most relevant thread I found to post the following question and hope anyone of you could help me out.
The alphaEqn.H file has evolved from 2.2 version to 2.3 version. I can understand the 2.2 version of the code but I cant follow most part of the 2.3 version of alphaEqn. Can you help me out?

What does the following code indicate?
fvScalarMatrix alpha1Eqn
(
fv::EulerDdtScheme<scalar>(mesh).fvmDdt(alpha1)
+ fv::gaussConvectionScheme<scalar>
(
mesh,
phi,
upwind<scalar>(mesh, phi)
).fvmDiv(phi, alpha1)
);

The later part is from the older version for advection of alpha with velocity and relative velocity terms.

I tried to use the older version into this but alpha is no more bounded between 0 and 1.

Thanks;
Saideep
Saideep is offline   Reply With Quote

Old   April 16, 2016, 21:00
Default
  #7
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 430
Rep Power: 16
santiagomarquezd will become famous soon enough
Hi Saideep, the new version of alphaEqn is based on operator splitting techniques to obtain a semi-implicit solver. The code you posted corresponds to the first-order bounded predictor of alpha which is lately corrected by a high-order corrector:

http://www.openfoam.org/version2.3.0/multiphase.php

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
http://www.cimec.org.ar
santiagomarquezd 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
alphaEqn.H in twoPhaseEulerFoam cheng1988sjtu OpenFOAM Bugs 15 May 1, 2016 16:12
Small time step in interFoam Andrea_85 OpenFOAM 34 August 22, 2013 22:24
compressibilty effect in alphaEqn.H nimasam OpenFOAM 20 May 25, 2012 08:57
Vof method in interFoam Andrea_85 OpenFOAM 10 February 17, 2011 05:41


All times are GMT -4. The time now is 05:43.