poisson/elastic_poisson/elastic_poisson.cc File Reference

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

Classes
class	ElasticFishBackElement
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[])

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.

{
 
  //Set up the problem with "elastic" fish back
 bool fix_position=false;
 RefineableFishPoissonProblem<ELEMENT> problem(fix_position,directory_name);
  
 // 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);  
 
 // 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 load on the fish back 
 problem.load()=-0.3;
 
 // Loop for different fish shapes
 //-------------------------------
 
 // Number of steps
 unsigned nstep=5;
 if (CommandLineArgs::Argc>1) nstep=1;
 
 // Increment in load
 double dp=0.6/double(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 load on ring
   problem.load()+=dp;
  } 
}

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_elastic");
 
 
}

References Flag_definition::setup().