sutherland_viscosity_unit.C File Reference

#include <iostream>
#include <cmath>
#include <limits>
#include "antioch/sutherland_viscosity.h"

Go to the source code of this file.

Functions
template<typename Scalar >
int	test_viscosity (const Scalar mu, const Scalar mu_exact, const Scalar tol)
template<typename Scalar >
int	tester ()
int	main ()

Function Documentation

int main

(

)

Definition at line 99 of file sutherland_viscosity_unit.C.

{
  return (tester<double>() ||
          tester<long double>() ||
          tester<float>());
}

template<typename Scalar >

int test_viscosity	(	const Scalar	mu,
		const Scalar	mu_exact,
		const Scalar	tol
	)

Definition at line 42 of file sutherland_viscosity_unit.C.

Referenced by tester().

{
  using std::abs;

  int return_flag = 0;

  const Scalar rel_error = abs( (mu - mu_exact)/mu_exact);

  if( rel_error  > tol )
    {
      std::cerr << "Error: Mismatch in viscosity" << std::endl
                << "mu(T)    = " << mu << std::endl
                << "mu_exact = " << mu_exact << std::endl
                << "rel_error = " << rel_error << std::endl
                << "tol = " << tol << std::endl;
      return_flag = 1;
    }

  return return_flag;
}

template<typename Scalar >

int tester

(

)

Definition at line 65 of file sutherland_viscosity_unit.C.

References Antioch::SutherlandViscosity< CoeffType >::reset_coeffs(), and test_viscosity().

{
  const Scalar mu_ref = 1.0e-3L;
  const Scalar T_ref = 300.0L;

  Antioch::SutherlandViscosity<Scalar> mu(mu_ref,T_ref);

  std::cout << mu << std::endl;

  const Scalar T = 1521.2L;

  // bc with scale=40 gives
  const Scalar mu_exact = .0325778060534850406481862157435995107036L;

  int return_flag = 0;

  const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 10;

  return_flag = test_viscosity( mu(T), mu_exact, tol );
  
  const Scalar mu_ref2 = 3.14159e-3L;
  const Scalar T_ref2 = 420.42L;

  mu.reset_coeffs(mu_ref2,T_ref2);

  // bc with scale=40 gives
  const Scalar mu_exact2 = .0959985656417205050367745642313443587197L;

  return_flag = test_viscosity( mu(T), mu_exact2, tol );

  return return_flag;
}