#include <RadiativeTransportEquationSystem.h>

Inheritance diagram for sierra::nalu::RadiativeTransportEquationSystem:

Public Member Functions
	RadiativeTransportEquationSystem (EquationSystems &equationSystems, const int quadratureOrder, const bool activateScattering, const bool activateUpwind, const bool deactivateSucv, const bool externalCoupling)

virtual	~RadiativeTransportEquationSystem ()

void	register_nodal_fields (stk::mesh::Part *part)

void	register_edge_fields (stk::mesh::Part *part)

void	register_element_fields (stk::mesh::Part *part, const stk::topology &theTopo)

void	register_interior_algorithm (stk::mesh::Part *part)

void	register_wall_bc (stk::mesh::Part *part, const stk::topology &theTopo, const WallBoundaryConditionData &wallBCData)

void	initialize ()

void	predict_state ()

void	solve_and_update ()
	Assemble the LHS and RHS and perform linear solve for prescribed number of iterations. More...

void	set_current_ordinate_info (const int k)

void	initialize_intensity ()

void	compute_bc_intensity ()

void	compute_radiation_source ()

bool	system_is_converged ()

double	provide_scaled_norm ()

double	provide_norm ()

void	zero_out_fields ()

void	zero_irradiation ()

void	assemble_boundary_area ()

void	assemble_fields ()

void	assemble_irradiation ()

void	normalize_irradiation ()

void	compute_div_norm ()

void	copy_ordinate_intensity (const ScalarFieldType &fromField, const ScalarFieldType &toField)

void	get_current_ordinate_info (double &weight, double *Sk) const

void	get_current_ordinate (double *Sk) const

double	get_stefan_boltzmann () const

ScalarFieldType *	get_intensity () const

void	create_quadrature_set ()

Public Member Functions inherited from sierra::nalu::EquationSystem
	EquationSystem (EquationSystems &eqSystems, const std::string name="no_name", const std::string eqnTypeName="no_eqn_type_name")

virtual	~EquationSystem ()

void	set_nodal_gradient (const std::string &dofName)

virtual void	initial_work ()

virtual void	populate_derived_quantities ()

virtual void	pre_iter_work ()
	Perform setup tasks before entering the solve and update step. More...

virtual void	post_iter_work ()
	Perform setup tasks after he solve and update step. More...

virtual void	post_iter_work_dep ()
	Deprecated post iteration work logic. More...

virtual void	assemble_and_solve (stk::mesh::FieldBase *deltaSolution)

virtual void	provide_output ()

virtual void	pre_timestep_work ()

virtual void	reinitialize_linear_system ()

virtual void	post_adapt_work ()

virtual void	dump_eq_time ()

virtual double	provide_norm_increment ()

virtual void	register_inflow_bc (stk::mesh::Part *part, const stk::topology &theTopo, const InflowBoundaryConditionData &inflowBCData)

virtual void	register_open_bc (stk::mesh::Part *part, const stk::topology &theTopo, const OpenBoundaryConditionData &openBCData)

virtual void	register_symmetry_bc (stk::mesh::Part *part, const stk::topology &theTopo, const SymmetryBoundaryConditionData &symmetryBCData)

virtual void	register_periodic_bc (stk::mesh::Part partMaster, stk::mesh::Part partSlave, const stk::topology &theTopoMaster, const stk::topology &theTopoSlave, const PeriodicBoundaryConditionData &periodicBCData)

virtual void	register_non_conformal_bc (stk::mesh::Part *part, const stk::topology &theTopo)

virtual void	register_overset_bc ()

virtual void	create_constraint_algorithm (stk::mesh::FieldBase *theField)

virtual void	register_surface_pp_algorithm (const PostProcessingData &theData, stk::mesh::PartVector &partVector)

virtual void	register_initial_condition_fcn (stk::mesh::Part *part, const std::map< std::string, std::string > &theNames, const std::map< std::string, std::vector< double > > &theParams)

virtual void	evaluate_properties ()

PecletFunction *	create_peclet_function (const std::string dofName)

virtual void	load (const YAML::Node &node)

Simulation *	root ()

EquationSystems *	parent ()

void	report_invalid_supp_alg_names ()

void	report_built_supp_alg_names ()

bool	supp_alg_is_requested (std::string name)

bool	supp_alg_is_requested (std::vector< std::string >)

bool	nodal_src_is_requested ()

void	update_iteration_statistics (const int &iters)

bool	bc_data_specified (const UserData &, std::string &name)

UserDataType	get_bc_data_type (const UserData &, std::string &name)

std::string	get_bc_function_name (const UserData &, std::string &name)

std::vector< double >	get_bc_function_params (const UserData &, std::string &name)

std::vector< std::string >	get_bc_function_string_params (const UserData &, std::string &name)

