Hi Xujun,
I am working together with Wolfgang to implement XFEM in deal.ii.
Soon there will be a kind of tutorial available from the website of deal.
You need a special contruction to manage the enriched dofs. In the
tutorial we will explain how to use them for different kind of
discontinuities.
Which kind of application do you have in mind?
Thomas
> Hi there,
>
> Is there anyone working on xfem? I browsed the online tutorial, but
> unfortunately couldn't find any implementation of XFEM in deal II. I have
> written XFEM, Fast marching, level set codes with Matlab, but not familiar
> with deal ii library yet. anyone can give me some hints? I read some
> examples recently, but still have some points unclear
>
> 1. By justifying the signs of level set function in each element, I can
> easily identify the standard element and enriched element. But how can I
> change the settings of dof's in these elements. For example, in the
> enriched elements, additional dof's exist due to the enrichment. Also, the
> local_dof_indices are changed in this process.
>
> 2. We can achieve loops over each active cell to calculate local stiffness
> matrix using DoFHandler<dim>:: active_cell_iterator. My question is can I
> do the loop operation over each nodes and also find the corresponding
> supporting elements around it? Any connection functions from nodes to
> elements? for example, 3x3 Q1 element, totally 4x4 nodes, I would like to
> do the following loops
> for i=1:total_num_nodes
> support_element_number = some_function( i );
> n_support_element = size( support_element_number );
> for j=1:n_support_element
> ... operations in current j-th element
> ...
> end
> end
>
>
> 3. If I assembled two sets of global matrix, say, [ A ], [ B ] and rhs
> {f},
> respectively
> [ A ]*{u} + [ B ]*{a} = {f}
> [ B ]_T *{a} = {g}
> how can I combine them as a big matrix K ? for example,
> K = [ A B ]
> [ BT 0 ]
> F = {f}
> {g}
>
> Thank you very much.
> _______________________________________________
> dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
>
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
