antioch_wilke_evaluator_regression.C File Reference

#include "grins_config.h"
#include <iomanip>
#include <limits>
#include <vector>
#include "grins/antioch_wilke_transport_mixture.h"
#include "grins/antioch_wilke_transport_evaluator.h"
#include "grins/cached_values.h"
#include "libmesh/getpot.h"

Include dependency graph for antioch_wilke_evaluator_regression.C:

Go to the source code of this file.

Functions
int	test_generic (const libMesh::Real value, const libMesh::Real value_reg, const std::string &name)
template<typename Viscosity >
int	test_mu (const libMesh::Real mu)
template<typename Thermo , typename Conductivity >
int	test_k (const libMesh::Real k)
template<typename Thermo , typename Viscosity , typename Conductivity , typename Diffusivity >
int	test_D (const std::vector< libMesh::Real > &D)
template<>
int	test_mu< Antioch::MixtureViscosity< Antioch::BlottnerViscosity< libMesh::Real > > > (const libMesh::Real mu)
template<typename Thermo , typename Viscosity , typename Conductivity , typename Diffusivity >
int	test_evaluator (const GetPot &input)
int	main (int argc, char *argv[])

Function Documentation

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 191 of file antioch_wilke_evaluator_regression.C.

{
  // Check command line count.
  if( argc < 2 )
    {
      // TODO: Need more consistent error handling.
      std::cerr << "Error: Must specify input file." << std::endl;
      exit(1);
    }

  GetPot input( argv[1] );

  int return_flag = 0;

  std::cout << std::endl <<  "Running StatMesh, Blottner, Eucken, Constant Lewis regression test." << std::endl;
  return_flag = test_evaluator<Antioch::StatMechThermodynamics<libMesh::Real>,
    Antioch::MixtureViscosity<Antioch::BlottnerViscosity<libMesh::Real> >,
    Antioch::EuckenThermalConductivity<Antioch::StatMechThermodynamics<libMesh::Real> >,
    Antioch::ConstantLewisDiffusivity<libMesh::Real> >(input);

  return return_flag;
}

template<typename Thermo , typename Viscosity , typename Conductivity , typename Diffusivity >

int test_D

(

const std::vector< libMesh::Real > &

D

)

template<typename Thermo , typename Viscosity , typename Conductivity , typename Diffusivity >

int test_evaluator

(

const GetPot &

input

)

Definition at line 115 of file antioch_wilke_evaluator_regression.C.

References GRINS::CachedValues::add_quantity(), GRINS::AntiochWilkeTransportEvaluator< Thermo, Viscosity, Conductivity, Diffusivity >::D(), GRINS::AntiochWilkeTransportEvaluator< Thermo, Viscosity, Conductivity, Diffusivity >::k(), GRINS::Cache::MASS_FRACTIONS, GRINS::Cache::MIXTURE_DENSITY, GRINS::AntiochWilkeTransportEvaluator< Thermo, Viscosity, Conductivity, Diffusivity >::mu(), GRINS::AntiochWilkeTransportEvaluator< Thermo, Viscosity, Conductivity, Diffusivity >::mu_and_k(), GRINS::CachedValues::set_values(), GRINS::CachedValues::set_vector_values(), GRINS::AntiochChemistry::species_name(), and GRINS::Cache::TEMPERATURE.

{
  GRINS::AntiochWilkeTransportMixture<Thermo,Viscosity,Conductivity,Diffusivity> mixture(input);

  GRINS::AntiochWilkeTransportEvaluator<Thermo,Viscosity,Conductivity,Diffusivity> evaluator(mixture);

  const libMesh::Real T = 1000;

  const libMesh::Real rho = 1.0e-3;

  const unsigned int n_species = 5;

  std::vector<libMesh::Real> Y(n_species,0.2);

  GRINS::CachedValues cache;

  cache.add_quantity(GRINS::Cache::TEMPERATURE);
  std::vector<double> Tqp(1,T);
  cache.set_values(GRINS::Cache::TEMPERATURE, Tqp);

  cache.add_quantity(GRINS::Cache::MIXTURE_DENSITY);
  std::vector<double> rhoqp(1,rho);
  cache.set_values(GRINS::Cache::MIXTURE_DENSITY, rhoqp);

  cache.add_quantity(GRINS::Cache::MASS_FRACTIONS);
  std::vector<std::vector<double> > Yqp(1,Y);
  cache.set_vector_values(GRINS::Cache::MASS_FRACTIONS, Yqp);

  libMesh::Real mu = evaluator.mu( cache, 0 );

  libMesh::Real k = evaluator.k( cache, 0 );

  libMesh::Real mu2 = 0.0;
  libMesh::Real k2 = 0.0;

  evaluator.mu_and_k( cache, 0, mu2, k2 );

  std::vector<libMesh::Real> D(n_species,0.0);

  evaluator.D( cache, 0, D );

  std::cout << std::scientific << std::setprecision(16)
            << "mu = " << mu << std::endl;

  std::cout << std::scientific << std::setprecision(16)
            << "k = " << k << std::endl;

  for( unsigned int i = 0; i < n_species; i++ )
    {
      std::cout << std::scientific << std::setprecision(16) 
                << "D(" << mixture.species_name(i) << ") = " << D [i] << std::endl;
    }

  int return_flag = 0;

  int return_flag_temp = 0;

  return_flag_temp = test_mu<Viscosity>( mu );
  if( return_flag_temp != 0 ) return_flag = 1;

  return_flag_temp = test_k<Thermo,Conductivity>( k );
  if( return_flag_temp != 0 ) return_flag = 1;

  return_flag_temp = test_mu<Viscosity>( mu2 );
  if( return_flag_temp != 0 ) return_flag = 1;

  return_flag_temp = test_k<Thermo,Conductivity>( k2 );
  if( return_flag_temp != 0 ) return_flag = 1;
  
  return_flag_temp = test_D<Thermo,Viscosity,Conductivity,Diffusivity>( D );
  if( return_flag_temp != 0 ) return_flag = 1;
 
  return return_flag;
}

int test_generic	(	const libMesh::Real	value,
		const libMesh::Real	value_reg,
		const std::string &	name
	)

Definition at line 43 of file antioch_wilke_evaluator_regression.C.

Referenced by test_mu< Antioch::MixtureViscosity< Antioch::BlottnerViscosity< libMesh::Real > > >().

{
  int return_flag = 0;

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

  const double rel_error = std::fabs( (value - value_reg)/value_reg );

  if( rel_error > tol )
    {
      return_flag = 1;
      std::cout << "Mismatch in "+name  << std::endl
                << name+" = " << value << std::endl
                << name+"_reg = " << value_reg << std::endl
                << "rel_error = " << rel_error << std::endl;
    }
  
  return return_flag;
}

template<typename Thermo , typename Conductivity >

int test_k

(

const libMesh::Real

k

)

template<typename Viscosity >

int test_mu

(

const libMesh::Real

mu

)

template<>

int test_mu< Antioch::MixtureViscosity< Antioch::BlottnerViscosity< libMesh::Real > > >

(

const libMesh::Real

mu

)

Definition at line 74 of file antioch_wilke_evaluator_regression.C.

References test_generic().

{
  double mu_reg = 4.5123309407810213e-05;

  return test_generic(mu,mu_reg,"mu");
}