adaptive_bubble_in_channel.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "solid.h"
#include "constitutive.h"
#include "fluid_interface.h"
#include "meshes/triangle_mesh.h"

Classes
class	oomph::MyTaylorHoodElement< DIM >
	Overload TaylorHood element to modify output. More...
class	oomph::FaceGeometry< MyTaylorHoodElement >
class	oomph::FaceGeometry< FaceGeometry< MyTaylorHoodElement > >
class	BubbleInChannelProblem< ELEMENT >
	Problem class to simulate inviscid bubble propagating along 2D channel. More...

Namespaces
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).
	Problem_Parameter
	Namespace for Problem Parameter.

Functions
int	main (int argc, char **argv)
	Driver code for moving bubble problem. More...

Variables
double	Problem_Parameter::Radius = 0.25
	Initial radius of bubble. More...
double	Problem_Parameter::Volume = -MathematicalConstants::Pi*Radius*Radius
double	Problem_Parameter::Inflow_veloc_magnitude = 0.0

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver code for moving bubble problem.

{
  
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Define possible command line arguments and parse the ones that
 // were actually specified
 
 // Validation?
 CommandLineArgs::specify_command_line_flag("--validation");
 
 // Parse command line
 CommandLineArgs::parse_and_assign(); 
 
 // Doc what has actually been specified on the command line
 CommandLineArgs::doc_specified_flags();
 
 // Create generalised Hookean constitutive equations
 Problem_Parameter::Constitutive_law_pt = 
  new GeneralisedHookean(&Problem_Parameter::Nu);
 
 // Open trace file
 Problem_Parameter::Trace_file.open("RESLT/trace.dat"); 
 
 // Increase precision of output
 Problem_Parameter::Trace_file.precision(20);
 
 // Open norm file
 Problem_Parameter::Norm_file.open("RESLT/norm.dat");
 
 
 // Create problem in initial configuration
 BubbleInChannelProblem<ProjectableTaylorHoodElement<MyTaylorHoodElement> > 
  problem;  
 
 
 // Output the problem's state with the bubble in its
 // initial polygonal representation
 problem.doc_solution();
 
 // Before starting the time-integration we want to "inflate" it to form 
 // a proper circular bubble. We do this by setting the inflow to zero
 // and doing a steady solve (with one adaptation)
 Problem_Parameter::Inflow_veloc_magnitude=0.0;
 
 problem.steady_newton_solve(1);
 
 // If all went well, this should show us a nice circular bubble
 // in a stationary fluid
 problem.doc_solution();
 
 // Initialise timestepper
 double dt=0.025;
 problem.initialise_dt(dt);
 
 // Perform impulsive start from current state
 problem.assign_initial_values_impulsive();
 
 
 // Now switch on the inflow and re-assign the boundary conditions
 // (Call to complete_problem_setup() is a bit expensive given that we
 // we only want to set the inflow velocity but who cares -- it's just
 // a one off.
 Problem_Parameter::Inflow_veloc_magnitude=1.0;
 problem.complete_problem_setup();
 
 
 // Solve problem on fixed mesh
 unsigned nstep=6;
 if (CommandLineArgs::command_line_flag_has_been_set("--validation"))
  {
   nstep=2;
   oomph_info << "Remeshing after every second step during validation\n";
  }
 for (unsigned i=0;i<nstep;i++)
  {
   // Solve the problem
   problem.unsteady_newton_solve(dt);    
   problem.doc_solution();
  } // done solution on fixed mesh
 
 
 // Now do a proper loop, doing nstep timesteps before adapting/remeshing
 // and repeating the lot ncycle times
 unsigned ncycle=1000;
 if (CommandLineArgs::command_line_flag_has_been_set("--validation"))
  {
   ncycle=1;
   oomph_info << "Only doing one cycle during validation\n";
  }
 
 // Do the cycles
 for(unsigned j=0;j<ncycle;j++)
  {       
   // Allow up to one level of refinement for next solve
   unsigned max_adapt=1;
 
   //Solve problem a few times
   for (unsigned i=0;i<nstep;i++)
    {     
     // Solve the problem
     problem.unsteady_newton_solve(dt,max_adapt,false); 
 
 
     // Build the label for doc and output solution
     std::stringstream label;
     label << "Adaptation " <<j << " Step "<< i;
     problem.doc_solution(label.str());
 
     // No more refinement for the next nstep steps
     max_adapt=0;
    }
 
  }
 
 
} //End of main

References oomph::CommandLineArgs::command_line_flag_has_been_set(), Problem_Parameter::Constitutive_law_pt, oomph::CommandLineArgs::doc_specified_flags(), i, Problem_Parameter::Inflow_veloc_magnitude, j, UniformPSDSelfTest::label, Problem_Parameter::Norm_file, Problem_Parameter::Nu, oomph::oomph_info, oomph::CommandLineArgs::parse_and_assign(), problem, Flag_definition::setup(), oomph::CommandLineArgs::specify_command_line_flag(), and Problem_Parameter::Trace_file.