61 double FLongitudTotal;
64 nmTipoMallado FTipoMallado;
67 double *FDiametroTubo;
75 double *FDerLinArea12;
81 nmTipoTransCal FTipoTransCal;
84 double *FCoefTurbulencia;
152 nmTipoCalculoEspecies FCalculoEspecies;
153 nmCalculoGamma FCalculoGamma;
154 double *FComposicionInicial;
155 double **FFraccionMasicaEspecie;
156 double **FFraccionMasicaEspecie1;
157 double **FFraccionMasicaCC;
164 double *FVelocidadCC;
183 double *FCourantLocal;
192 double FTiempoMedSUM;
193 double FControlResMed;
195 int FNumDistSensores;
197 double FAnguloTotalCiclo;
198 double FRegimenFicticio;
205 int FNumeroParesdeCanales;
206 nmTipoCanal FTipoCanal;
209 int FTipoSeccionCanal;
210 nmSeccionCanal FSeccionCanal;
211 double *FLadoCanalTramo;
212 double *FLadoMenorCanalTramo;
213 double *FLadoMayorCanalTramo;
216 double *FLadoCanalD12;
217 double *FLadoCanalS12;
221 double *FVelocidadPared;
222 double *FViscosidadDinamica;
223 double *FEspesorSoot;
225 int FNodoInicialFuente;
226 int FNodoFinalFuente;
227 int FNodoInicialFuenteCE;
228 double FDistanciaInterpolacion;
233 double **FTasaFraccionMasicaEspecie;
234 double **FFraccionMasicaSalida;
238 double *FSupEspecifica;
241 double FBetaExtremoCerrado;
242 double FLandaExtremoCerrado;
243 double FEntropiaExtremoCerrado;
251 void CalculoPuntosMalla();
253 double InterpolaTubo(
double vizq,
double vder,
double axid,
double xif);
255 void Transforma1Area(
double v,
double a,
double p,
double **U,
double area,
double Gamma,
double Gamma1,
256 double *Yespecie,
int i);
258 void Transforma2Area(
double *v,
double *a,
double *p,
double **U,
double Area,
double Gamma,
double Gamma1,
259 double *Yespecie,
int i);
261 void Transforma3Area(
double **Ufct,
double **U,
double Area,
double Gamma,
double Gamma1,
double Gamma6,
int i);
263 void Transforma4Area(
double **U1,
double **Ufctd,
double Area,
double Gamma,
double Gamma1,
double Gamma3,
264 double Gamma4,
double Gamma6,
int i);
266 double Maximo(
double x,
double y);
268 double Minimo(
double x,
double y);
270 void LaxWendroffArea();
272 void FluxCorrectedTransport();
274 void CalculaFlujo(
double **U,
double **W,
double *Gamma,
double *Gamma1,
int Nodos);
276 void CalculaFuente1(
double **U,
double **V1,
double *Gamma,
double *Gamma1,
int Nodos);
278 void CalculaFuente1Area(
double **U,
double **V1,
double *Area,
double *Gamma1,
int Nodos);
280 void CalculaFuente2(
double **U,
double **V2,
double *diame,
double *CoefTurbulencia,
double *hi,
double *rho,
281 double *TempParedTubo,
double *Gamma,
double *Rmezcla,
double *Gamma1,
double *LadoCanal,
double *VelocidadPared,
282 double *ViscosidadDinamica,
double *q_reac1,
double *q_reac2,
double **FraccionMasicaSalida,
double *H0Pared,
int Nodos,
283 double *EspesorSoot,
double *Eficiencia,
double *Rreg1,
double *Rreg2,
double *SupEspecifica,
double *LongitudVC,
284 double **FraccionMasica);
286 void CalculaFuente2Area(
double **U,
double **V2,
double *diame,
double *CoefTurbulencia,
double *hi,
double *rho,
287 double *TempParedTubo,
double *Gamma,
double *Rmezcla,
double *Gamma1,
double *LadoCanal,
double *VelocidadPared,
288 double *ViscosidadDinamica,
double *q_reac1,
double *q_reac2,
double **FraccionMasicaSalida,
double *H0Pared,
int Nodos,
289 double *Area,
double *EspesorSoot,
double *Eficiencia,
double *Rreg1,
double *Rreg2,
double *SupEspecifica,
290 double *LongitudVC,
double **FraccionMasica);
292 double DerLinF(
double d1,
double d2,
double xref);
294 double DerLinFArea(
double area1,
double area2,
double xref);
296 void Colebrook(
double rug,
double dia,
double *f,
double Re);
298 void TransmisionCalor(
double tgas,
double lado,
double *q,
double CoefTurbulencia,
double hi,
double rho,
double Tw);
300 void TransformaContorno(
double *L,
double *B,
double *E,
double *a,
double *v,
double *p,
int modo,
double Gamma1,
301 double Gamma3,
double Gamma4,
double Gamma5);
303 double CalculaNIT(
double a,
double v,
double p,
double d,
double Gamma,
double R);
307 void Calculo_Entropia(
double *entropia,
double *velocidadp,
int ind,
double dist,
int signo,
double DeltaTiempo,
308 int indiceCC, nmPipeEnd TipoExtremoTubo);
310 void Calculo_Caracteristica(
double *caracteristica,
double *velocidadp,
double *asonidop,
int ind,
double dist,
311 int signo,
double entropia,
double DeltaTiempo, nmPipeEnd TipoExtremoTubo);
319 void CalculaMatrizJacobiana();
321 void TVD_Estabilidad();
323 void TVD_Limitador();
325 void DimensionaTVD();
329 double Limita(
double r);
331 void NodoTerminoFuente();
342 TCanalDPF(
int NumeroEspecies,
int j, nmTipoCalculoEspecies CalculoEspecies, nmCalculoGamma CalculoGamma,
bool HayEGR,
343 TDPF *DPF,
int numerocanal,
int NumeroDPF);
347 void LeeDatosGeneralesCanal(
const char *FileWAM, fpos_t &filepos,
int NodoIzquierdo,
int NodoDerecho);
349 void LeeDatosGeometricosCanal(
const char *FileWAM, fpos_t &filepos);
351 void IniciaVariablesFundamentalesCanalDPF();
353 void EstabilidadMetodoCalculo();
355 void CalculaVariablesFundamentales();
359 void ReduccionFlujoSubsonico();
361 void ReduccionFlujoSubsonicoFCT();
363 void LeeResultadosMediosCanalDPF(
int NumResMedios,
char *FileWAM, fpos_t &filepos);
365 void CabeceraResultadosMedios(std::stringstream& medoutput,
stEspecies *DatosEspecies)
const;
367 void ImprimeResultadosMedios(std::stringstream& medoutput)
const;
369 void LeeResultadosInstantaneosCanalDPF(
int NumResInstantaneos,
char *FileWAM, fpos_t &filepos);
371 void CabeceraResultadosInstantaneos(std::stringstream& insoutput,
stEspecies *DatosEspecies)
const;
373 void ImprimeResultadosInstantaneos(std::stringstream& insoutput)
const;
375 void CalculaResultadosMedios(
double theta);
377 void CalculaResultadosInstantaneos();
379 void SalidaGeneralTubos(
stEspecies *DatosEspecies)
const;
381 void AjustaPaso(
double Intervalo);
385 double Interpola_Entropia(nmPipeEnd TipoExtremoTubo,
double DeltaTiempo);
387 double Interpola_Caracteristica(
double entropia,
int signo,
int extremo,
double DeltaTiempo, nmPipeEnd TipoExtremoTubo);
393 void ActualizaPropiedadesGas();
395 void CalculaCoeficientePeliculaExterior(
double pamb);
397 void CalculaCoeficientePeliculaInterior();
399 int GetNumeroDeposito() {
400 return FNumeroDeposito;
403 double GetDiametro(
int i) {
404 return FDiametroTubo[i];
407 double GetFraccionMasicaInicial(
int i) {
408 return FComposicionInicial[i];
411 double GetFraccionMasicaCC(
int j,
int i) {
412 return FFraccionMasicaCC[j][i];
415 double GetViscosidadDinamica(
int i) {
416 return FViscosidadDinamica[i];
423 double getNodoIzq() {
427 double getNodoDer() {
431 void putTime1(
double valor) {
438 void putTime0(
double valor) {
442 double getDeltaTime() {
445 void putDeltaTime(
double valor) {
457 double getLongitudTotal() {
458 return FLongitudTotal;
461 nmTipoTransCal getTipoTransCal() {
462 return FTipoTransCal;
465 double getCoefAjustFric() {
466 return FCoefAjusFric;
469 int getNumeroCanales() {
470 return FNumeroParesdeCanales;
473 double getMallado() {
478 double GetDensidad(
int i) {
481 double GetArea(
int i) {
484 double GetVelocidad(
int i) {
485 return FVelocidad0[i];
488 double GetAsonido(
int i) {
492 double GetPresion(
int i) {
496 double GetCoefTurbulencia(
int i) {
497 return FCoefTurbulencia[i];
500 double getTemperaturaInicial() {
504 double getPresionInicial() {
508 double getVelocidadMedia() {
512 double GetYespecie(
int j,
int i) {
513 return FFraccionMasicaEspecie[j][i];
516 double GetYespecieCC(
int j,
int i) {
517 return FFraccionMasicaCC[j][i];
520 double GetVelocidadCC(
int i) {
521 return FVelocidadCC[i];
524 double GetDensidadCC(
int i) {
525 return FDensidadCC[i];
528 double GetAreaCC(
int i) {
532 double GetGamma(
int i) {
536 double GetRMezcla(
int i) {
540 double GetCp(
int i) {
544 double GetCv(
int i) {
548 int getNumeroEspecies() {
549 return FNumeroEspecies;
556 double Gethext(
int i) {
560 double GetH0Pared(
int i) {
565 nmFormulacionLeyes getFormulacionLeyes() {
566 return FMod.FormulacionLeyes;
569 double GetU0(
int i,
int j) {
573 int getNodoInicialFuente() {
574 return FNodoInicialFuente;
577 int getNodoFinalFuente() {
578 return FNodoFinalFuente;
581 double getDistanciaInterpolacion() {
582 return FDistanciaInterpolacion;
585 double GetLadoCanal(
int i) {
586 return FLadoCanal[i];