Yes exactly. The only difference is in weak form of my equation which should be multiplied by 2*pi*r which I define r_value = fe_values_eta.quadrature_point(q_point)[0]; In fact, I should solve my problem on a sphere (circle). According to your noteworthy comments, the only thing that I have to do is change the geometry and use higher order mapping for accuracy. Another thing is that consider a half of a circle in r-z coordinate system and the problem is axisymmetric. what should be my boundary condition for rotation axis (z-coordinate)? For the r-coordinate, I can use no-normal flux boundary conditions or just limit the displacement normal to the surface, however, for rotation axis I am really confused.
On Wednesday, October 19, 2016 at 1:50:49 PM UTC-5, Timo Heister wrote: > > > Yes You are right. I mean curved domain. My question is should i use > codim > > instead of mapping because of the fact that there is only two tutorials > in > > dealii which solved a differential equation on c curved domain according > to > > No, codim problems solve systems on surfaces embedded in a higher > dimensional space. This has nothing to do with solving on a domain > that is not a rectangle. There are other examples that use a different > domain, see step-32 for example. > > > the fact that my derived weak expressions are in cylindrical coordinate > > systems? > > Are you sure you want to describe your problem in cylindrical > coordinates? While possible, this probably only makes sense if you > have additional symmetries in your system that you want to exploit, > but you will have many other issues with that (quadrature etc.). > > If you just want to solve on a disk or sphere, you just need to change > the geometry and potentially use a higher order mapping for accuracy. > > > > > Thanks, > > Benhour > > > > > > On Wednesday, October 19, 2016 at 1:04:04 PM UTC-5, Timo Heister wrote: > >> > >> > Let me put it in a such way. I want to solve my equation on a curved > >> > edge. > >> > According to dealii library, I should use MappingQ whenever I have > >> > nonlinear > >> > domain. Am I correct? > >> > >> I don't think "nonlinear" is the word you are looking for. Do you mean > >> a curved boundary? > >> > >> > In addition, Should I replace all dim in defining my > >> > tensors to spacedim? It should be noted that dim = spacedim - 1. > >> > >> Using a mapping (see step-10, step-11, etc.) has nothing to do with > >> codim 1 or 2 problems (where dim is not equal to spacedim, see > >> step-38). > > > -- > Timo Heister > http://www.math.clemson.edu/~heister/ > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