virtual void	post_converged_work ()

Public Attributes
const int	quadratureOrder_

const bool	activateScattering_

const bool	activateUpwind_

const bool	deactivateSucv_

const bool	externalCoupling_

ScalarFieldType *	intensity_

ScalarFieldType *	currentIntensity_

ScalarFieldType *	intensityBc_

ScalarFieldType *	emissivity_

ScalarFieldType *	transmissivity_

ScalarFieldType *	environmentalT_

ScalarFieldType *	iTmp_

ScalarFieldType *	dualNodalVolume_

VectorFieldType *	coordinates_

ScalarFieldType *	temperature_

VectorFieldType *	radiativeHeatFlux_

ScalarFieldType *	divRadiativeHeatFlux_

ScalarFieldType *	radiationSource_

ScalarFieldType *	scalarFlux_

ScalarFieldType *	scalarFluxOld_

ScalarFieldType *	absorptionCoeff_

ScalarFieldType *	scatteringCoeff_

VectorFieldType *	edgeAreaVec_

ScalarFieldType *	irradiation_

ScalarFieldType *	bcTemperature_

ScalarFieldType *	assembledBoundaryArea_

AlgorithmDriver *	bcIntensityAlgDriver_

bool	isInit_

int	ordinateDirections_

std::vector< double >	Sn_

std::vector< double >	weights_

std::vector< double >	currentSn_

double	currentWeight_

double	stefanBoltz_

double	systemL2Norm_

double	nonLinearResidualSum_

double	firstNonLinearResidualSum_

std::vector< stk::mesh::Part * >	interiorPartVec_

std::vector< stk::mesh::Part * >	bcPartVec_

Public Attributes inherited from sierra::nalu::EquationSystem
EquationSystems &	equationSystems_

Realm &	realm_

std::string	name_

std::string	userSuppliedName_

const std::string	eqnTypeName_

int	maxIterations_

double	convergenceTolerance_

SolverAlgorithmDriver *	solverAlgDriver_

double	timerAssemble_

double	timerLoadComplete_

double	timerSolve_

double	timerMisc_

double	timerInit_

double	timerPrecond_

double	avgLinearIterations_

double	maxLinearIterations_

double	minLinearIterations_

int	nonLinearIterationCount_

bool	reportLinearIterations_

bool	firstTimeStepSolve_

bool	edgeNodalGradient_

std::vector< AuxFunctionAlgorithm * >	bcDataAlg_

std::vector< Algorithm * >	bcDataMapAlg_

std::vector< Algorithm * >	copyStateAlg_

LinearSystem *	linsys_

size_t	num_graph_entries_

std::vector< Algorithm * >	propertyAlg_

std::vector< AlgorithmDriver * >	preIterAlgDriver_
	List of tasks to be performed before each solve_and_update of this equation system. More...

std::vector< AlgorithmDriver * >	postIterAlgDriver_
	List of tasks to be performed after each solve_and_update of this equation system. More...

Constructor & Destructor Documentation

sierra::nalu::RadiativeTransportEquationSystem::RadiativeTransportEquationSystem	(	EquationSystems &	equationSystems,
		const int	quadratureOrder,
		const bool	activateScattering,
		const bool	activateUpwind,
		const bool	deactivateSucv,
		const bool	externalCoupling
	)

sierra::nalu::RadiativeTransportEquationSystem::~RadiativeTransportEquationSystem ( )

virtual

Member Function Documentation

void sierra::nalu::RadiativeTransportEquationSystem::assemble_boundary_area ( )

References assembledBoundaryArea_, bcPartVec_, sierra::nalu::Realm::bulk_data(), sierra::nalu::Realm::get_buckets(), sierra::nalu::Realm::hasPeriodic_, sierra::nalu::Realm::meta_data(), sierra::nalu::Realm::periodic_field_update(), and sierra::nalu::EquationSystem::realm_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::assemble_fields ( )

References currentSn_, currentWeight_, sierra::nalu::Realm::get_buckets(), intensity_, interiorPartVec_, sierra::nalu::Realm::meta_data(), radiativeHeatFlux_, sierra::nalu::EquationSystem::realm_, and scalarFlux_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::assemble_irradiation ( )

References bcPartVec_, currentSn_, currentWeight_, sierra::nalu::Realm::get_buckets(), sierra::nalu::MasterElementRepo::get_surface_master_element(), intensity_, irradiation_, sierra::nalu::Realm::meta_data(), sierra::nalu::MasterElement::nodesPerElement_, sierra::nalu::MasterElement::numIntPoints_, sierra::nalu::EquationSystem::realm_, sierra::nalu::Realm::realmUsesEdges_, sierra::nalu::MasterElement::shape_fcn(), and sierra::nalu::MasterElement::shifted_shape_fcn().

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::compute_bc_intensity ( )

