@@ -325,7 +325,7 @@ mod tests {
325325 use approx:: assert_relative_eq;
326326 use feos_core:: { Contributions , DensityInitialization , PhaseEquilibrium , State } ;
327327 use ndarray:: arr1;
328- use quantity:: si:: { BAR , KELVIN , METER , PASCAL , RGAS , SECOND } ;
328+ use quantity:: si:: { BAR , KELVIN , METER , PASCAL , RGAS } ;
329329
330330 #[ test]
331331 fn ideal_gas_pressure ( ) {
@@ -378,7 +378,7 @@ mod tests {
378378 fn vle_pure_t ( ) {
379379 let e = Rc :: new ( Pets :: new ( argon_parameters ( ) ) ) ;
380380 let t = 300.0 * KELVIN ;
381- let vle = PhaseEquilibrium :: pure_t ( & e, t, None , Default :: default ( ) ) ;
381+ let vle = PhaseEquilibrium :: pure ( & e, t, None , Default :: default ( ) ) ;
382382 if let Ok ( v) = vle {
383383 assert_relative_eq ! (
384384 v. vapor( ) . pressure( Contributions :: Total ) ,
@@ -388,30 +388,6 @@ mod tests {
388388 }
389389 }
390390
391- #[ test]
392- fn critical_point ( ) {
393- let e = Rc :: new ( Pets :: new ( argon_parameters ( ) ) ) ;
394- let t = 300.0 * KELVIN ;
395- let cp = State :: critical_point ( & e, None , Some ( t) , Default :: default ( ) ) ;
396- if let Ok ( v) = cp {
397- assert_relative_eq ! ( v. temperature, 375.1244078318015 * KELVIN , epsilon = 1e-8 )
398- }
399- }
400-
401- #[ test]
402- fn speed_of_sound ( ) {
403- let e = Rc :: new ( Pets :: new ( argon_parameters ( ) ) ) ;
404- let t = 300.0 * KELVIN ;
405- let p = BAR ;
406- let m = arr1 ( & [ 1.0 ] ) * MOL ;
407- let s = State :: new_npt ( & e, t, p, & m, DensityInitialization :: None ) . unwrap ( ) ;
408- assert_relative_eq ! (
409- s. speed_of_sound( ) ,
410- 245.00185709137546 * METER / SECOND ,
411- epsilon = 1e-4
412- )
413- }
414-
415391 #[ test]
416392 fn mix_single ( ) {
417393 let e1 = Rc :: new ( Pets :: new ( argon_parameters ( ) ) ) ;
@@ -437,48 +413,4 @@ mod tests {
437413 epsilon = 1e-12
438414 )
439415 }
440-
441- fn viscosity ( ) -> EosResult < ( ) > {
442- let e = Rc :: new ( Pets :: new ( argon_parameters ( ) ) ) ;
443- let t = 300.0 * KELVIN ;
444- let p = BAR ;
445- let n = arr1 ( & [ 1.0 ] ) * MOL ;
446- let s = State :: new_npt ( & e, t, p, & n, DensityInitialization :: None ) . unwrap ( ) ;
447- assert_relative_eq ! (
448- s. viscosity( ) ?,
449- 0.00797 * MILLI * PASCAL * SECOND ,
450- epsilon = 1e-5
451- ) ;
452- assert_relative_eq ! (
453- s. ln_viscosity_reduced( ) ?,
454- ( s. viscosity( ) ? / e. viscosity_reference( s. temperature, s. volume, & s. moles) ?)
455- . into_value( )
456- . unwrap( )
457- . ln( ) ,
458- epsilon = 1e-15
459- ) ;
460- Ok ( ( ) )
461- }
462-
463- fn diffusion ( ) -> EosResult < ( ) > {
464- let e = Rc :: new ( Pets :: new ( argon_parameters ( ) ) ) ;
465- let t = 300.0 * KELVIN ;
466- let p = BAR ;
467- let n = arr1 ( & [ 1.0 ] ) * MOL ;
468- let s = State :: new_npt ( & e, t, p, & n, DensityInitialization :: None ) . unwrap ( ) ;
469- assert_relative_eq ! (
470- s. diffusion( ) ?,
471- 0.01505 * ( CENTI * METER ) . powi( 2 ) / SECOND ,
472- epsilon = 1e-5
473- ) ;
474- assert_relative_eq ! (
475- s. ln_diffusion_reduced( ) ?,
476- ( s. diffusion( ) ? / e. diffusion_reference( s. temperature, s. volume, & s. moles) ?)
477- . into_value( )
478- . unwrap( )
479- . ln( ) ,
480- epsilon = 1e-15
481- ) ;
482- Ok ( ( ) )
483- }
484416}
0 commit comments