Hi Felix,
> Okay. When I read you e-mail I sincerely hoped it wasn't the problem. But it
> was. Now the results match.
That’s great that this was an easy fix! I’m glad that you can now move along
with your problem :-)
> I feel so so stupid. Such a rookie mistake. I already made this mistake 3
Okay. When I read you e-mail I sincerely hoped it wasn't the problem. But
it was. Now the results match. I feel so so stupid. Such a rookie mistake.
I already made this mistake 3 years ago and forgot about it.
Thank you very much for taking the time to look through this part of the
code. I own
Felix,
Jean-Paul's answer might already have solved your problem. Anyway:
I'm not sure to understand how to manufacture a solution for my problem. I
> have read step 7 a few times and I understand but for a coupled system like
> mine, I don't see how to do it.
> Should I choose a solution with u
Dear Felix,
I don’t claim to have read through your comprehensive documentation of your
issue in any detail, but one observation of your code is that you have a couple
of suspect terms in your local matrix assembly. In particular, the terms of
this form
> 1/2/phiscale*(u_j_y-v_j_x)
do not compu
One final remark for today, if I remove the two asymmetrical terms in my
equation so it becomes :
*Then all three codes gives back the same solution...*
So I thought it would come from my assembly of this part of the system but
I have checked and checked again... It should not come from my wea
>
> If the reference solution is 0, you can't see if you are missing a
> constant factor in your equation even when the rate of convergence looks
> correct.
>
I understand thank you !
> Yes, but then you also have to be careful with the boundary conditions. If
> you have an analytical solut
Felix,
> Is the rate of convergence as expected?
>>
>
> I believe so. I did a linear regression of the convergence graphs in
> logscale and got a value of :
> -7.6 for the speed : I believe it corresponds to the -8 from second order
> finete element in 2D
> -3.2 for the pressure : I beliebe it co
> I have tried a simple testcase for my problem.
>> I don't know if it's suitable as a testcase but if I use X=0 in my
>> equations above, I know the solution.
>>
>> The solutions is u = v = 0, \phi = (1-x^2 -y^2) and the pressure is then
>> P = 8 \phi_lin^2 (1-x²-y²)
>> It fits with my bound
Felix,
in addition to what has been suggested before, think also about the fact that
you don't know which of the two solutions you have is correct. Before you draw
any conclusions, you need to find a way to establish which one is the one that
is correct and which one is not.
Best
W.
On 7/