Thank you for your reply.
My problem seems to be that constraints not only constraint the boundaries
but also constrain the interior.
Maybe I should take one field at a time to debug.
Best
Y.
在2024年6月11日星期二 UTC+8 06:56:36<Wolfgang Bangerth> 写道：
>
>
> On 6/1/24 21:24, Tiny Y wrote:
> > The solution for one of my physical field models is not changing in any
> > way. I made a 2D model before and was able to calculate the results, but
> > when I switched to a 3D model I had the above problem, my mesh was
> > generated with
> > GridGenerator::cylinder_shell(triangulation,1000,r1,r2,160,10), the
> > boundary conditions were applied at r1=0.05 with a temperature of -53°C
> > The boundary condition imposed at r2=0.45 is 0°C, and the initial
> > condition is 0°C. However the temperature is constant at all moments
> > except time=0 when it is zero。图片1.jpg
>
> Tiny Y:
> I don't really know what it is that is going wrong without seeing the
> whole code, but I think that what you are doing in your function that
> sets boundary ids is to ask whether for a given face center (x,y) you have
> x^2 + y^2 < 10^{-4}*r
> Are you sure that there are face centers that satisfy this criterion?
> Similarly for the other conditions in the following else if blocks. This
> is where I would start to investigate.
>
> Separately, you could write that condition much simpler as
> if (cell->face(face_number)->center().norm() < 1e-4 * r)
>
> Best
> W.
>
--
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/2d529953-bd61-496b-957f-d31b496f5f88n%40googlegroups.com.