References bcPartVec_, bcTemperature_, emissivity_, environmentalT_, sierra::nalu::Realm::get_buckets(), get_stefan_boltzmann(), intensityBc_, irradiation_, sierra::nalu::Realm::meta_data(), sierra::nalu::EquationSystem::realm_, and transmissivity_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::compute_div_norm ( )

References absorptionCoeff_, divRadiativeHeatFlux_, sierra::nalu::Realm::get_buckets(), get_stefan_boltzmann(), interiorPartVec_, sierra::nalu::Realm::l2Scaling_, sierra::nalu::Realm::meta_data(), sierra::nalu::EquationSystem::realm_, scalarFlux_, scalarFluxOld_, sierra::nalu::NaluEnv::self(), systemL2Norm_, and temperature_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::compute_radiation_source ( )

References absorptionCoeff_, externalCoupling_, sierra::nalu::Realm::get_buckets(), get_stefan_boltzmann(), interiorPartVec_, sierra::nalu::Realm::meta_data(), radiationSource_, sierra::nalu::EquationSystem::realm_, and temperature_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::copy_ordinate_intensity	(	const ScalarFieldType &	fromField,
		const ScalarFieldType &	toField
	)

References sierra::nalu::Realm::bulk_data(), sierra::nalu::field_copy(), sierra::nalu::Realm::get_activate_aura(), sierra::nalu::Realm::meta_data(), and sierra::nalu::EquationSystem::realm_.

Referenced by initialize_intensity(), set_current_ordinate_info(), and solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::create_quadrature_set ( )

References sierra::nalu::Realm::meta_data(), sierra::nalu::NaluEnv::naluOutputP0(), ordinateDirections_, quadratureOrder_, sierra::nalu::EquationSystem::realm_, sierra::nalu::NaluEnv::self(), Sn_, and weights_.

Referenced by RadiativeTransportEquationSystem().

void sierra::nalu::RadiativeTransportEquationSystem::get_current_ordinate ( double * Sk ) const

References currentSn_, sierra::nalu::Realm::meta_data(), and sierra::nalu::EquationSystem::realm_.

Referenced by sierra::nalu::AssembleRadTransEdgeSolverAlgorithm::execute(), sierra::nalu::AssembleRadTransEdgeUpwindSolverAlgorithm::execute(), sierra::nalu::AssembleRadTransElemSolverAlgorithm::execute(), sierra::nalu::AssembleRadTransWallSolverAlgorithm::execute(), sierra::nalu::RadTransFemElemSuppAlg::setup(), and sierra::nalu::RadTransSupgElemSuppAlg::setup().

void sierra::nalu::RadiativeTransportEquationSystem::get_current_ordinate_info	(	double &	weight,
		double *	Sk
	)		const

References currentSn_, currentWeight_, sierra::nalu::Realm::meta_data(), and sierra::nalu::EquationSystem::realm_.

ScalarFieldType * sierra::nalu::RadiativeTransportEquationSystem::get_intensity ( ) const

References intensity_.

Referenced by sierra::nalu::AssembleRadTransEdgeSolverAlgorithm::execute(), sierra::nalu::AssembleRadTransEdgeUpwindSolverAlgorithm::execute(), sierra::nalu::AssembleRadTransElemSolverAlgorithm::execute(), sierra::nalu::AssembleRadTransWallSolverAlgorithm::execute(), sierra::nalu::RadTransBlackBodyNodeSuppAlg::setup(), sierra::nalu::RadTransFemElemSuppAlg::setup(), and sierra::nalu::RadTransSupgElemSuppAlg::setup().

double sierra::nalu::RadiativeTransportEquationSystem::get_stefan_boltzmann ( ) const

References sierra::nalu::Realm::get_stefan_boltzmann(), and sierra::nalu::EquationSystem::realm_.

Referenced by compute_bc_intensity(), compute_div_norm(), compute_radiation_source(), and initialize_intensity().

void sierra::nalu::RadiativeTransportEquationSystem::initialize ( )

virtual

Reimplemented from sierra::nalu::EquationSystem.

References sierra::nalu::LinearSystem::finalizeLinearSystem(), sierra::nalu::SolverAlgorithmDriver::initialize_connectivity(), sierra::nalu::EquationSystem::linsys_, and sierra::nalu::EquationSystem::solverAlgDriver_.

void sierra::nalu::RadiativeTransportEquationSystem::initialize_intensity ( )

References copy_ordinate_intensity(), sierra::nalu::Realm::get_buckets(), get_stefan_boltzmann(), intensity_, interiorPartVec_, sierra::nalu::Realm::meta_data(), ordinateDirections_, sierra::nalu::EquationSystem::realm_, and temperature_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::normalize_irradiation ( )

