31 #include "TCCExtremoCerrado.h"
37 TCCExtremoCerrado::TCCExtremoCerrado(nmTypeBC TipoCC,
int numCC, nmTipoCalculoEspecies SpeciesModel,
int numeroespecies,
38 nmCalculoGamma GammaCalculation,
bool ThereIsEGR) :
39 TCondicionContorno(TipoCC, numCC, SpeciesModel, numeroespecies, GammaCalculation, ThereIsEGR) {
47 TCCExtremoCerrado::~TCCExtremoCerrado() {
49 delete[] FTuboExtremo;
56 void TCCExtremoCerrado::AsignaTubos(
int NumberOfPipes,
TTubo **Pipe) {
61 FTuboExtremo[0].Pipe = NULL;
63 while(FNumeroTubosCC < 1 && i < NumberOfPipes) {
64 if(Pipe[i]->getNodoIzq() == FNumeroCC) {
65 FTuboExtremo[FNumeroTubosCC].Pipe = Pipe[i];
66 FTuboExtremo[FNumeroTubosCC].TipoExtremo = nmLeft;
67 FCC = &(FTuboExtremo[FNumeroTubosCC].Beta);
68 FCD = &(FTuboExtremo[FNumeroTubosCC].Landa);
73 if(Pipe[i]->getNodoDer() == FNumeroCC) {
74 FTuboExtremo[FNumeroTubosCC].Pipe = Pipe[i];
75 FTuboExtremo[FNumeroTubosCC].TipoExtremo = nmRight;
76 FCC = &(FTuboExtremo[FNumeroTubosCC].Landa);
77 FCD = &(FTuboExtremo[FNumeroTubosCC].Beta);
78 FNodoFin = FTuboExtremo[FNumeroTubosCC].Pipe->getNin() - 1;
86 FFraccionMasicaEspecie =
new double[FNumeroEspecies - FIntEGR];
87 for(
int i = 0; i < FNumeroEspecies - FIntEGR; i++) {
91 }
catch(exception &N) {
92 std::cout <<
"ERROR: TCCExtremoCerrado::AsignaTubos en la condicion de contorno: " << FNumeroCC << std::endl;
93 std::cout <<
"Tipo de error: " << N.what() << std::endl;
101 void TCCExtremoCerrado::CalculaCondicionContorno(
double Time) {
103 double FraccionMasicaAcum = 0.;
108 for(
int j = 0; j < FNumeroEspecies - 2; j++) {
110 FraccionMasicaAcum += FFraccionMasicaEspecie[j];
112 FFraccionMasicaEspecie[FNumeroEspecies - 2] = 1. - FraccionMasicaAcum;
114 FFraccionMasicaEspecie[FNumeroEspecies - 1] = FTuboExtremo[0].Pipe->
GetFraccionMasicaCC(FIndiceCC, FNumeroEspecies - 1);
115 }
catch(exception &N) {
116 std::cout <<
"ERROR: TCCExtremoCerrado::CalculaCondicionesContorno en la condicion de contorno: " << FNumeroCC <<
118 std::cout <<
"Tipo de error: " << N.what() << std::endl;
126 #pragma package(smart_init)