On 1/6/21 9:12 AM, Karthi wrote:

Question (1) Should I use VectorTools::project or VectorTools::interpolate for initializing the variables?
In practice, it does not matter very much, but interpolate() is much cheaper than project().


Both the initial conditions are not giving the desired steady state solution (possibly there is a bug that is not related with the initial condition). However, I want to make sure I understand how to set up initial condition generated using random values.

Have you tried to figure out *how specifically* it is wrong?


Question (2)Though, I set the initial condition only to one of the three solution variables, in the output I see all three variables initialized (with values different from zero) at time zero. I copied ‘output’ from tutorial step-21. Maybe the output is not generated at time zero or I made a mistake in the set up InitialValues?

I don't have access to your code, so I can't tell. The implementation of your vector_value() function looks correct to me. If you don't know which time step the function is called, then it's probably worth just running a single time step and looking at the output after that step.


Question (3)Since, I have a system of three variables (each with FE_Q<dim>(1) finite element space) in my function InitialValues I only defined InitialValues::vector_value and not InitialValues::value. Did I do right?

I think so. You would have gotten an error message (assuming you are working in debug mode) if you had done that wrong.

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/cafd3ca0-0514-a2be-9290-b219d89d14f6%40colostate.edu.

Reply via email to