/[ascend]/trunk/models/johnpye/fprops/fluids/nitrogen.c
ViewVC logotype

Diff of /trunk/models/johnpye/fprops/fluids/nitrogen.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2758 by jpye, Fri Mar 21 08:00:18 2014 UTC revision 2759 by jpye, Fri Mar 21 11:25:55 2014 UTC
# Line 171  const ThermalConductivityData thcond_nit Line 171  const ThermalConductivityData thcond_nit
171          .k_star = 1e-3          .k_star = 1e-3
172          ,.T_star = 126.192          ,.T_star = 126.192
173          ,.rho_star = 11.1839          ,.rho_star = 11.1839
174          ,.v1=&(visc_carbondioxide.data.v1)          ,.v1=&(visc_nitrogen.data.v1)
175          ,.eps_over_k = 98.94          ,.eps_over_k = 98.94
176          ,.nc = 6          ,.nc = 6
177          ,.ct=(const ThCondCSTerm[]){          ,.ct=(const ThCondCSTerm[]){
# Line 205  EosData eos_nitrogen = { Line 205  EosData eos_nitrogen = {
205      ,FPROPS_HELMHOLTZ      ,FPROPS_HELMHOLTZ
206      ,.data = {.helm = &helmholtz_data_nitrogen}      ,.data = {.helm = &helmholtz_data_nitrogen}
207      ,.visc = &visc_nitrogen      ,.visc = &visc_nitrogen
208      ,.thcond = &thcond_nitrogen;      ,.thcond = &thcond_nitrogen
209  };  };
210    
211    
# Line 221  EosData eos_nitrogen = { Line 221  EosData eos_nitrogen = {
221    
222  #include "../ideal_impl.h"  #include "../ideal_impl.h"
223  #include "../visc.h"  #include "../visc.h"
224    #include "../thcond.h"
225  #include "../test.h"  #include "../test.h"
226  #include <math.h>  #include <math.h>
227  #include <stdio.h>  #include <stdio.h>
# Line 310  fprintf(stderr,"mu(T=%f, rho=%f) = %e (t Line 311  fprintf(stderr,"mu(T=%f, rho=%f) = %e (t
311    
312      fprintf(stderr,"done\n");      fprintf(stderr,"done\n");
313    
314        //--------------------------------------------------------------------------
315        fprintf(stderr,"Testing thermal conductivity values from REFPROP 8.0\n");
316        thcond_prepare(P, &thcond_nitrogen, &err);
317        ASSERT(FPROPS_NO_ERROR==err);
318        ASSERT(V != NULL);
319    
320        double k;
321    #define THCOND_TEST(T__1,RHO__1,K__1,TOL__1) \
322        S = fprops_set_Trho(T__1, RHO__1, P, &err); \
323        k = fprops_k(S,&err); \
324        fprintf(stderr,"k(T=%f, rho=%f) = %e (target: %e)\n",S.T,S.rho,k,K__1); \
325        ASSERT(FPROPS_NO_ERROR==err); \
326        ASSERT(fabs(k - K__1)<TOL__1);
327    
328        THCOND_TEST(100,0,          9.2775e-3, 0.00005e-3);
329        THCOND_TEST(300,0,          25.936e-3, 0.0005e-3);
330        THCOND_TEST(100,25,         10.309e-3, 0.0005e-3);
331        THCOND_TEST(200,10,         18.545e-3, 0.0005e-3);
332        THCOND_TEST(300,5,          26.085e-3, 0.0005e-3);
333        THCOND_TEST(126.195,11.180, 12.132e-3, 0.0005e-3);
334    
335        //--------------------------------------------------------------------------
336      fprintf(stderr,"CONSISTENCY TESTS (of test data): u, T, s, a... ");      fprintf(stderr,"CONSISTENCY TESTS (of test data): u, T, s, a... ");
337      for(i=0; i<n; ++i){      for(i=0; i<n; ++i){
338          u = td[i].u*1e3;          u = td[i].u*1e3;

Legend:
Removed from v.2758  
changed lines
  Added in v.2759

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22