Dear Wolfgang, Thank you for your response. I am working on phase-stability in materials, for smooth initial conditions it produces the correct phase but for random initialization it produces a phase which I am not expecting. The phases are every sensitive to the initial-conditions used. If I happen to figure out the issue, I will post back my findings.
Cheers, Karthi. On Friday, January 8, 2021 at 5:48:01 AM UTC+1 Wolfgang Bangerth wrote: > 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/037afdb3-824c-4cf1-b9fa-4c1969c1b451n%40googlegroups.com.
