Thank you very much for your answer.
So for the first function
VectorTools::point_value(dof_handler,solution,point,values);
I can specify the observed point in real coordiante, solution vector and
dof_handler to get the values right?
For the second one, assuming that I am using QGauss for
Thank you very much for your answer.
So for the first function
VectorTools::point_value(dof_handler,solution,point,values);
I can specify the observed point in real coordiante, solution vector and
dof_handler to get the values, right?
For the second one, assuming that I am using QGauss for
Dear Pham,
When you use the FEValuesExtractors in conjunction with the FEValues class
using the operator[] like you’ve done here
> fe_values[vector_re]
then what is returned is an FEValuesViews type object. Since your vector_re
extractor is of type FEValuesExtractors::vector, the returned
Hi Wolfgang, Sorry for the reply - for some reason i didn't get the
notification.
What i mean by correct is that I do have two uncoupled equations, but that
it is correct in that I had already verified it using what you suggested
(by using exact solutions - MMS).
By it fails, I mean that the
Note that I think this was the issue.
You would need to put local_rhs = 0 within the part below to apply the
Dirichlet condition. Please do let me know if you think this is wrong but I
am getting the correct output. Thanks
for (unsigned int face_no=0;
face_no::faces_per_cell;
++face_no)
if
On 3/14/19 1:28 AM, Phạm Ngọc Kiên wrote:
>
> So for the first function
> VectorTools::point_value(dof_handler,solution,point,values);
> I can specify the observed point in real coordiante, solution vector and
> dof_handler to get the values, right?
Correct.
> For the second one, assuming
On 3/13/19 8:18 PM, Jaekwang Kim wrote:
>
>
> I wonder now ...
>
> 1. why smaller scale mesh only results such error,
You have two terms in your matrix, one from the advection term and one
from the reaction term. If you scale the mesh by a factor of c, the
height of the shape functions used
On 3/14/19 9:13 AM, jane@jandj-ltd.com wrote:
> Note that I think this was the issue.
>
> You would need to put local_rhs = 0 within the part below to apply the
> Dirichlet condition. Please do let me know if you think this is wrong
> but I am getting the correct output. Thanks
That seems
Evgeny,
> I found one problem with inp file - in line "1 0 quad 1 2 3 4". I
> guess the second number should be nonzero here. I have try to set 1 and
> get better results. Now I can see velocity changes around the
> input-output bounds.
Is this an input file you created yourself, or
Hi all,
as an addition to this post, For Dirichlet condition that are in the weak
form, would this have to be done?
And for example, for an inhomogeneous normal component of the normal stress
condition in step-22 for Stokes (1st of the partial bc discussion in the
tutorial), do you also have
Hi Wolfgang,
I would think that it might be useful in the tutorial!!
Please do let me know if you have any suggestions on a perhaps analogous
question for step-22 in my other reply to this thread. I'm making the
analogy between the dirichlet boundary condition implementation and the
normal
Hi Jane,
> I would think that it might be useful in the tutorial!!
The same applies here: I deal with too much email to be clear about what
"it" replies to :-) If you think "it" would be useful, can you make a
concrete suggestion?
> Please do let me know if you have any suggestions on a
On 3/14/19 11:42 AM, jane@jandj-ltd.com wrote:
>
> as an addition to this post, For Dirichlet condition that are in the
> weak form, would this have to be done?
>
> And for example, for an inhomogeneous normal component of the normal
> stress condition in step-22 for Stokes (1st of the
I highly appreciate all of your answers. It becomes clearly for me to look
through the library.
Let me show all of my function and my ideas below.
I set the point for evaluating its solution to be on the vertex of a cell,
and using the following function to get the solution
template
void
On 3/14/19 5:02 PM, jane@jandj-ltd.com wrote:
>
> I think just a comment in the bit on how to implement the dirichlet bc in the
> weak form would be sufficient - something to say 'In the case of an
> inhomogeneous boundary condition, you would need to set local_rhs = 0 before
> adding the
Thanks for sharing insight!
On Thursday, March 14, 2019 at 12:18:30 PM UTC-5, Wolfgang Bangerth wrote:
>
> On 3/13/19 8:18 PM, Jaekwang Kim wrote:
> >
> >
> > I wonder now ...
> >
> > 1. why smaller scale mesh only results such error,
>
> You have two terms in your matrix, one from the
On 3/14/19 5:50 PM, jane@jandj-ltd.com wrote:
>
> I've tried looking a cell->face(face_no)->at_boundary() type input within the
> cell iterator and by using something like cell->vertex(v)(1) but I can't
> quite
> get my head around how to find the y coordinate of the top boundary for every
17 matches
Mail list logo
