|
[Sponsors] |
November 19, 2015, 05:00 |
Y1iEqn in reactingTwoPhaseEulerFoam
|
#1 |
Member
Victor Koppejan
Join Date: May 2015
Posts: 40
Rep Power: 10 |
Hi everyone,
I'm playing around with the new reactingTwoPhaseEulerFoam solver. I'd like to adapt to include laminar diffusivities per species. After looking around I think I now have found several options for this 1. Specify a Lewis number for each specie. 2. Rewrite the species transp. eqn to include a mueff / Sc number, and define Sc for each specie. 3. Rewrite the species transport eqn to laminar form (as in scalarTransportFoam) and add define diffusivities for each specie. One of the problems I have is that I can't seem to find the species transport eqn. In the YEqns.H file it is defined as: Code:
forAll(Y1, i) { tmp<fvScalarMatrix> Y1iEqn(phase1.YiEqn(Y1[i])); if (Y1iEqn.valid()) { Y1iEqn = ( Y1iEqn == *massTransfer[Y1[i].name()] + fvOptions(alpha1, rho1, Y1[i]) ); Y1iEqn->relax(); Y1iEqn->solve(mesh.solver("Yi")); } } Btw, I know the Le number is correctly used in all the masstransfer models (apart from a few typos), but I really want to include species diffusivity in the Laplacian term. Thanks for your help. Victor |
|
November 21, 2015, 02:24 |
|
#2 |
Member
Mattia de\' Michieli Vitturi
Join Date: Mar 2009
Posts: 50
Rep Power: 17 |
Hi Victor,
you can find the equation here: applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseModel/MultiComponentPhaseModel/MultiComponentPhaseModel.C Code:
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) |
|
January 16, 2020, 06:35 |
3 terms of the scalar transport eq
|
#3 | |
Member
Rosario Arnau
Join Date: Feb 2017
Location: Spain
Posts: 57
Rep Power: 9 |
Hi Foamers,
I took a look at the transport equation and I'm confused about some terms of it. Quote:
1. Code:
- fvm::Sp(this->continuityError(), Yi) Code:
+ fvc::ddt(residualAlpha_*rho, Yi) Code:
- fvm::ddt(residualAlpha_*rho, Yi) Thanks, |
||
|
|