Hi,
how can I remesh my domain and update the equation system?
I have tried something like
Mesh mesh (2);
{
std::vector<Hole*> holes;
PolygonHole hole_1(Point(0.1, 0.1), // center
radius_hole, // radius
8); // n. points
MeshTools::Generation::build_delaunay_square(mesh,
3,
3,
0,
size_x,
0,
size_y,
TRI3,
&holes );
}
EquationSystems equation_systems (mesh);
...
...
...
TIME LOOP {
{
//some chemical calculations
// and solving the equation system for diffusion
// Create the vector of Hole*
std::vector<Hole*> holes;
holes.push_back(&hole_1);
// Attach the list of holes to the triangulator object
//t.attach_hole_list(&holes);
mesh.clear();
MeshTools::Generation::build_delaunay_square (mesh,
3,
3,
0,
size_x,
0,
size_y,
TRI3,
&holes );
equation_systems.reinit();
}
}
New holes are created in areas where chemical reation has been completed
and nothing can happen any more - this changes at every time step.
Furthermore, for my simplified model the holes should have circular
shapes
thanks,
Robert
------------------------------------------------------------------------------
Oracle to DB2 Conversion Guide: Learn learn about native support for PL/SQL,
new data types, scalar functions, improved concurrency, built-in packages,
OCI, SQL*Plus, data movement tools, best practices and more.
http://p.sf.net/sfu/oracle-sfdev2dev
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
