Re: [deal.II] Geometric Conservation Law

2020-06-16 Thread Alexander Cicchino
Thank you for responding Wolfgang Bangerth. The GCL condition comes from the discretized scheme satisfying free-stream preservation. I will demonstrate this for 2D below, (can be interpreted for spectral, DG, finite difference, finite volume etc): Consider the conservation law: \frac{\partial

Re: [deal.II] Help with ParticleHandler

2020-06-16 Thread Wolfgang Bangerth
On 6/16/20 1:46 PM, Andrew Davis wrote: and I have gotten what I expect.  I have also tried attaching the particle properties using: unsigned int npart = 0; for( auto iter=particleHandler.begin(); iter!=particleHandler.end(); ++iter, ++part ) {   dealii::ArrayView

Re: [deal.II] Geometric Conservation Law

2020-06-16 Thread Wolfgang Bangerth
Alexander, I am wondering if anybody has also found that the inverse of the Jacobian from FE Values, with MappingQGeneric does not satisfy the Geometric Conservation Law (GCL), in the sense of: Kopriva, David A. "Metric identities and the discontinuous spectral element method on

Re: [deal.II] Geometric Conservation Law

2020-06-16 Thread Martin Kronbichler
Dear Alex, This has been on my list of things to implement and verify with deal.II over a range of examples for quite a while, so I'm glad you bringing the topic up. It is definitely true that our way to define Jacobians does not take those identities into account, but I believe we should add

[deal.II] Geometric Conservation Law

2020-06-16 Thread Alexander Cicchino
Hello, I am wondering if anybody has also found that the inverse of the Jacobian from FE Values, with MappingQGeneric does not satisfy the Geometric Conservation Law (GCL), in the sense of: Kopriva, David A. "Metric identities and the discontinuous spectral element method on curvilinear

[deal.II] Help with ParticleHandler

2020-06-16 Thread Andrew Davis
I have a question about how to set the quantities for dealii::Particles that are stored in a dealii::ParticleHandler. I have successfully created a dealii::ParticleHandler and randomly placed particles. I have also successfully solved a PDE, whose solution is stored in the vector solution.

[deal.II] Re: Help with ParticleHandler

2020-06-16 Thread Andrew Davis
Sorry for the second post---it might also be worth noting that unsigned int part = 0; for( auto iter=particleHandler.begin(); iter!=particleHandler.end(); ++iter, ++part ) { std::vector quantities(ncomponents); for( unsigned int i=0; iset_properties(quantities); } results in this