|
[Sponsors] |
![]() |
![]() |
#1 |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14 ![]() |
Hi Foamers,
I want to modify the interFoam solver where the alpha Eqn is like: Eqn 1: ![]() The fluid Cell contains fluid phase 1 and phase 2, along with the solid phase, here \epsilon_f stands for the volume fraction occupied by fluid (V_phase1 + V_phase2). \epsilon_f filed is also defined in createFields.H Code:
volScalarField voidfractionFluid ( IOobject ( "voidfractionFluid", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); Eqn 2: ![]() the implementation code of alpha equation: Code:
for (int aCorr=0; aCorr<nAlphaCorr; aCorr++) { surfaceScalarField phiAlpha ( fvc::flux ( phi, alpha1, alphaScheme ) + fvc::flux ( -fvc::flux(-phir, scalar(1) - alpha1, alpharScheme), alpha1, alpharScheme ) ); MULES::explicitSolve(alpha1, phi, phiAlpha, 1, 0); rhoPhi = phiAlpha*(rho1 - rho2) + phi*rho2; Cheers, David Long Last edited by keepfit; October 2, 2013 at 12:21. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 452
Rep Power: 24 ![]() |
Hi, I guess it could be done using the source term included in MULES.
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 |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Mahdi
Join Date: Jul 2012
Posts: 53
Rep Power: 14 ![]() |
Hi Santiago
I followed your PhD thesis where you explain how MULE works. Concerning this problem how the source terms should look like? because MULE needs Su and Sp source terms to solve it implicitly or explicitly, but the question is how to derive Su and Sp from this source term on the right hand side. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Przemek
Join Date: Jun 2011
Posts: 249
Rep Power: 16 ![]() |
__________________
best regards pblasiak |
|
![]() |
![]() |
![]() |
Tags |
alphaeqn, mules, source term |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Using MULES to solve transport equation with sharp interface | derkermit | OpenFOAM Programming & Development | 3 | February 27, 2014 06:37 |
interFoam: Solve gives better results than MULES | hrushi.397 | OpenFOAM Running, Solving & CFD | 6 | September 10, 2013 12:57 |
How to solve energy eqn with enthalpy terms? | sabitha | FLUENT | 0 | April 6, 2006 04:58 |
Alpha to be phased out | Heiko Gerhauser | Main CFD Forum | 17 | July 25, 2001 14:40 |
Determining alpha and beta for porous baffle | Liaqat Khan | Siemens | 1 | October 27, 2000 04:44 |