On 9/5/24 01:07, 杨荣伟 wrote:
We implemented the work based on Step26 (We really appreciated
that). Now the problem we encountered is that there still exist
unexpected errors on the temperature field in Deal II solution (see
the figures (circles) enclosed in attachment ). We suspect that the
errors arise from the large magnitude of thermal conductivity (K_T in
the equations in the attachment), which is equal to 10^6, which will
render the stiffness matrix being errous. Does anyone can help us to
figure out the critical errors and make any suggestions?
At the end of the day, you will have to develop the skills to debug
these sorts of problems -- this will not be the last time you wonder
about the correctness of your program :-)
In your case, I think you are using an adaptive mesh. Make the problem
simpler: Does the solution look correct if you are using a uniformly
refined mesh? If yes, are you using a higher-order element (e.g.,
quadratic shape functions)? If so, go back to linear elements and make
sure the solution is correct. If so, you can try quadratic elements
again, but make sure that what you are visualizing is an adequate
representation of the actual solution -- see
https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#in-my-graphical-output-the-solution-appears-discontinuous-at-hanging-nodes
There are many more directions one can go into, but these are perhaps
good starting points already.
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/95b639b8-7bed-4937-8e87-a3041519510d%40colostate.edu.
