Thanks Wolfgang, I'll try that.
Pietro
________________________________
From: Wolfgang Bangerth <[email protected]>
To: [email protected]
Sent: Monday, April 2, 2012 8:46 PM
Subject: Re: [deal.II] How to set some mesh cells as non-conducting or non-cells
> 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
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
