thermo.cc File Reference

#include "generic.h"
#include "solid.h"
#include "unsteady_heat.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	QThermalPVDElement< DIM >
class	ThermalProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char **argv)
	Driver code for 2D Thermoelasticity problem. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver code for 2D Thermoelasticity problem.

{
 
 // "Big G" Linear constitutive equations:
 Global_Physical_Variables::Constitutive_law_pt = 
  new GeneralisedHookean(&Global_Physical_Variables::Nu,
                         &Global_Physical_Variables::E);
 
 //Construct our problem
 ThermalProblem<QThermalPVDElement<2> > problem;
 
 //Number of quasi-steady steps
 unsigned n_steps = 11;
 //If we have additional command line arguemnts, take fewer steps
 if(argc > 1) {n_steps = 2;}
 for(unsigned i=0;i<n_steps;i++)
  {
   //Increase the thermal expansion coefficient
   Global_Physical_Variables::Alpha = 0.1*i;
   
   //Perform a single steady newton solve
   problem.newton_solve();
   //Document the solution
   problem.doc_solution();
  }
} // end of main

References Global_Physical_Variables::Alpha, Global_Physical_Variables::Constitutive_law_pt, Global_Physical_Variables::E, i, oomph::Global_Physical_Variables::Nu, and problem.