00001 #ifndef MULTIPHASE_H
00002 #define MULTIPHASE_H
00003
00004 #include "genfile.h"
00005
00006 class multiph
00007 {
00008 public:
00009 multiph();
00010 ~multiph();
00011
00012 void matcond (matrix &d,long ri,long ci,long ipp);
00013 void matcap (double &c,long ri,long ci,long ipp);
00014
00015 void matcond1d (matrix &d,long ri,long ci,long ipp);
00016 void matcond2d (matrix &d,long ri,long ci,long ipp);
00017 void matcond3d (matrix &d,long ri,long ci,long ipp);
00018
00019 void matcond2 (matrix &d,long ri,long ci,long ipp);
00020
00021 void matcond1d_2 (matrix &d,long ri,long ci,long ipp);
00022 void matcond2d_2 (matrix &d,long ri,long ci,long ipp);
00023 void matcond3d_2 (matrix &d,long ri,long ci,long ipp);
00024
00025 void rhs_volume (matrix &d,long ri,long ci,long ipp);
00026 void rhs1d1 (matrix &d,long ri,long ci,long ipp);
00027 void rhs2d1 (matrix &d,long ri,long ci,long ipp);
00028 void rhs3d1 (matrix &d,long ri,long ci,long ipp);
00029
00030 double transmission_transcoeff(double trc,long ri,long ci,long nn,long bc,long ipp);
00031 double transmission_transcoeff(double trc,long ri,long ci,long nn,long bc,long ipp,int flag);
00032 double transmission_nodval(double nodval,double trc2,long ri,long ci,long nn,long bc,long ipp);
00033 double transmission_flux(double nodval,double trc2,long ri,long ci,long nn,long bc,long ipp);
00034
00035 void gaspress_check(double pc,double &pg,double t,long ipp);
00036 void cappress_check(double &pc,double pg,double t,long ipp);
00037 void cappress_stop(double &pc,double pg,double t,long ipp);
00038
00039 void values_correction (vector &nv);
00040
00041 double get_kcc(double pc,double pg,double t,long ipp);
00042 double get_capcc(double pc,double pg,double t,long ipp);
00043 double get_kcg(double pc,double pg,double t,long ipp);
00044 double get_capcg(double pc,double pg,double t,long ipp);
00045 double get_kct(double pc,double pg,double t,long ipp);
00046 double get_capct(double pc,double pg,double t,long ipp);
00047
00048 double get_kgg(double pc,double pg,double t,long ipp);
00049 double get_capgg(double pc,double pg,double t,long ipp);
00050 double get_kgc(double pc,double pg,double t,long ipp);
00051 double get_capgc(double pc,double pg,double t,long ipp);
00052 double get_kgt(double pc,double pg,double t,long ipp);
00053 double get_capgt(double pc,double pg,double t,long ipp);
00054
00055 double get_ktt1(double pc,double pg,double t,long ipp);
00056 double get_ktt2(double pc,double pg,double t,long ipp);
00057 double get_captt(double pc,double pg,double t,long ipp);
00058 double get_ktg(double pc,double pg,double t,long ipp);
00059 double get_captg(double pc,double pg,double t,long ipp);
00060 double get_ktc(double pc,double pg,double t,long ipp);
00061 double get_captc(double pc,double pg,double t,long ipp);
00062
00063 double get_ktt2a(double pc,double pg,double t,long ipp);
00064 double get_ktt2b(double pc,double pg,double t,long ipp);
00065 double get_ktt2c(double pc,double pg,double t,long ipp);
00066 double get_ktt2d(double pc,double pg,double t,long ipp);
00067
00068 double get_fc1(double pc,double pg,double t,long ipp);
00069 double get_fg(double pc,double pg,double t,long ipp);
00070 double get_ft1(double pc,double pg,double t,long ipp);
00071
00072 double get_transmission_transcoeff_cc(double pc,double pg,double t,long bc,long ipp);
00073 double get_transmission_transcoeff_cc(double pc,double pg,double t,long bc,long ipp,int flag);
00074 double get_transmission_nodval_cc(double bv,double pc,double pg,double t,long bc,long ipp);
00075 double get_transmission_flux_cc(double bv,double pc,double pg,double t,long bc,long ipp);
00076
00077 double get_transmission_transcoeff_tt(double pc,double pg,double t,long bc,long ipp);
00078 double get_transmission_nodval_tt(double bv,double trr,double pc,double pg,double t,long bc,long ipp);
00079 double get_transmission_flux_tt(double bv,double trr,double pc,double pg,double t,long bc,long ipp);
00080
00081 double heat_rate(double stime,double time);
00082 double iso_fire (double time, double t0, double tfirestart);
00083
00084 double get_othervalue(long compother,long ipp,double *r);
00085 void print_othervalue_name(FILE *out,long compother);
00086
00087 private:
00088
00089 double ma;
00090 double mw;
00091 double gasr;
00092
00093 double scale_pc,scale_pg,scale_t;
00094 };
00095
00096 #endif