31 #include "TCCExtremoAnecoico.h"
37 TCCExtremoAnecoico::TCCExtremoAnecoico(nmTypeBC TipoCC,
int numCC, nmTipoCalculoEspecies SpeciesModel,
38 int numeroespecies, nmCalculoGamma GammaCalculation,
bool ThereIsEGR) :
39 TCondicionContorno(TipoCC, numCC, SpeciesModel, numeroespecies, GammaCalculation, ThereIsEGR) {
47 TCCExtremoAnecoico::~TCCExtremoAnecoico() {
49 delete[] FTuboExtremo;
56 void TCCExtremoAnecoico::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 FCC = &(FTuboExtremo[FNumeroTubosCC].Beta);
67 FCD = &(FTuboExtremo[FNumeroTubosCC].Landa);
68 FTuboExtremo[FNumeroTubosCC].TipoExtremo = nmLeft;
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: TCCExtremoAnecoico::AsignaTubos en la condicion de contorno: " << FNumeroCC << std::endl;
93 std::cout <<
"Tipo de error: " << N.what() << std::endl;
101 void TCCExtremoAnecoico::CalculaCondicionContorno(
double Time) {
103 double FraccionMasicaAcum = 0.;
105 *FCD = FTuboExtremo[0].Entropia;
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);
116 }
catch(exception &N) {
117 std::cout <<
"ERROR: TCCExtremoAnecoico::CalculaCondicionesContorno en la condicion de contorno: " << FNumeroCC <<
119 std::cout <<
"Tipo de error: " << N.what() << std::endl;
127 #pragma package(smart_init)