two_d_unsteady_heat_2adapt.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);
 
 
 // Specify duration of the simulation
 double t_max=0.04; // 2.0;
 
 // Build problem: Pass pointer to source function and initial timestep
 RefineableUnsteadyHeatProblem<RefineableQUnsteadyHeatElement<2,3> >
  problem(&TanhSolnForUnsteadyHeat::get_source);
   
 // 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 or the "suggested next dt" from the restarted run
 double dt=problem.next_dt();
 
 cout << "Doing first timestep with dt: " << dt << std::endl;
 
 // 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();
  }
 
 // Write header for trace file
 problem.write_trace_file_header();
 
 
 // Timestepping loop
 while (problem.time_pt()->time()<t_max)
  {
   // Take timestep with temporal and spatial adaptivity
   double dt_new=
   problem.doubly_adaptive_unsteady_newton_solve(dt,problem.epsilon_t(),
                                                 max_adapt,first);
   cout << "Suggested new dt: " << dt_new << std::endl;
   dt=dt_new;
   
   // Store for restart
   problem.next_dt()=dt_new;
 
   // 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 -- too scared that the interpolation error will 
   // wipe out any further gains...
   max_adapt=1;
   
   //Output solution
   problem.doc_solution();
 
  }
 
 
}; // end of main

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