algebraic_free_boundary_poisson.cc File Reference

#include "generic.h"
#include "poisson.h"
#include "meshes/fish_mesh.h"
#include "circle_as_generalised_element.h"

Classes
class	RefineableFishPoissonProblem< ELEMENT >

Namespaces
	ConstSourceForPoisson
	Namespace for const source term in Poisson equation.

Functions
void	ConstSourceForPoisson::get_source (const Vector< double > &x, double &source)
	Const source function. More...
template<class ELEMENT >
void	demo_fish_poisson (const string &directory_name)
template<class ELEMENT >
void	demo_elastic_fish_poisson (const string &directory_name)
int	main (int argc, char *argv[])

Variables
double	ConstSourceForPoisson::Strength =1.0
	Strength of source function: default value 1.0. More...

Function Documentation

demo_elastic_fish_poisson()

template<class ELEMENT >

void demo_elastic_fish_poisson

(

const string &

directory_name

)

Demonstrate how to solve "elastic" 2D Poisson problem in deformable fish-shaped domain with mesh adaptation.

{
 
 // Loop over all cases
 for (unsigned i_case=0;i_case<5;i_case++)
  //unsigned i_case=1;
  {
   std::cout << "[CR residuals] " << std::endl;
   std::cout << "[CR residuals]==================================================" 
             << std::endl;
   std::cout << "[CR residuals] " << std::endl;
   //Set up the problem with "elastic" fish back
   bool fix_position=false;
   RefineableFishPoissonProblem<ELEMENT> problem(fix_position,
                                                 directory_name,
                                                 i_case);
   
   // Doc refinement targets
   problem.fish_mesh_pt()->doc_adaptivity_targets(cout);
   
   // Do some uniform mesh refinement first
   //--------------------------------------
   problem.refine_uniformly();
   problem.refine_uniformly();
   
   
   // Initial value for load on fish back
   problem.load()=0.0;
   
   // Solve/doc
   unsigned max_solve=2; 
   problem.newton_solve(max_solve);
   problem.doc_solution();
  }
 
 
}

References problem.

demo_fish_poisson()

template<class ELEMENT >

void demo_fish_poisson

(

const string &

directory_name

)

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// Demonstrate how to solve 2D Poisson problem in deformable fish-shaped domain with mesh adaptation.

{
 
 // Set up the problem with prescribed displacement of fish back
 bool fix_position=true;
 RefineableFishPoissonProblem<ELEMENT> problem(fix_position,directory_name,0);
 
 // Doc refinement targets
 problem.fish_mesh_pt()->doc_adaptivity_targets(cout);
  
 // Do some uniform mesh refinement first
 //--------------------------------------
 problem.refine_uniformly();
 problem.refine_uniformly(); 
 
 // Initial value for the vertical displacement of the fish's back
 problem.y_c()=-0.3;
 
 // Loop for different fish shapes
 //-------------------------------
 
 // Number of steps
 unsigned nstep=5;
 
 // Increment in displacement
 double dyc=0.6/double(nstep-1);
 
 // Valiation: Just do one step
 if (CommandLineArgs::Argc>1) nstep=1;
 
 for (unsigned istep=0;istep<nstep;istep++)
  {    
   // Solve/doc
   unsigned max_solve=2; 
   problem.newton_solve(max_solve);
   problem.doc_solution();
   
   //Increment counter for solutions 
   problem.doc_info().number()++;    
   
   // Change vertical displacement
   problem.y_c()+=dyc;
  } 
}

References oomph::CommandLineArgs::Argc, and problem.

main()

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

Driver for "elastic" fish poisson solver with adaptation. If there are any command line arguments, we regard this as a validation run and perform only a single step.

{
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Shorthand for element type
  typedef AlgebraicElement<RefineableQPoissonElement<2,3> > ELEMENT;
 
  // Compute solution of Poisson equation in various domains
  demo_fish_poisson<ELEMENT>("RESLT");
 
  // Compute "elastic" coupled solution directly
  demo_elastic_fish_poisson<ELEMENT>("RESLT_coupled"); 
 
}

References Flag_definition::setup().