Thanks, I understand what you meant!! I fixed it. Add up this.. can I ask your intuition for my another problem?
I am using deal-ii to analyze Non-newtonian fluids, which varies viscosity. For example, my viscosity eta is function of (u) field. I am solving this with iteration method.. (Picard) After coding, I am testing my code with method of manufacture. For squared domain, I implemented all boundary condition with my assumed solution and impose source term also to numerically generate my assumed solution however, the problem I got is.... also follow While my velocity field is qualitatively similar to manufactured solution, the pressure field is not following the exact form as seen below figure. It is very strange that......I can approach u-soltution form while my pressure does not.... In addition , any mesh refinement would not make more accurate solution now. If you were me, what would you check more to find out what is problem? <https://lh3.googleusercontent.com/-mOzFBuSfPGc/WJt7ZXct4rI/AAAAAAAAA_E/0UZ-rvkW6HMaGLnDA0uqM6cO51203vTbwCLcB/s1600/problem.png> 2017년 2월 8일 수요일 오후 1시 40분 46초 UTC-6, Wolfgang Bangerth 님의 말: > > On 02/08/2017 12:35 PM, Jaekwang Kim wrote: > > */_local_rhs(i) += fe_values.shape_value(i,q) *_/* > > > > */_ rhs_values[q] *_/* > > > > */_ fe_values.JxW(q);_/* > > > > rhs_values[q] is a Vector<double>. You need to say which component of it > you want to multiply with. (Hint: The correct component is > fe.system_to_component_index(i).first). > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > <javascript:> > 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
