[deal.II] Interlacing update ghost and assemble

2021-09-30 Thread Praveen C
Dear all I use la::d::Vector in a continuous nodal FE code. I want to mix update ghost and assemble so that I dont wait for update ghost to finish solution.update_ghost_values_start() Assemble on all locally owned cells that do not need ghost values solution.update_ghost_values_finish() Assemble

[deal.II] Implementing a fiber reinforcement in a biphasic model

2021-09-30 Thread Matthew Rich
Hi all, I am trying to add some anisotropy to my model via fibers. These fibers would only be active in tension and have different orientations depending on spatial location in the model. My plan of attack was to simply store and extra value at the quadrature points that had a direction

[deal.II] Re: Issues related to usage of parallel::fullydistributed::triangulation

2021-09-30 Thread 'peterrum' via deal.II User Group
Just replace `dealii::parallel::fullydistributed::Triangulation` by `dealii::Triangulation` within `Check_meshPeriodicity()` . I think that should work. https://www.dealii.org/developer/doxygen/deal.II/namespaceTriangulationDescription_1_1Utilities.html#af575881c2cf233fe6f85d1a3a65a73f6 and

[deal.II] Issues related to usage of parallel::fullydistributed::triangulation

2021-09-30 Thread Aditya Pratap Singh
Hi, I am working with dealii@9.2.0 for one of my research projects and one of the applications includes reading from a large mesh file. So I was using the parallel::distributed::triangulation in my code which was working perfectly but now as the size of my mesh files has increased more, I have to

[deal.II] Solving Problems Involving Multiple Triangulations

2021-09-30 Thread ah.s...@gmail.com
Hello all, I would like to solve an elastic problem in a thick slab supported by a truss structure, as shown in the attached image. Looking at this problem, this is how I would discretize the domain - Use 1D truss elements for the struts shown in green; therefore, one needs to use a

Re: [deal.II] Dirichlet BC on block matrices

2021-09-30 Thread Мария Бронзова
So, there are two boundary integrals in the formulation and I am trying to implement them for the case of fixed displacement boundary condition. The first boundary integral falls to zero for such a BC, as no displacement variation is possible in this case. So, I am implementing the second

[deal.II] Apply boundary conditions using the method of lines

2021-09-30 Thread bob bill
Hello, I'm trying to solve the 2D heat equation with the method of lines, instead of Rothe's method as explained in the tutorials, as I'd like to use ARKODE to march in time. Assuming that the *boundary data is 0* and the initial data is *x*y*(1-x)*(1-y)*, after writing the weak formulation, I

Re: [deal.II] Dirichlet BC on block matrices

2021-09-30 Thread Wolfgang Bangerth
On 9/30/21 2:48 AM, Мария Бронзова wrote: Thank you a lot for the reply! I am trying to implement the Biot equations for the propagation of acoustic waves in porous media. The equations are attached in the "Biot .PNG" file. There are two boundary integrals there, containing information on

Re: [deal.II] Small suggestion to improve GridIn::read_unv()

2021-09-30 Thread Wolfgang Bangerth
On 9/30/21 6:02 AM, vachanpo...@gmail.com wrote: I found this thread which highlights one small issue in importing a unv file (generated by Salome): https://groups.google.com/g/dealii/c/iWnLv4v0V1M/m/38jm5Ck8WXYJ. If I understand correctly, dealii expects the file to start with "-1" and

[deal.II] Small suggestion to improve GridIn::read_unv()

2021-09-30 Thread vachanpo...@gmail.com
Dear all, I found this thread which highlights one small issue in importing a unv file (generated by Salome): https://groups.google.com/g/dealii/c/iWnLv4v0V1M/m/38jm5Ck8WXYJ. If I understand correctly, dealii expects the file to start with "-1" and then have the section "2411", instead of all

Re: [deal.II] Dirichlet BC on block matrices

2021-09-30 Thread Мария Бронзова
Dear Wolfgang, Thank you a lot for the reply! I am trying to implement the Biot equations for the propagation of acoustic waves in porous media. The equations are attached in the "Biot .PNG" file. There are two boundary integrals there, containing information on work, done by external forces

Re: [deal.II] Re: function face_component_to_system_index() for Discontinuous Galerkin ?

2021-09-30 Thread Sylvain Mathonnière
*About writing the function :* I must confess I did not try to write the function myself but I did check to see how it was implemented and it seems like it is just retrieving the value from a table. I currently only have the user version of deal.II (v 9.2.0), If I wanted to write this