uns_adapt_3d.cc File Reference

#include "generic.h"
#include "navier_stokes.h"
#include "oomph_tetgen.h"
#include "meshes/tetgen_mesh.h"
#include "meshes/refineable_tetgen_mesh.h"

Classes
class	SphericalTetMeshFacetedSurface
	TetMeshFacetedSurface that defines inner boundary. More...
class	CubicTetMeshFacetedSurface
	TetMeshFacetedSurface that defines outer boundary. More...
class	FallingBlockProblem< ELEMENT >
	Micky mouse problem. More...

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
int	main (int argc, char *argv[])
	Demonstrate how to solve Poisson problem. More...

Variables
double	Global_Parameters::Visc_Ratio = 2.0
double	Global_Parameters::Box_width = 3.0
double	Global_Parameters::Box_length = 10.0

Function Documentation

main()

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

Demonstrate how to solve Poisson problem.

{
 // Label for output
 DocInfo doc_info;
 
 // Output directory
 doc_info.set_directory("RESLT");
  
 // Number of output points per edge
 unsigned nplot=2;
 
 //Output trace file
 std::ofstream dissipation("RESLT/diss.dat");
 
 // Do the problem with quadratic elements
 //---------------------------------------
 {
  FallingBlockProblem<ProjectableCrouzeixRaviartElement<
  TCrouzeixRaviartElement<3> > > problem;
  //problem.snap_onto_sphere();
 
  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()++;
    //Output the dissipation
    dissipation << " " << Global_Parameters::Re << "  "
                << problem.get_dissipation() << std::endl;
 
    //Increase the Reynolds number
    Global_Parameters::Re += 0.1;
   }
 
  dissipation.close();
 }
 
 
}

References n, oomph::DocInfo::number(), problem, Global_Parameters::Re, and oomph::DocInfo::set_directory().