Re: [deal.II] Interest in a step-12-like DG tutorial but for the advection-diffusion equation?

2020-07-10 Thread Timo Heister
Hi Bruno, I introduced FEInterfaceValues because I was frustrated by the lack of DG examples and as a consequence DG users in the community and I realized that a more intuitive way to teach it was necessary. Of course we need more than 1 or 2 examples to do that. I have a lot of unfinished stuff o

Re: [deal.II] Interest in a step-12-like DG tutorial but for the advection-diffusion equation?

2020-07-10 Thread Wolfgang Bangerth
Bruno, I really enjoy the way step-12 is written, especially because of the use of the FEInterfaceValue class. In my opinion it makes low level stuff very easily accessible, yet it is relatively easy to understand. However, I am under the impression that there are no similar test for elliptic

Re: [deal.II] step-11 boundary element

2020-07-10 Thread Wolfgang Bangerth
On 7/9/20 8:50 PM, Alex wrote: Thank you. Do you have any recommendation for reference notes or books which explains the math behind mapping? I can see some inĀ Mapping< dim, spacedim > Class Template Reference. I'm not good with what FE book talks about which, but most FE books will have sect

Re: [deal.II] Interest in a step-12-like DG tutorial but for the advection-diffusion equation?

2020-07-10 Thread Bruno Blais
Dear Timo, A step-12-like version of Step-39 for the Laplace equation would be very interesting. I would be glad to help on making one of those. The same could be said about the Stokes equations. Ideally, once that is done, i would like to try and make one for stabilized Navier-Stokes for both

Re: [deal.II] Mesh-induced elastic anisotropy and distorting the quad. points as a way to palliate it

2020-07-10 Thread David F
I have a 2D system for which I create the stiffness tensor of an isotropic material, but for each finite element I create it with a different shear modulus. The shear modulus is random for each element (I use an exponential distribution, but any distribution leads to the same behavior as long as

Re: [deal.II] Mesh-induced elastic anisotropy and distorting the quad. points as a way to palliate it

2020-07-10 Thread Wolfgang Bangerth
On 7/10/20 9:15 AM, David F wrote: I have a 2D system for which I create the stiffness tensor of an isotropic material, but for each finite element I create it with a different shear modulus. The shear modulus is random for each element (I use an exponential distribution, but any distribution l

Re: [deal.II] step-7 remarks on superconvergence and quadrature

2020-07-10 Thread Wolfgang Bangerth
Praveen, https://github.com/cpraveen/fembook/blob/master/deal.II/ex04/demo.cc

Re: [deal.II] Geometric Conservation Law

2020-07-10 Thread Alexander Cicchino
Dear Martin or whom it may concern, I have solved the previous problem and am confirming that the conservative curl form has now been implemented and passes 2 complicated tests for GCL on symmetric and non-symmetric curvilinear grids for different polynomial degrees. Turns out that we did not n

Re: [deal.II] Interest in a step-12-like DG tutorial but for the advection-diffusion equation?

2020-07-10 Thread luca.heltai
If you want to take some more inspiration using mesh_loop, ScratchData, and FEInterfaceData, you could take a look at this: https://github.com/dealii/dealii/blob/master/tests/meshworker/scratch_data_08.cc where the FEInterfaceValues which is now included in ScratchData is used to solve SIPG on

Re: [deal.II] Interest in a step-12-like DG tutorial but for the advection-diffusion equation?

2020-07-10 Thread Timo Heister
> where the FEInterfaceValues which is now included in ScratchData Yes, I was going to suggest using that after step-12. > Ps: I am planning to extend FEInterfaceValues to be able to extract function > values and gradients, so that it can also be used easily with the automatic > differentiation

Re: [deal.II] Mesh-induced elastic anisotropy and distorting the quad. points as a way to palliate it

2020-07-10 Thread Sebastian Stark
Hi David, I'm neither an expert, nor do I know the literature well, but looking on your pictures, I think, the situations you are studying are geometrically anisotropic. Just plot the distribution of angles the faces of your inhomogenities make with the x-axis. For the quad-case, you'll get two d