How to apply a point source condition (point boundary condition) in Discontinuous Galerkin methods

[Jiaqi Zhang]

Dear All,

I was trying to use DG to solve a Hamilton-Jacobi equation, and the problem has a point source condition

for example, the center of a unit square pointCenter = 0.

First I thought of using

ZeroSourcePoint.add_line(cell->vertex_dof_index(v,0));

but later found that DG elements aren't defined at the vertices of the mesh in FAQ.

So how do I impose this point condition in DEAL II?

Thanks in advance!

Best,

Jiaqi

[Bruno Turcksin]

Jiaqi,

to create a point source, you can use VectorTools::create_point_source_vector (http://dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html#ad03b858b1a3b59003a76f6224e67efc7)

Best,

Bruno

[Jiaqi ZHANG]

Thanks Bruno, I should have searched a little bit more.

Best,

Jiaqi

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

--

Jiaqi Zhang

Graduate student

Department of Mathematics

Virginia Tech