uns_adapt_3d_fs.cc File Reference

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

Classes
class	oomph::RefineableTetgenMesh< ELEMENT >
class	oomph::RefineableSolidTetgenMesh< ELEMENT >
class	RisingBubbleProblem< ELEMENT >
	Micky mouse problem. More...

Namespaces
	Global_Parameters
	Namespace for global parameters.
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).

Functions
Vector< double >	Global_Parameters::G (3)
	Gravity. More...
int	main (int argc, char *argv[])
	Demonstrate how to solve Poisson problem. More...

Variables
double	Global_Parameters::Ca = 1.0
	Set the Capillary number. More...
double	Global_Parameters::Volume = -16.0*atan(1.0)*Radius*Radius*Radius/3.0
double	Global_Parameters::ReInvFr = 0.0
	The Reynolds INverse Froude number. More...

Function Documentation

main()

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

Demonstrate how to solve Poisson problem.

{
 
 //Set the direction of gravity
 Global_Parameters::G[0] = 0.0;
 Global_Parameters::G[1] = 0.0;
 Global_Parameters::G[2] = -1.0;
 
 // Label for output
 DocInfo doc_info;
 
 // Output directory
 doc_info.set_directory("RESLT");
  
 // Number of output points per edge
 unsigned nplot=2;
 
 // Do the problem with quadratic elements
 //---------------------------------------
 {
  RisingBubbleProblem<ProjectableTaylorHoodElement<
  PseudoSolidNodeUpdateElement<TTaylorHoodElement<3>,
   TPVDElement<3,3> > > > problem;
  problem.snap_onto_sphere();
  problem.Fluid_mesh_pt->set_lagrangian_nodal_coordinates();
  problem.doc_solution(5,doc_info);
  doc_info.number()++;
 
  Global_Parameters::Volume -= 0.1255;
 
  problem.steady_newton_solve();
  problem.doc_solution(5,doc_info);
  doc_info.number()++;
 
  //Assume it's always been this way
  double dt = 0.1;
  problem.Fluid_mesh_pt->set_lagrangian_nodal_coordinates();
  problem.assign_initial_values_impulsive(dt);
  //Turn on gravity
  Global_Parameters::ReInvFr = 1.0;
 
  for(unsigned t=0;t<10;t++)
   {
    problem.unsteady_newton_solve(dt);
    problem.doc_solution(5,doc_info);
    doc_info.number()++;
    //Reset the lagrangian
    problem.Fluid_mesh_pt->set_lagrangian_nodal_coordinates();
   }
 
  //Now can I adapt
  problem.adapt();
  
  problem.unsteady_newton_solve(dt);
  problem.doc_solution(5,doc_info);
  doc_info.number()++;
  //Reset the lagrangian
  problem.Fluid_mesh_pt->set_lagrangian_nodal_coordinates();
 
  exit(1);
 
  for(unsigned n=0;n<1;++n)
   {
    // Solve the problem
    problem.steady_newton_solve(1);
    
    //Output solution with 5 points per edge
    nplot=5;
    problem.doc_solution(nplot,doc_info);
    
    //Increment counter for solutions 
    doc_info.number()++;
    
    Global_Parameters::Re += 0.1;
   }
 
 }
 
 
}

References Global_Parameters::G, n, oomph::DocInfo::number(), problem, Global_Parameters::Re, Global_Parameters::ReInvFr, oomph::DocInfo::set_directory(), plotPSD::t, and Global_Parameters::Volume.