References assembledBoundaryArea_, bcPartVec_, sierra::nalu::Realm::bulk_data(), sierra::nalu::Realm::get_buckets(), sierra::nalu::Realm::hasPeriodic_, irradiation_, sierra::nalu::Realm::meta_data(), sierra::nalu::Realm::periodic_field_update(), and sierra::nalu::EquationSystem::realm_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::predict_state ( )

virtual

Reimplemented from sierra::nalu::EquationSystem.

double sierra::nalu::RadiativeTransportEquationSystem::provide_norm ( )

virtual

Reimplemented from sierra::nalu::EquationSystem.

References nonLinearResidualSum_.

double sierra::nalu::RadiativeTransportEquationSystem::provide_scaled_norm ( )

virtual

Reimplemented from sierra::nalu::EquationSystem.

References firstNonLinearResidualSum_, and nonLinearResidualSum_.

void sierra::nalu::RadiativeTransportEquationSystem::register_edge_fields ( stk::mesh::Part * part )

virtual

Reimplemented from sierra::nalu::EquationSystem.

References edgeAreaVec_, sierra::nalu::Realm::meta_data(), sierra::nalu::EquationSystem::realm_, and sierra::nalu::Realm::realmUsesEdges_.

void sierra::nalu::RadiativeTransportEquationSystem::register_element_fields	(	stk::mesh::Part *	part,
		const stk::topology &	theTopo
	)

virtual

Reimplemented from sierra::nalu::EquationSystem.

void sierra::nalu::RadiativeTransportEquationSystem::register_interior_algorithm ( stk::mesh::Part * part )

virtual

Reimplemented from sierra::nalu::EquationSystem.

References activateScattering_, activateUpwind_, sierra::nalu::SolutionOptions::elemSrcTermsMap_, sierra::nalu::INTERIOR, interiorPartVec_, sierra::nalu::EquationSystem::realm_, sierra::nalu::Realm::realmUsesEdges_, sierra::nalu::Realm::solutionOptions_, sierra::nalu::EquationSystem::solverAlgDriver_, sierra::nalu::SolverAlgorithmDriver::solverAlgMap_, sierra::nalu::SRC, sierra::nalu::Algorithm::supplementalAlg_, and sierra::nalu::SolutionOptions::useConsolidatedSolverAlg_.

void sierra::nalu::RadiativeTransportEquationSystem::register_nodal_fields ( stk::mesh::Part * part )

virtual

void sierra::nalu::RadiativeTransportEquationSystem::register_wall_bc	(	stk::mesh::Part *	part,
		const stk::topology &	theTopo,
		const WallBoundaryConditionData &	wallBCData
	)

virtual

Reimplemented from sierra::nalu::EquationSystem.

References assembledBoundaryArea_, sierra::nalu::EquationSystem::bcDataAlg_, bcPartVec_, bcTemperature_, sierra::nalu::Emissivity::emissivity_, emissivity_, sierra::nalu::WallUserData::emissivity_, sierra::nalu::EnvironmentalT::environmentalT_, environmentalT_, sierra::nalu::WallUserData::environmentalT_, sierra::nalu::Realm::initCondAlg_, intensityBc_, irradiation_, sierra::nalu::WallUserData::isInterface_, sierra::nalu::Realm::meta_data(), sierra::nalu::EquationSystem::realm_, sierra::nalu::Realm::realmUsesEdges_, sierra::nalu::EquationSystem::solverAlgDriver_, sierra::nalu::SolverAlgorithmDriver::solverAlgMap_, sierra::nalu::Temperature::temperature_, sierra::nalu::UserData::temperature_, sierra::nalu::UserData::tempSpec_, sierra::nalu::Transmissivity::transmissivity_, transmissivity_, sierra::nalu::WallUserData::transmissivity_, sierra::nalu::WallBoundaryConditionData::userData_, and sierra::nalu::WALL.

void sierra::nalu::RadiativeTransportEquationSystem::set_current_ordinate_info ( const int k )

References copy_ordinate_intensity(), currentIntensity_, currentSn_, currentWeight_, intensity_, sierra::nalu::Realm::meta_data(), sierra::nalu::EquationSystem::realm_, Sn_, and weights_.

Referenced by solve_and_update().

void sierra::nalu::RadiativeTransportEquationSystem::solve_and_update ( )

virtual

Assemble the LHS and RHS and perform linear solve for prescribed number of iterations.

This method is invoked in EquationSystems::solve_and_update method as shown below

pre_iter_work();
// Iterate over all equation systems
for (auto eqsys: equationSystems_) {
  eqsys->pre_iter_work();
  eqsys->solve_and_update();             //<<<< Assemble and solve system
  eqsys->post_iter_work();
}
post_iter_work();