On 02/09/2017 04:45 AM, Jaekwang Kim wrote:
The manufactured solution of velocities has not satisfies the continuity
equation.
(i.e. manufactured solution does not satisfy, div.u=0)...
Yes, that would do it :-)
Best
W.
--
Dr. Bangerth
Thank you,
I just fixed up what was wrong...
in most cases, it is usually easy problems.
The manufactured solution of velocities has not satisfies the continuity
equation.
(i.e. manufactured solution does not satisfy, div.u=0)...
Jaekwang Kim
--
The deal.II project is
Thank you for your advice.
I would make the problem simpler. Do you get the right solution if you
> have a constant viscosity? If you iterate, do you get the solution after
> one iteration?
>
Yes, I checked this.
After one iteration, I get the solution after one iteration when constant
On 02/08/2017 01:15 PM, Jaekwang Kim wrote:
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
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..
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. You need to say which component of it
you
Hi all,
I was testing my code, with method of manufactured solution.
I have a source term, which is vector value, defined as ..
template
class RightHandSide : public Function
{
public:
RightHandSide () : Function(dim+1) {}
virtual void vector_value (const Point ,