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

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

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

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

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

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.

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