Example X 1 : 3D particle sedimentation with sharp interface

In this example we study the sedimentation of a particle in a fluid. This case is based on the case E4 of the experiment realized by Cate et al (click here for a link to the article describing this case). This case model the sedimentation of a particle (D=1.5 cm) in a container of 10x10x16 cm. The fluid choice and the relative density of the particle and the fluid lead to approximated maximal Reynold number of RE=31.9.

We first defined the container using a subdivided hyper rectangle from the deal.ii mesh generation tool.

#---------------------------------------------------
# Mesh
#---------------------------------------------------
subsection mesh
    set type                 = dealii
    set grid type            = subdivided_hyper_rectangle
    set grid arguments       = 5,8,5: 0,0,0 : 10 , 16 ,10 : true
    set initial refinement   = 1
end

We applied no slip boundary condition on all the side except at the top where we let the boundary open.

# --------------------------------------------------
# Boundary Conditions
#---------------------------------------------------
subsection boundary conditions
set number                  = 5
    subsection bc 0
	set id = 0
        set type              = noslip
    end
    subsection bc 1
	set id = 1
        set type              = noslip
    end

    subsection bc 2
	set id = 2
        set type              = noslip
    end
    subsection bc 4
	set id = 4
        set type              = noslip
    end
    subsection bc 5
	set id = 5
        set type              = noslip
    end

end