antioch_mixture_averaged_transport_evaluator_regression.C File Reference

#include "grins_config.h"
#include <iomanip>
#include <limits>
#include <vector>
#include "grins/antioch_mixture_averaged_transport_mixture.h"
#include "grins/antioch_mixture_averaged_transport_evaluator.h"
#include "grins/cached_values.h"
#include "grins/materials_parsing.h"
#include "grins/physics_naming.h"
#include "grins/antioch_mixture_averaged_transport_mixture_builder.h"
#include "libmesh/getpot.h"

Include dependency graph for antioch_mixture_averaged_transport_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::BlottnerViscosity< libMesh::Real > > (const libMesh::Real mu)
template<typename KineticsThermo , 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 176 of file antioch_mixture_averaged_transport_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::CEACurveFit<libMesh::Real>,
                               Antioch::StatMechThermodynamics<libMesh::Real>,
                               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 KineticsThermo , typename Thermo , typename Viscosity , typename Conductivity , typename Diffusivity >

int test_evaluator

(

const GetPot &

input

)

Definition at line 118 of file antioch_mixture_averaged_transport_evaluator_regression.C.

References GRINS::AntiochMixtureAveragedTransportMixtureBuilder::build_mixture(), GRINS::AntiochEvaluator< Antioch::CEACurveFit< libMesh::Real >, Thermo >::cp(), GRINS::AntiochMixtureAveragedTransportEvaluator< KineticsThermoCurveFit, Thermo, Viscosity, Conductivity, Diffusivity >::mu_and_k_and_D(), and GRINS::AntiochEvaluator< Antioch::CEACurveFit< libMesh::Real >, Thermo >::R_mix().

{
  GRINS::AntiochMixtureAveragedTransportMixtureBuilder builder;

  libMesh::UniquePtr<GRINS::AntiochMixtureAveragedTransportMixture<KineticsThermo,Thermo,Viscosity,Conductivity,Diffusivity> >
    mixture_ptr = builder.build_mixture<KineticsThermo,Thermo,Viscosity,Conductivity,Diffusivity>
    (input,"TestMaterial");


  const GRINS::AntiochMixtureAveragedTransportMixture<KineticsThermo,Thermo,Viscosity,Conductivity,Diffusivity> &
    mixture = *mixture_ptr;

  GRINS::AntiochMixtureAveragedTransportEvaluator<KineticsThermo,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);

  libMesh::Real p0 = rho*T*evaluator.R_mix(Y);
  libMesh::Real mu = 0.0;
  libMesh::Real k = 0.0;
  std::vector<libMesh::Real> D(n_species,0.0);

  evaluator.mu_and_k_and_D( T, rho, evaluator.cp(T,p0,Y), Y, mu, k, 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_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 46 of file antioch_mixture_averaged_transport_evaluator_regression.C.

Referenced by test_mu< 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::BlottnerViscosity< libMesh::Real > >

(

const libMesh::Real

mu

)

Definition at line 77 of file antioch_mixture_averaged_transport_evaluator_regression.C.

References test_generic().

{
  double mu_reg = 4.5123309407810213e-05;

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