two_d_adv_diff_adapt2.cc File Reference

#include "generic.h"
#include "advection_diffusion.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	RefineableAdvectionDiffusionProblem< ELEMENT >

Namespaces
	TanhSolnForAdvectionDiffusion

Functions
void	TanhSolnForAdvectionDiffusion::tanh_profile (const Vector< double > &x, Vector< double > &u)
	Tanh profile for assignment of boundary conditons as a Vector. More...
void	TanhSolnForAdvectionDiffusion::tanh_profile (const Vector< double > &x, double &u)
	Tanh profile for assignment of boundary conditons as a Vector. More...
void	TanhSolnForAdvectionDiffusion::source_function (const Vector< double > &x_vect, double &source)
	Source function required to make the solution above an exact solution. More...
void	TanhSolnForAdvectionDiffusion::wind_function (const Vector< double > &x, Vector< double > &wind)
	Wind. More...
int	main ()
	Driver code for 2D AdvectionDiffusion problem. More...

Function Documentation

main()

int main

(

)

Driver code for 2D AdvectionDiffusion problem.

{
 
 //Set up the problem
 //------------------
 
 // Create the problem with 2D nine-node refineable elements from the
 // RefineableQuadAdvectionDiffusionElement family. Pass pointer to 
 // source and wind function. 
 RefineableAdvectionDiffusionProblem<RefineableQAdvectionDiffusionElement<2,
  3> > problem(&TanhSolnForAdvectionDiffusion::source_function,
               &TanhSolnForAdvectionDiffusion::wind_function);
 
 // Create label for output
 //------------------------
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 // Step number
 doc_info.number()=0;
 
 // Check if we're ready to go:
 //----------------------------
 cout << "\n\n\nProblem self-test ";
 if (problem.self_test()==0) 
  {
   cout << "passed: Problem can be solved." << std::endl;
  }
 else 
  {
   throw OomphLibError("Self test failed",
                       OOMPH_CURRENT_FUNCTION,
                       OOMPH_EXCEPTION_LOCATION);
  }
 
 // Set the orientation of the "step" to 45 degrees
 TanhSolnForAdvectionDiffusion::TanPhi=1.0;
 
 // Choose a large value for the steepness of the "step"
 TanhSolnForAdvectionDiffusion::Alpha=50.0; 
 
 // Doc (default) refinement targets
 //----------------------------------
 problem.mesh_pt()->doc_adaptivity_targets(cout);
 
 
 
 // Solve/doc the problem on the initial, very coarse mesh
 //-------------------------------------------------------
 
 // Solve the problem
 problem.newton_solve();
 
 //Output solution
 problem.doc_solution(doc_info);
 
 //Increment counter for solutions 
 doc_info.number()++; // end of Solve/doc very coarse mesh
 
 // Now do (up to) four rounds of fully automatic adapation in response to 
 //-----------------------------------------------------------------------
 // error estimate
 //---------------
 unsigned max_solve=4;
 for (unsigned isolve=0;isolve<max_solve;isolve++)
  {
   // Adapt problem/mesh
   problem.adapt(); 
         
   // Re-solve the problem if the adaptation has changed anything
   if ((problem.mesh_pt()->nrefined()  !=0)||
       (problem.mesh_pt()->nunrefined()!=0))
    {
     problem.newton_solve();
    }
   else
    {
     cout << "Mesh wasn't adapted --> we'll stop here" << std::endl;
     break;
    }
   
   //Output solution
   problem.doc_solution(doc_info);
   
   //Increment counter for solutions 
   doc_info.number()++;
  } // end (up to) four rounds of fully automatic adapation
 
 
} //end of main

References TanhSolnForAdvectionDiffusion::Alpha, oomph::DocInfo::number(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, problem, oomph::DocInfo::set_directory(), TanhSolnForAdvectionDiffusion::source_function(), TanhSolnForAdvectionDiffusion::TanPhi, and TanhSolnForAdvectionDiffusion::wind_function().