Re: [deal.II] what is the different between VectorTools::point_value and fe_values.get_function_values()?

On 3/13/19 8:21 PM, Phạm Ngọc Kiên wrote: > I am testing my codes for output the solution at a point in my numerical > model. > I saw in the library 2 ways to do this task. The first one is to use > VectorTools::point_value() function and the second one is > fe_values.get_function_values(). The

[deal.II] what is the different between VectorTools::point_value and fe_values.get_function_values()?

Hi colleagues, I am testing my codes for output the solution at a point in my numerical model. I saw in the library 2 ways to do this task. The first one is to use VectorTools::point_value() function and the second one is fe_values.get_function_values(). For the latter function, I am using the c

[deal.II] ILU preconditioner and mesh-scale dependency?

Hi I had a trouble in solving scalar transport equation using DG elements. For example, I was solving for \lambda in the following equation where *u *is also a numerical solution [image: Screen Shot 2019-03-13 at 9.08.17 PM.png] The thing is I do not get any problem when my mesh file is sc

Re: [deal.II] Meeting some questions using constraints to apply boundary values

> I am coding for my graduation design, and hoping to achieve a coupling > geomechanics and fluid model. > I use FE_RT for velocities and FE_Q for displacement of the rock. > > In order to valid my model and programming, I use 1 D Terzaghi. > But my boundary conditions appears have some problem

Re: [deal.II] Re: About convergence using DG

On 3/13/19 9:18 AM, Apurva Tiwari wrote: > > So far I'd been using an RK3 integrator that I had written myself.  To > rectify > the problem I faced (mentioned in the mail earlier),  I thought of, was using > the built-in explicit RK4 integrator TimeStepping::ExplicitRungeKutta >

Re: [deal.II] About convergence using DG

On 3/12/19 5:34 AM, Apurva Tiwari wrote: > I'm trying to get appropriate convergence rates for a 2D scalar advection > problem. The strategy is quite standard. I start with a coarse grid, > successively refine it globally and generate the convergence table as done in > step-7. As can be seen in

Re: [deal.II] How to apply first kind of boundary conditions to one component of a vector

On 3/10/19 7:39 AM, chenxi.edu1...@gmail.com wrote: > > In order to check if this code is right, I want to use a model with > analytical > solution to verify it. And I need to make the displacement in y direction to > be zero. > At first I use fe.component_mask(displacement), > But this works f

[deal.II] Re: About convergence using DG

Hello! So far I'd been using an RK3 integrator that I had written myself. To rectify the problem I faced (mentioned in the mail earlier), I thought of, was using the built-in explicit RK4 integrator TimeStepping::ExplicitRungeKutta

Re: [deal.II] issues with manifold id (at least in codim-1) revealed by PR 7775

I believe the first point is a bug of the read_ucd function. I opened an issue for this (https://github.com/dealii/dealii/issues/7802) The second point, however, is a feature. What should you do when two material ids meet? I think it is the responsability of the user to make sure that the corr