Hi Dan!
That's not so difficult.
Mesh: you can greate a *.inp file containing the mesh data.
Boundary conditions around the cylinder:
You can use this piece of code from me (for 2D):
GridIn<dim> grid_in;
grid_in.attach_triangulation (triangulation);
std::ifstream input_file("mesh.inp");
Assert (dim==2, ExcInternalError());
grid_in.read_ucd (input_file);
// Center of the cylinder:
Point<dim> p(x ,y);
double radius = radius_of_cylinder;
static const HyperBallBoundary<dim> boundary(p,radius);
// Boundary color of the cylinder is 80
triangulation.set_boundary (80, boundary);
Cheers,
Thomas
--
++----------------------------------------++
Thomas Wick
University of Heidelberg
AG Numerik, Room 213
Im Neuenheimer Feld 293
69120 Heidelberg, Germany
Phone: +49(0)6221 / 54-5449
Email: [email protected]
www: http://numerik.uni-hd.de/~twick/
++----------------------------------------++
--
On Thu, 7 Oct 2010, Keith Brauss wrote:
Hello all,
I am interested in solving a incompressible fluid flow (in a piper or
rectangular duct domain) around a cylinder
contained within the pipe or duct domain. I was wondering if anybody
could point me in the right direction concerning making the mesh and
implementing the boundary conditions on the cylinder. Thanks.
Dan Brauss
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
