Hello,
I'm using getfem to solve Stokes problem in a cube (as in tests/stokes.cc),
but I want to put a sphere inside the cube, as if it was a solid particle
in a volume of fluid, so I don't want to solve stokes inside the sphere.
I'm using level-sets to cut the mesh and the option INTEGRATE_OUTSIDE to
solve the problem only outside the region cut by the level set function.
How can I set no-slip boundary conditions on the sphere, i.e. Dirichlet
condition?
Is there any way to export the cut mesh in vtk? I tried with
getfem::vtk_export vtkp("stokes.vtk");
vtkp.exporting(mfls);
vtkp.write_mesh();
where mfls is a mesh_fem_level_set object, but loading the stokes.vtk file
in paraview I see the whole mesh, and no sign of the cut.
Thank you in advance.
Irene
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
