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
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 c
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 adv
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 OutputToo
Hi all,
I am solving a problem where one of the functions in the assembly process
requires the knowledge of the boundary coordinate of the mesh, which moves
at each timestep.
My solution should be only varying in the y-direction (it is indeed with
the current output) and I have a known linear
For step-20,
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 cell contributions for the boundary condition'.
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 per
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 part
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 com
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 t
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 ar
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 i
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 r
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 tha
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 (ce
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
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 obj
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 s
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 solving
19 matches
Mail list logo
