There's a piece(s) of copper wire in a brain. I'd like to set the cells immediately surrounding this wire (except for the top and bottom surfaces) including any splinters as non-conducting cells, i.e., I'd like to set the diffusion coefficient (scalar/tensor) to zero.
Is the model still well-posed in that case? > Is there
a way to do this in deal.II? I'm thinking of the FE_Nothing class but I'm not certain that's the best way to go.
You can use the material_indicator field of every cell to mark all such cells once at the beginning of the program. You'd then only assemble on cells for which cell->material_indicator() only returns certain values.
Best W. ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
