On 10/21/21 7:33 AM, Мария Бронзова wrote:

I believe I answered you privately, but I cannot find my email anywhere. In case you didn"t receive it, I am sending such a test case here again. The testcase is a box porous material model with 8 cells to be able to interpolate BC's and to see the problematics clearlier. The confusing thing might be the Time function, as I am working with the frequency domain, but I simply used it artificially as a counter for frequencies to generate a series of output files.

The problem appears in the pressure boundary condition, as I am trying to force (250,0) complex-valued boundary values on certain faces. After solving the system, the Dirichlet boundary values are not held, so (250,0) becomes somewhat else.

Please, let me know, if some other information is required!

Maria -- the program is still too large. For debugging these kinds of problems, it is often useful to strip everything that is not necessary out. Set all of your coefficients to one (i.e., omit them from the code altogether), remove the loop over frequencies, remove the postprocessing steps other than graphical output, etc. I bet that you can reduce your current 724 lines to 250 this way, and it will be much easier to see what could possibly go wrong in your code then.

Best
 W.


--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/

--
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/33d2f98c-6756-f7c3-af72-f5a6469b4389%40colostate.edu.

Reply via email to