Commit 7f064794 authored by gfourestey's avatar gfourestey
Browse files

flux imposition inserted

parent 18be5aee
...@@ -330,8 +330,8 @@ struct Cylinder::Private ...@@ -330,8 +330,8 @@ struct Cylinder::Private
const Real& /*z*/, const Real& /*z*/,
const ID& id ) const const ID& id ) const
{ {
if (id == 3) // if (id == 3)
return 1.; return 10.;
return 0.; return 0.;
} }
...@@ -424,13 +424,14 @@ Cylinder::run() ...@@ -424,13 +424,14 @@ Cylinder::run()
bcH.addBC( "Cylinder", CYLINDER, Essential, Full, uZero, 3 ); bcH.addBC( "Cylinder", CYLINDER, Essential, Full, uZero, 3 );
// bcH.addBC( "Slipwall", SLIPWALL, Essential, Full, uZero , 3 ); // bcH.addBC( "Slipwall", SLIPWALL, Essential, Full, uZero , 3 );
#endif #endif
#ifdef TUBE20_MESH_SETTINGS #ifdef TUBE20_MESH_SETTINGS
//BCFunctionBase unormal( d->get_normal() ); //BCFunctionBase unormal( d->get_normal() );
//cylinder //cylinder
// bcH.addBC( "Inlet", INLET, Essential, Full, uPois, 3 ); //bcH.addBC( "Inlet", INLET, Natural, Full, uIn, 3 );
bcH.addBC( "Inlet", INLET, Essential, Full, uPois, 3 ); bcH.addBC( "Inlet", INLET, Flux, Full, uOne, 3);
bcH.addBC( "Outlet", OUTLET, Natural, Full, uZero, 3 ); bcH.addBC( "Outlet", OUTLET, Natural, Full, uZero, 3 );
//bcH.addBC( "Wall", WALL, Natural, Full, uNormal, 3 ); //bcH.addBC( "Wall", WALL, Natural, Full, uNormal, 3 );
//bcH.addBC( "Wall", WALL, Natural, Full, uNormal, 3 ); //bcH.addBC( "Wall", WALL, Natural, Full, uNormal, 3 );
...@@ -439,6 +440,8 @@ Cylinder::run() ...@@ -439,6 +440,8 @@ Cylinder::run()
bcH.addBC( "RingOut", RINGOUT, Essential, Full, uZero, 3 ); bcH.addBC( "RingOut", RINGOUT, Essential, Full, uZero, 3 );
#endif #endif
int numLM = 1;
DataNavierStokes<RegionMesh3D<LinearTetra> > dataNavierStokes( dataFile ); DataNavierStokes<RegionMesh3D<LinearTetra> > dataNavierStokes( dataFile );
partitionMesh< RegionMesh3D<LinearTetra> > meshPart(*dataNavierStokes.mesh(), *d->comm); partitionMesh< RegionMesh3D<LinearTetra> > meshPart(*dataNavierStokes.mesh(), *d->comm);
...@@ -472,15 +475,19 @@ Cylinder::run() ...@@ -472,15 +475,19 @@ Cylinder::run()
UInt totalPressDof = pFESpace.map().getMap(Unique)->NumGlobalElements(); UInt totalPressDof = pFESpace.map().getMap(Unique)->NumGlobalElements();
if (verbose) std::cout << "Total Velocity Dof = " << totalVelDof << std::endl; if (verbose) std::cout << "Total Velocity Dof = " << totalVelDof << std::endl;
if (verbose) std::cout << "Total Pressure Dof = " << totalPressDof << std::endl; if (verbose) std::cout << "Total Pressure Dof = " << totalPressDof << std::endl;
if (verbose) std::cout << "Calling the fluid constructor ... "; if (verbose) std::cout << "Calling the fluid constructor ... ";
bcH.setOffset("Inlet", totalVelDof + totalPressDof);
Oseen< RegionMesh3D<LinearTetra> > fluid (dataNavierStokes, Oseen< RegionMesh3D<LinearTetra> > fluid (dataNavierStokes,
uFESpace, uFESpace,
pFESpace, pFESpace,
*d->comm); *d->comm,
numLM);
EpetraMap fullMap(fluid.getMap()); EpetraMap fullMap(fluid.getMap());
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment