two_d_unsteady_heat_ALE.cc File Reference

#include "generic.h"
#include "unsteady_heat.h"
#include "meshes/quarter_circle_sector_mesh.h"

Classes
class	MyEllipse
class	RefineableUnsteadyHeatProblem< ELEMENT >
	Unsteady heat problem in deformable ellipse domain. More...

Namespaces
	TanhSolnForUnsteadyHeat

Functions
double	TanhSolnForUnsteadyHeat::step_position (const double &time)
	Position of step (x-axis intercept) More...
void	TanhSolnForUnsteadyHeat::get_exact_u (const double &time, const Vector< double > &x, Vector< double > &u)
	Exact solution as a Vector. More...
void	TanhSolnForUnsteadyHeat::get_exact_u (const double &time, const Vector< double > &x, double &u)
	Exact solution as a scalar. More...
void	TanhSolnForUnsteadyHeat::get_source (const double &time, const Vector< double > &x, double &source)
	Source function to make it an exact solution. More...
void	TanhSolnForUnsteadyHeat::prescribed_flux_on_fixed_y_boundary (const double &time, const Vector< double > &x, double &flux)
	Flux required by the exact solution on a boundary on which y is fixed. More...
int	main (int argc, char *argv[])

Function Documentation

main()

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

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// Demonstrate how to solve an unsteady heat problem in deformable domain with mesh adaptation. Command line arguments specify the name of the restart file.

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Build problem
 RefineableUnsteadyHeatProblem<RefineableQUnsteadyHeatElement<2,3> >
  problem(&TanhSolnForUnsteadyHeat::get_source);
   
 // Specify duration of the simulation
// double t_max=3.0;
 
 // Set targets for spatial adaptivity
 problem.bulk_mesh_pt()->max_permitted_error()=0.001;
 problem.bulk_mesh_pt()->min_permitted_error()=0.0001;
 
 // First timestep?
 bool first=true;
 
 // Max. number of spatial adaptations per timestep. Allow plenty
 // of adaptations at first timestep as the initial conditions
 // can be reset "exactly" from without any interpolation error.
 unsigned max_adapt=10; 
 
 // Set IC
 problem.set_initial_condition();
 
 // Initial timestep: Use the one used when setting up the initial
 // condition
 double dt=problem.time_pt()->dt();
 
 // If restart: The first step isn't really the first step,
 // i.e. initial condition should not be re-set when 
 // adaptive refinement has been performed. Also, limit
 // the max. number of refinements per timestep to the
 // normal value straightaway.
 if (CommandLineArgs::Argc==2)
  {
   first=false;
   max_adapt=1;
  }
 // If no restart, refine mesh uniformly before we get started
 else
  {
   problem.refine_uniformly();
   problem.refine_uniformly();
  }
 
 //Output solution
 problem.doc_solution();
 
 // find number of steps
 unsigned nstep = 6; //unsigned(t_max/dt);
 
 // Timestepping loop
 for (unsigned istep=0;istep<nstep;istep++)
  {
   // Take timestep 
   problem.unsteady_newton_solve(dt,max_adapt,first);
   
   // Now we've done the first timestep -- don't re-set the IC
   // in subsequent steps
   first=false;
   
   // Reduce the number of spatial adaptations to one per 
   // timestep
   max_adapt=1;
   
   //Output solution
   problem.doc_solution();
 
  }
 
 
}; // end of main

References oomph::CommandLineArgs::Argc, TanhSolnForUnsteadyHeat::get_source(), problem, and Flag_definition::setup().